r 

n 

AD-A128  643  CCD  MATRIX  PROCESSORS)  HONEYWELL  SYSTEMS  AND 
CENTER  MINNEAPOLIS  MN  P  C  ROBERTS  ET  AL  API 
82SRC59  RADC-TR-82-294  F 1 9628-80-C-0 1 54 

UNCLASSIFIED 

RESEARCH 

83 

/G  9/2  NL 

■ 

* 

□ 

■ 

■ 

■ 

■ 

■ 

■ 

■ 

■ 

■ 

! _ 

■ 

II 

■■ 

■■ 

li 

1 

■■ 

■■ 

1 

■■ 

■■ 

■ 

■ 

1 _ _ 

1 


I 


* 


xs 


l 


f 

6 

1 

I 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BURiAU  01  STANDARDS  1963  A 


V  *  a 


tV-'.-'Wjsrf  * 


UNCLASS IFIED 


SECURITY  CLASSIFICATION  OF  THIS  P  AGE  (Whit  D«atnlill« 


REPORT  DOCUMENTATION  PAGE 


READ  WSTRUCTtONS 
BEFORE  COMPLETING  FORM 


t.  report  number 
RADC-TR-82-294 


J.  OO VT  ACCESSION  NO 


1  RECIPIENT1*  CATALOG  NUMBER 

n _ 


4.  TITLE  (and  Submit) 


»■  type  of  report  a  period  COVEREO 


CCD  MATRIX  PROCESSOR 


Final  Technical  Report 
25  Aug  80-1  Mar  82 

*  PERFORMING  OTG.  REPORT  NUMBER 


T.  AUTHORfO 


» 


82SRC59 _ 

.  CONTRACT  OR  GRANT  NUMBER/^ 


P.  C.  T.  Roberts  R.  Fitch 

J.  D.  Joseph  B.  R.  Hanzal 

J.  A.  Hoschette  D.  B.  Yanke 

».  PERFORMING  ORGANIZATION  NAME  ANO  ADDRESS 

Honeywell  Systems  and  Research  Center 

2600  Ridgway  Parkway 

Minneapolis  MN  55413 _ 

H.  CONTROLLING  OFFICE  NAME  ANO  AODRESS 


F19628-80-C-0154 


10  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  A  WORK  UNIT  NUMBERS 

61102F 

2305J135 

II.  REPORT  OATE 


Rome  Air  Development  Center  (ESED) 

Hanscom  AFB  MA  01731 

14.  MONITORING  AGENCY  NAME  A  AOORESVlf  diTfSUnt  horn  Controlling  Otllca) 


April  1983 

IS.  NUMBER  OF  PAGES 

234 _ 

IS.  SECURITY  CLASS,  (o I  ihl»  r«Rort) 


Same 


14.  DISTRIBUTION  STATEMENT  la  I  tbit  R»porlj 


UNCLASSIFIED 

15«.  DECLASSIFICATION' DOWN  GRADING 
. SCHEDULE 

N/A _ 


Approved  for  public  release;  distribution  unlimited. 

I?.  DISTRIBUTION  STATEMENT  (oJ  the  obotropt  entered  In  Block  20,  It  different  from  Report) 

Same 


IB.  SUPPLEMENTARY  NOTES 

RADC  Project  Engineer :  Paul  W.  Pellegrini  (ESED) 


It-  KEY  WOHOl  (Continue  on  roooeee  aide  It  nocooooey  end  tdorttfy  by  block  number) 

Parallel  image  processor 
CCD  matrix  array 
Image  processing  algorithms 
Threshold,  edge  enhancement 

20  ABSTRACT  (Continue  on  roooroo  aide  It  mcmmtt  m4  I  don  tiff  by  block  number) 

>A  new  CCD  Matrix  Array  has  been  fabricated  including  a  Microprocessor 
Host-CPU  and  Controller  Exerciser  Electronics  Unit.  Each  unit  cell  of 
the  array  performs  analog  arithmetic  and  logic  operations  on  input  data 
matrix  elements  in  the  Single  Instruction  Multiple  Data  architecture. 
Ultra-high  signal  and  image  processing  throughput  is  demonstrated. 


00  ,  1473  EOlTIONOF  I  NOV  A*  I*  OBSOLETE  UNCLASSIFIED  \ 

1ECURITV  CLASSIFICATION  OF  TNI*  PAOE  fAto*  Ama  M«W) 


I 

H 


CONTENTS 


Section  Page 

I  INTRODUCTION  1 

Overview  of  the  Contract  1 

Synopsis  of  the  Report  5 

List  of  Items  Delivered  to  RADC  5 

Reporting  5 

Hardware  5 

Documentation  6 

II  DETAIL  OF  TASKS  PERFORMED  7 

Task  I-- Chip  Exerciser  Design  and  Fabrication  7 

I/O  Board  Addressing  7 

Timing  Generator  Board  9 

Timing  Generator  Multibus  Commands  9 

Timing  Generator  2K  x  48  Bit  Memory  17 

Double-Pipeline  Branching  20 

Timing  Generator  Memory  Map  21 

Condition  Codes  (CCj_o)  21 

I/O  Buffer  Board  21 

I/O  Buffer  Board  Multibus  Commands  22 

Task  2— Initial  CCD  Characterization  27 

Original  PIP  Array  and  Cell  Design  27 

Initial  Measurements  32 

Task  3 — Design  Iteration  and  Fabrication  of  16  x  16  Array  56 


iii 


CONTENTS  (continued) 


Section 


III 


Page 


Design  Modifications  56 

Improved  Functionality  65 

Task  4 — Final  16  x  16  CCD  Characterization  72 

Chip  Tests  and  Wafer  Tests  of  #3022  CCDs  72 

System  Tests  with  Original  #2214  CCDs  89 

Task  5->-System  Architecture  and  Software  Definition  94 

Target  Acquisition  and  Tracking  Algorithms  95 

Processor  Architectures  for  Signal  and  Image 
Processing  99 

Matrix  Processor  Control  102 

PIP  Architecture  104 

A  Matrix  Processor  Example  107 

Summary  of  Example  118 

PIP  Modeling  122 

Prepared  Architectural  Improvements  125 

Task  6 — Initial  Design  of  32  x  32  Array  Cell  130 

Task  7 — Design  Optimization  Review  With  RADC  137 

PROJECTIONS  OF  FUTURE  DEVELOPMENTS  OF  MATRIX  PROCESSOR 

CCDs  WITH  PARALLEL  I/O  139 

Sensor/Processor  139 

Software  and  Control  140 


iv 


CONTENTS  (concluded) 


Section  Page 

IV  SUMMARY  AND  CONCLUSIONS  143 

APPENDIX  A.  HP41C  POCKET  COMPUTER  LISTING  145 

APPENDIX  B.  HP9845C  DESKTOP  COMPUTER  BASIC  PROGRAM  157 

APPENDIX  C.  SOFTWARE  USER'S  GUIDE  177 

REFERENCES 


Accession  For 

NTIS  GRA&I 

DTIC  TAB  ' 

Unannounced 

Just J fication 

E 

□ 

Br 

E 1  "  ;  ion./ 

Avr> -3  '  ---  Cod 

03 

i/or 


* .  -  t  1 

w->  % 


Figure 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 


LIST  OF  ILLUSTRATIONS 

Page 


Completed  Unit  Together  with  Associated  Cameras,  vtr, 

and  Monitor  8 

Block  Diagram  of  Timing  Generator  Board  11 

PIP  Clock  Generation  Format  18 

Branch  Address  Format  18 

I/O  Buffer  Board  23 

Photograph  of  original  #2214  CCD  Chip  28 

Calma  Plot  of  the  2214  Cell,  Together  with  a  Photograph 
of  the  Actual  Cell  (Poly  I  is  shaded)  29 

Regenerator  Operation  31 

Clock  Waveform  Timing  Diagram  for  One  Primitive  (Phases 

not  shown  are  at  rest  in  initial  state.)  33 

Leakage  Residual  Buildup  at  Frame  Period  of  23  msec, 

25°C  (2214)  34 

Input  Fill  and  Spill  to  Sense  7-amp  Performance  Using 
Serial  I/O  Registers  Only  35 

XY  Display  Using  the  Routine:  El  E2  E3  LI  R1  L2  Rl  L2  U2 
SU  SU  SU  SL  SL  SL  L2  M2  35 

Video  Waveform  at  Top  Line  of  6  x  6  Test  Pattern  (Zero 
reference  at  center  of  screen.)  36 

Typical  Video  Output  After  Buffer  38 

Zero  Data  Leakage  Measurement:  Shift  Right  (2214)  40 

Zero  Data  Leakage  Measurement:  Shift  Up  (2214)  41 

Integration  Mode  Leakage  of  STl:  Room  Temperature  42 

Histogram  of  Data  in  Figure  17  43 


i 

I 


vii 


LIST  OF  ILLUSTRATIONS  (continued) 


Figure  Page 

19  Integration  Mode  Leakage  of  ST1:  Cooled  ~0°C  44 

20  Histogram  of  Data  in  Figure  19  45 

21  Integration  Mode  Leakage  of  ST2 :  Room  Temperature  46 

22  Histogram  of  Data  in  Figure  21  47 

23  Integration  Mode  Leakage  of  ST2:  Cooled  ~0°C  48 

24  Integration  Mode  Leakage  of  ST3:  Room  Temperature  49 

25  Histogram  of  Data  in  Figure  24  50 

26  Integration  Mode  Leakage  of  ST3i  Room  Temperature  51 

27  Histogram  of  Data  in  Figure  26  52 

28  Transfer  Inefficiency  Effects  in  I/O  Registers  54 

29  Magnified  Scale  of  Data  in  Figure  28  55 

30  Transfer  Inefficiency  Effects  in  I/O  and  Array  Registers  57 

31  Magnified  Scale  of  Data  in  Figure  30  58 

32  Modified  Regenerator  Gate  Structure  Showing  Reduced 

Loading,  Increased  Sensitivity,  and  C-Clock  Screening  60 

33  CALMA  Plot  of  3022  Cell  61 

34  CALMA  Plot  of  2214  Cell  62 

35  CALMA  Plot  of  Metal  Tracks  on  3022  Cell  63 

36  CALMA  Plot  of  Metal  Tracks  on  2214  Cell  64 

37  Revised  40-Pin  DIP  Pin  Assignments  for  *3022  Chip  66 

38  Shift  Right  (SR)  Clock  Sequence  with  Termination  Phase  T 

as  Now  Implemented  67 


viii 


LIST  OP  ILLUSTRATIONS  (continued) 


Figure 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 


Page 


Shift  Up  (SU)  Clock  Sequence  with  Termination  Phase  T 

as  Now  Implemented  68 

Sequence  of  Primitives  for  #3022  to  Generate  Roberts  Cross, 

R,  and  Save  Original  Image,  a,  for  Further  Processing  71 

Photograph  of  Entire  3022  Chip  73 

Leakage  Current  with  Zero  Data  (25%,  47  msec)  7  5 

Leakage  Current  with  Zero  Data,  Without  Erasing  Stores  75 

Performance  with  Data  (70%  full  well,  0%  fat  zero) 

Showing  Transfer  Efficiency  77 

Performance  with  Data,  Without  Erasing  Stores  77 

Loading  and  Unloading  ST1  78 

Loading  and  Unloading  ST2  78 

Loading  and  Unloading  ST3  79 

Swapping  Data  (ST2,  ST3) ;  Software  Error  Suspected  80 

Leakage  from  ST1  (25°C,  47  msec)  81 

Leakage  from  ST2  (25°C,  47  msec)  81 

Leakage  from  ST3  (25°C,  47  msec)  82 

Regenerator  Performance — Output  Charge  84 

Regenerator  Performance — Input  Charge  84 

Regenerator  Performance  Showing  Linearity  85 

Physical  Basis  for  Multiplier  Gain  Control  (Patent 

applied  for:  Honeywell  #a4109674-US)  86 

Multiplier  Performance  (MG2  High)  88 

Multiplier  Performance  (MG2  Low)  88 


ix 


Figure 


LIST  OF  ILLUSTRATIONS  (concluded) 


SEM  Photograph  of  3022  Cell 

ST3  Metal  Bue  and  Contact  to  Cross-Over  Screen  on  3022 
Chip 

RSNI  to  Poly  I  Contact  on  3022  Chip 
FLIR  Image  To  Be  Processed  by  2214  Chip 

Processed  Images  Using  Threshold  and  Subtract  Subroutines 

Functional  Block  Diagram  of  Generic  Automatic  Target 
Screening/Cueing  Approaches 

Five  Techniques  for  Computing  Two-Dimensional  Filters  in 
PIP 

PIP  Cell  Model 

Ring  Store  Architecture 

Clock  Phase  Two-Metal  Interconnect  CALMA  Plot  of 
Second-Generation  Matrix  Processor 

Field  Cut  and  Two-Poly  Layer  CAUiA  Plot  of  Second- 
Generation  Matrix  Processor  (Notes  Racetrack  internal 
to  cell) 

Composite  CALMA  Plot  of  Second-Generation  Matrix 
Processor 


LIST  OF  TABLES 


Table  Page 

1  Multibus  I/O  Port  Addresses  10 

2  Clock  Assignments  19 

3  Condition  Codes  22 

4  CCD  Matrix  Processor  Primitives  (Chip  #3022)  :  Level  1 

Commands  69 

5  2214  Matrix  Processor  Mnemonics  for  Function 

Implementation  104 

6  PIP  Functions  Needed  for  Implementing  Tracker  and 

Acquisition  Algorithms  108 

7  Processing  Sequence/Missile  Seeker  109 

8  Nonseparable  Seeker  Algorithms,  Template/Centroid/ 

Tracking  Implementation  116 

9  Separable  Seeker  Algorithms,  Template/Centroid/ 

Tracking  Implementation  119 

10  PIP  Execution  Times  for  Missile  Seeker  Algorithms  (2214)  121 

11  Model  of  the  PIP  Intracell  Instruction  Set  126 


12  Execution  Times  for  Seeker  Algorithms  (New  Cell  Functions)  130 


xi 


SECTION  I 


INTRODUCTION 


OVERVIEW  OF  THE  CONTRACT 


The  objectives  of  the  work  described  in  this  final  report  on  contract 
F19628-80-C-0154  are  straightforward: 

o  Analyze  by  simulation  and  measurement  the  performance  capabilities 
of  an  ultra-high  throughput  CCD  matrix  processor  employing  parallel 
analog  processor  cells 

o  Design  a  second-generation  CCD  matrix  processor  cell  capable  of 
executing  a  complete  set  of  image  processing  algorithms 

o  Construct  a  flexible,  programmable  Exerciser  unit  capable  of 
operating  the  CCD  matrix  processor  chip  in  all  modes  of 
functionality 

o  Deliver  the  15  best-effort  CCD  matrix  procrssor  chips  and  the 
Exerciser  unit,  complete  with  basic  software 

All  of  the  objectives  have  been  met.  The  contract  was  extended  as  per 
modification  P00003  requested  by  RADC  to  provide  extra  time  and  funding  for 
the  greatly  improved  capabilities  desired  for  the  Exerciser  unit.  The  unit 
now  contains  expansion  capabilities  to  meet  future  matrix  processor 
requirements  without  significant  hardware  modifications.  The  Exerciser  is  a 
microprocessor-controlled  unit  capable  of  exercising  not  only  parallel  image 
processing  (PIP)  algorithms,  but  also  a  broau  range  of  general  signal 
processing  functions  which  are  expressible  as  single  instruction,  multiple 
data  (SIMD)  sequences  of  primitive  operators. 
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Early  characterization  tasks  revealed  that  the  original  CCD  cell  and  array 
designs  were  deficient  in  some  key  respects.  After  consultation  with  radc 
(Contract  Monitor,  Mr.  Paul  Pellegrini),  it  was  decided  that  suitable 
redesign  of  the  mask  set  was  required  and  would  greatly  improve  the 
performance  of  the  final  CCDs  to  be  delivered  to  RADC. 

The  changes  involved  provision  of  charge  sink  diodes  at  the  originally 
unterminated  boundaries  of  the  matrix  and  also  a  redesign  of  the  structure 
of  the  floating  gate  amplifier  (PGA)  and  gain-control  circuitry 
(Multiplier) .  Funding  for  the  new  mask  set  was  provided  by  internal 
Honeywell  sources. 

No  architectural  changes  were  made  to  the  cell  design  at  that  time  since  a 
minimum-risk  improvement  was  desired:  functional  CCDs  from  the  original 
mask  set  had  been  obtained  before  the  start  of  the  contract.  However,  low 
yield  of  chips  has  continued  to  hamper  characterization  of  some  aspects  of 
the  cell  performance.  A  remaining  design  error  discovered  after  the  first 
wafer  processing  run  was  immediately  corrected.  However,  the  second  and 
third  wafer  runs  (12  wafers  each)  proved  to  be  substandard  due  to  high 
charge  transfer  inefficiency  and  also  to  random  short-circuit  metal  tracks. 
Additional  wafer  process  runs  have  been  performed  on  Honeywell  funding.  The 
fourth  run  has  yielded  excellent  transfer  inefficiency  data  and  verified 
that  the  design  modifications  are  effective  in  completely  removing  the 
termination  blemishes.  We  have  now  obtained  some  functional  devices,  but 
because  of  randomly  open  circuit  control  gates  in  the  cell  and  threshold 
nonuniformity,  there  are  no  perfect  chips  available.  There  are  no  known 
design  errors  and  failure  analysis  in  progress  to  identify  the  causes  of 
the  fabrication  problems.  Other  types  of  CCDs  fabricated  on  the  same 
facility  have  yielded  satisfactorily  good  area  arrays,  and  there  is  no 
obvious  explanation  for  the  random  defects.  Wafer  runs  and  probing  will 
continue  on  Honeywell  funding  until  a  satisfactory  result  is  obtained,  since 
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success  of  this  device  is  of  primary  importance  to  second-generation  chip 
viability.  Partially  functional  chips  have  been  diced  and  packaged  ready 
for  delivery. 

As  a  major  task  in  the  contract,  a  detailed  redesign  of  the  cell 
architecture  and  array  peripherals  was  performed.  The  basis  for  this 
redesign  was  the  determined  need  to  provide  a  minimum  of  six  storage  sites 
in  the  cell  interior  as  well  as  to  provide  a  coarse  scaling  function  (xl/5) 
in  addition  to  the  FGA.  Fully  parallel  I/O  was  found  to  be  a  primary 
requirement  for  even  small  array  numbers;  further,  a  fully  programmable 
random  access  readout  was  highly  significant  for  overall  system  throughput. 

A  major  breakthrough  was  achieved  by  simply  restructuring  the  storage 
locations  into  a  "racetrack"  or  ring  three-phase  CCD.  This  configuration 
provides  six  storage  locations  while  using  the  inherently  high  packing 
density  of  a  conventional  CCD  shift  register  and  requiring  only  five 
distinct  clock  phases  in  place  of  12  with  the  old  architecture.  With  this 
dramatic  packing  density  improvement  came  the  necessary  real  estate  to 
provide  not  only  the  desired  charge  splitter  (coarse-scaler )  but  also  a 
random  access  analog  output  FET  readable  from  the  periphery  of  the  array. 

All  the  earlier  limitations  on  throughput  have  now  been  overcome,  and  an 
elegantly  simple  solution  to  the  fully  parallel  sensor-to-host  CPU  has  been 
attained  in  concept.  It  is  even  feasible  to  read  data  into  the  array  from  a 
line-parallel  serial  input  concurrently  with  fully  parallel  pixel  data 
directly  into  the  cell.  This  opens  the  way  to  real-time  offset  compensation 
without  the  need  to  store  coefficients  in  the  cell;  an  auxiliary  storage 
device  (either  analog  MNOS  or  digital  ROM)  can-  be  read  once  every  frame. 
Depending  on  the  severity  of  offset  it  may  be  sufficient  to  use  a 
Shutterless  Compensation  Algorithm  (Reference  1) .  No  detailed  studies  were 
attempted  to  determine  feasibility  of  this  approach  because  of  limitations 
on  the  scope  of  the  contract. 


Software  complexity  was  dealt  with  in  two  ways.  First,  a  custom-designed 
8086  microprocessor-based  Exerciser  unit,  benefiting  from  coordinated 
development  efforts  on  other  associated  programs,  was  built  for  delivery  to 
RADC.  This  unit  is  in  practice  a  complete  subsystem  which  not  only  stores 
the  Controller  clock  waveform  primitives  in  RAM — 2K  x  48  in  six  chips  loaded 
from  EPROM — but  can  also  serve  as  the  host  to  interface  to  any  other  system 
through  the  RS232C  port.  A  built-in  expansion  provides  for  MATROX  RGB-256 
video  frame  memory,  allowing  16  grey  levels  NTSC  (or  PAL)  drive  to  a 
high-resolution  (256  x  256  x  4)  monitor  driven  from  a  MATROX  frame 
digitizer.  The  complete  unit  therefore  allows  algorithm  development  in  real 
time  using  any  standard  video  source  (camera,  VTR,  etc.  on  RS170  format)  and 
also  predetermined  test  patterns  from  EPROM/RAM.  To  modify  clock  waveforms 
requires  hook-up  to  a  Microprocessor  Development  system  (e.g.,  iMDS  235), 
but  to  use  the  Honeywell-supplied  primitives  a  simple  keyboard  input  is  all 
that  is  required.  A  Texas  Instruments  (TI)  "Silent  700"  is  supplied  with 
the  unit. 

The  second  aspect  of  software  management  concerns  the  CCD  architecture.  A 
"standard"  was  adopted  for  the  Chip  types  #2214  and  #3022  (firBt  generation) 
requiring  simple  concatenation  of  commands,  each  of  which  has  a  Default  or 
Initial  condition  as  the  first  and  last  clock  states  of  the  sequence.  The 
second-generation  matrix  processor  will  have  a  substantially  greater  degree 
of  concurrency  as  a  designed  feature,  and  so  a  new  "standard"  is  needed. 

The  Exerciser  unit  will  satisfy  both  generations,  since  concurrency  is 
invisible  to  the  user  insofar  as  a  new  set  of  primitives  is  all  that  is 
required  to  be  defined.  The  hierarchy  for  upper  levels  remains  unrestricted. 

Follow-on  development  is  discussed  in  Section  HI  of  this  report. 


SYNOPSIS  OF  THE  REPORT 

The  report  is  organized  in  order  of  the  seven  tasks  which  were  executed.  In 
many  instances,  these  tasks  were  executed  simultaneously.  The  data 
presented  in  each  of  the  task  details  is  not  chronological;  rather,  an 
attempt  is  made  to  order  the  information  to  make  logical  reading  possible. 

Much  has,  by  necessity,  been  omitted  from  the  implications  of  such  a 
powerful  device.  It  is  hoped,  however,  that  the  reader  will  be  left  with 
the  impression  of  a  highly  multi-disciplinary  development,  not  merely  a 
single  chip  design. 


Section  III  gives  some  background  of  Honeywell's  demonstrated  capabilities 
in  indium  bump  direct  sensor  interface  technology.  An  indication  is  given 
of  the  growth  potential  for  a  family  of  lesser  or  greater  capabilities  for 
smart  sensor/seeker  applications  in  the  near  term,  up  to  five  years  from  now. 

As  an  aid  to  brevity,  the  mnemonic  "PIP"  will  be  used  in  place  of  “CCD 
matrix  processor";  however,  this  should  not  be  construed  as  an  implied 
limitation  of  the  device  to  parallel  image  processing. 

LIST  OF  ITEMS  DELIVERED  TO  RADC 


5  quarterly  R6D  status  reports 


Hardware 


15  packaged  CCD  chips;  3022-3-8  *'s  27,  28, 

33,  36,  44,  45,  46,  50,  52, 
55,  58,  70,  72,  82,  91 


,  M,tc  ,,talti~  *ni 

Intel  SBC  86/12  operating  software  in  EPROM 

1  Texas  instruments  Silent  700  Terminal 


Documentation 


!  complete  ..t  ot  .Wicri*  ■»«*“«  “1M 

user  manual  for  Exerciser  unit 

MW.  blt-pettern  -pectnctl-  o<  cco  prUiti...  coeerlc  be.ic 
operation  of  the  matrix  processor  and  at  least  one  edge  ex  r 
(Roberts  Cross)  algorithm 

computer  prcrw  limbi*.  **  *“  “»"s£“  f'mCtl0n  ■=0»'"‘t,tl0,> 
with  constant  capacitance  load 

Final  report 
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SECTION  II 


DETAIL  OF  TASKS  PERFORMED 


TASK  I— CHIP  EXERCISER  DESIGN  AND  FABRICATION 

In  thi*  Mction  the  overall  description  of  the  custom  designed  circuits  for 
operating  tbe  #3022  PIP  is  presented. 

The  Exerciser  unit  comprises  one  custom- fabricated  metal  cabinet  mounted  on 
top  of  a  standard  Intel  660  chassis.  The  custom-fabricated  CCD  clock  driver 
and  signal  conditioning  circuits  are  mounted  on  a  separate  PC  board  mounted 
loosely  on  top  of  the  entire  unit.  This  configuration  was  chosen  as  the 
most  useful  for  access  to  all  the  individual  clock  voltage  control 
potent ioswters  and  was  also  the  least  costly  assembly.  Wire-wrap 
connections  are  used  extensively  throughout  the  unit  to  allow  repairs  and 
any  future  alterations  to  be  easily  implemented. 


A  photograph  of  the  completed  unit  together  with  associated  cameras,  VTR, 
and  monitor  is  shown  in  Figure  1  with  tbe  Intel  MDS  235  in  the  background. 

Circuit  diagrams  showing  parts  and  interconnections  are  detailed  in  the 
engineering  drawings  delivered  to  RADC. 

I/O  Board  Addressing 


In  the  exerciser,  all  boards  are  located  on  I/O  ports  of  tbe  SBC  86/12A 
board  and  are  addressed  by  the  eight  lower  bits  of  the  multibus  address 
lines}  the  upper  bits  are  arbitrary.  Excluding  any  I/O  addressed  functions 
on  the  86/12A  board,  there  are  four  boards  on  the  multibus:  the  RGB-256 
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frame  buffer,  the  FG01/8  digitizer,  the  timing  generator  board,  and  the  I/O 
buffer  board.  The  boards  are  addressed  by  the  5-bit  A?_3.  The  lower 
three  bits  f^2-o  are  used  f°r  functions  on  the  board.  Op  to  32  I/O  boards 
can  be  addressed.  The  addresses  of  these  four  boards  are  listed  in  Table  1. 

Timing  Generator  Board 

The  PIP  clock  generator  is  built  around  an  AMI2910  microprogrammable 
sequencer  capable  of  10  MHz  operation  and  having  a  versatile  programming  mix 
of  branching,  looping,  and  subroutine  instructions.  This  design  will  handle 
almost  all  future  timing  needs. 

A  block  diagram  of  the  timing  generator  board  is  shown  in  Figure  2.  The  key 
elements  are  the  2K  x  48,  90  nsec  random  access  memory  and  a  2910 
sequencer.  The  timing  generator  is  configured  as  a  classical 
double-pipelined  microprogrammable  sequencer.  The  memory  holds 
microinstructions  of  two  types:  clock  generation  and  sequencer  branch 
control.  Under  the  clock  generation  mode  of  operation,  the  memory  can 
simultaneously  generate  states  for  up  to  32  PIP  clocks.  Under  the  sequencer 
branch  control  mode  of  operation,  the  memory  clocks  are  used  for  looping, 
branching,  and  subroutine  calls. 

Timing  Generator  Multibus  Commands 

Following  is  a  brief  summary  of  the  multibus  commands  that  can  be  issued  to 
the  timing  generator  board. 


The  format  is: 


Figure  2.  Block  Diagram  of  Timing  Generator  Board 


The  command  below  loads  the  least  significant  byte  of  the  memory  address  via 
the  multibus  data  bus.  The  format  is: 

BIT 


The  address  is  stored  in  a  register  after  loading,  and  is  not  modified  by  any 
other  command  or  operation. 

The  command  below  loads  the  most  significant  byte  of  the  random  read/write 
memory  address  via  the  multibus: 


The  format  of  the  data  byte  is: 


Multibus 

Bit 

(Bits  8 — 15  Unused) 


The  address  is  stored  in  a  register  that  is  only  modified  by  this  command: 


I/O  Address 

R/W 

Command 

1A 

W  ' 

Write  to  Timing  Generator  Memory  (LSB) 

This  command  writes  the  least  significant  byte  of  data  on  the  multibus  data 
lines  to  the  timing  generator  memory  location  defined  by  18  and  19.  The  data 
byte  format  is: 


6  5  4  3  2  1  0 

7  Most  Significant  Bits 
of  the  11-Bit  Word  Address 


7  6  5  4  3  2  1  0 


Multibus 

Bit 

(Bits  8 — 15  Unused) 


Byte  of  Data  Written  to 
Timing  Generator  Memory 


The  memory  byte  address  is  always  that  specified  by  the  random  address  of  18 
and  19.  The  word  address  of  the  memory  depends  on  the  state  of  DIAG  in  the 
command  register.  If  DIAG  ■  1,  the  random  address  is  used.  If  DIAG  *  Or  the 
sequencer  address  is  used. 

The  command  below  reads  a  data  byte  from  the  timing  generator  memory  into  the 
least  significant  byte  of  the  multibus  data  bus  (the  most  significant  byte 
will  have  arbitrary  contents) : 


I/O  Address 

R/W 

Command 

1A 

R 

Read  Timing  Generator  Memory/Data  Byte 

The  data  byte  format  is: 


Multibus 

Bit 

(Bits  8 — 15  Unused) 


7  6  5  4  3  2  1  0 


Byte  of  Data  from  Timing 
Generator  Memory 


The  memory  address  is  the  same  as  for  the  write  instruction  of  this  command, 
except  when  DIAG  *  0,  the  last  instruction  executed  by  the  sequencer  can  be 
read  from  memory. 

The  command  below  writes  the  least  significant  byte  of  the  multibus  data  bus 
into  the  command  register  on  the  timing  generator  board: 
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I/O  Address 
IB 


R/W 

W 


_ Command _ 

Write  Command  Register  Byte 


The  format 

where 

RESET 

RUN 

SS 

ENDVR 

OIAG 

CLKRT 


J 


of  the  command  byte  is: 


7  6  5 

4 

3 

2 

1 

0 

Multibus 

Clock 

D 

E 

S 

R 

R 

Bit 

Speed 

I 

N 

S 

U 

E 

(Bits  8 — 15  Unused) 

2-0 

A 

D 

N 

S 

G 

V 

E 

R 

T 

*  0  resets  all  board  functions  (set  to  0  on  power-on) 

*  1  disables  reset  function 

*  1  sets  the  sequencer  to  run  mode 

0  halts  the  sequencer  (set  to  0  on  power-on) 

■  1  single-steps  the  sequencer  one  instruction  (RUN  *  0) 

-  0  no  single  step  (set  to  0  on  power-on) 

SS  must  be  toggled  to  step  multiple  instructions 

■  (not  used  at  present) 

*  1  selects  the  random  address  as  the  timing  generator  memory 

address  (RUN  must  equal  zero) 

*  0  selects  address  from  sequencer  as  timing  generator  address 

(set  to  0  on  power  on) 

-  000  clock  sequencer  runs  at  8  MHz 

-  001  4.16  MHZ 

-  010  1.78  MHz 

-  011  890  kHz 

-  100  445  kHz 
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101  223  kHz 


*  110  111  kHz 

*  111  55  kHz 

The  reset  true  condition  will  cause  a  jump-to-zero  instruction  to  be  put  on 
the  2910  instruction  lines. 

The  command  below  loads  the  least  significant  byte  of  the  multibus  data  bus 
into  the  least  significant  byte  of  the  11-bit  macro  address  register: 


I/O  Address 

R/W 

Command 

1C 

W 

Load  MACRO  Address 
LSB  (Byte) 

The  MACRO  address  is  a  word-address  only  and  has  no  byte  address.  The  format 
of  the  data  byte  is: 


7  6  5  4  3  2  1  0 


Multibus 

Bit 

Bits  8-15 
Unused) 


Least  Significant  Byte  of  the  11-bit 
MACRO  Address 


Bit 


This  address  can  be  continuously  written  to  the  board.  Wait  states  (until  the 
sequencer  uses  the  address)  are  automatically  generated  and  are  transparent  to 
the  user.  The  macro  address  LSB  can  only  be  changed  by  this  command.  If  the 
timing  generator  is  in  a  single-step  mode,  outputting  this  command  will 
generate  a  single-step  cycle. 


The  command  below  reads  onto  the  multibus  data  bus  the  least  significant  word 
current  address  on  the  timing  generator  for  memory: 
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I/O  Address 

R/W 

Command 

1C 

R 

Read  Least  Significant  Byte  of 
the  Timing  Generator 
Memory  Address 

The  source  of  the  address  is  determined  by  the  state  of  01 AG  in  the  command 
register.  If  DIAG  *  1,  then  the  address  is  the  word-address  portion  of  the 
random  address,  if  DIAG  -  0,  then  the  address  is  from  the  sequencer  (not  to 
be  confused  with  the  macro  address) .  The  format  of  the  data  to  the  multibus 
is  given  below: 

7  6  5  4  3  2  1  0  Bit 

Multibus  Least  Significant  Byte  of 

Bit  Current  Timing  Generation 

(Bits  8-15  Memory  Address 

Arbitrary  - - - — - — 


7  6  5  4  3  2  1  0 


Least  Significant  Byte  of 
Current  Timing  Generation 
Memory  Address 


The  command  below  loads  the  least  significant  byte  of  the  multibus  data  bus 
into  the  most  significant  bits  of  the  11-bit  macro  address: 


I/O  Address 

R/W 

Command 

ID 

W 

Load  MACRO  Address 

MSB 

The  format  of  the  data  byte  is: 


7  6  5  4  3  2  1  0  Bit 


Multibus 

UNUSED 

0 

3  Most  Significant 

Bit 

MACRO  Address 

(Bits  8-15 

.  1 

Bits 

Unused) 


Reserved 

for 

Expansion 
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This  address  can  be  continuously  written  to  the  board.  Wait  states  are 
automatically  generated.  The  macro  address  MSB  can  only  be  changed  by  this 
command. 

The  command  below  reads  onto  the  multibus  data  bus  the  most  significant  bits 
of  the  current  timing  generator  memory  word  address: 


I/O  Address 

R/W 

Command 

ID 

R 

Read  the  Most  Significant  Byte 
of  the  Timing  Generator 
Memory  Address 

The  source  of  the  address  is  determined  by  DIAG,  as  explained  in  the  LSB 
address  read  command  (100).  The  data  format  is: 

7  6  5  4  3  2  1  0  Bit 


Multibus 

UNUSED 

0 

3  Most 

Data 

(Zero) 

Significant 

(Bits  8-15 

Address  Bits 

Arbitrary 

Reserved 


Timing  Generator  2K  x  48  Bit  Memory 

There  are  two  48-bit  formats  for  the  microinstruction,  as  shown  in  Figures  3 
and  4.  Figure  3  is  the  clock  generation  format  and  Figure  4  is  the  branch 
address  format,  in  clock  generation,  32  bits  of  the  instruction  indicate 
clock  states  (C1  to  C32> .  The  specific  PIP  clock  assignments  in  this 
instruction  are  given  in  Table  2.  The  2910  instruction  (I3  Q)  must  not  be  a 
specified  address  branch  type  (but  it  can  be  a  stack  pull  or  subroutine 
return).  The  instruction  is  usually  a  continue.  The  condition  code  (CC^  Q) 
is  unused.  The  enable  bit  (EM),  when  1,  loads  the  clock  states  into  the  clock 
register  (to  go  to  external  devices)  on  the  next  sequencer  instruction.  If 
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16  1? 


=  SINGLE  PIPELINE  1  «  DOUBLE  PIPELINE 


TABLE  2.  CLOCK  ASSIGNMENTS 


3022 

PIP 

Clock 

PIP 

Mnemonic 

(for  programming) 

Default 

Binary 

State 

Timing  Generator 
Microinstruction 
Bit  position 

°1 

PHI 

H (high) 

16 

°2 

PH2 

L (low) 

17 

°3 

PH3 

L 

18 

°4 

PH4 

L 

19 

0 

A 

PHA 

L 

20 

°B 

PHB 

L 

21 

°C 

PHC 

L 

22 

0 

D 

PHD 

L 

23 

0 

X 

PHX 

H 

24 

°xs 

PHXS 

H 

25 

STl 

STl 

H 

26 

ST2 

ST2 

H 

27 

ST3 

ST3 

H 

28 

vs 

VS 

H 

29 

°RS 

PHRS 

L 

30 

FGR 

FGR 

H 

31 

RSNI 

RSNI 

L 

32 

RSI 

RSI 

L 

33 

°T 

PHT 

L 

34 

°XM 

PHXM 

L 

35 

md 

MD 

H 

36 

EN  is  0,  the  clock  register  will  not  be  loaded  (EN  is  almost  always  1  for 
this  instruction  format).  The  pipeline  mode  bit  (PLM) ,  when  0,  runs  the 
2910  as  a  single-pipeline  machine.  If  PLM  is  1,  the  2910  runs  as  a 
double-pipeline  machine.  The  only  constraint  to  switching  modes  is  that, 
when  going  from  single-  to  double-pipeline,  the  first  double-pipeline 
instruction  executes  twice. 

The  instruction  format  in  Figure  4  is  used  for  all  branching  or  looping 
instructions.  A  12-bit  address  (BAn_0)  is  available,  but  only  11  bits 
are  used  at  present  (i.e.,  BAj^  ■  0) .  The  instruction  can  be  any 

of  the  16  used  by  the  2910.  The  only  instruction  that  is  not  implemented  i 
a  standard  way  is  the  conditional  jump  vector  (CJV) ,  which  will  be  used  for 
macro  address  branching.  The  condition  code  on  a  CJV  should  always  be 
zero.  Also,  a  jump-to-zero  instruction  should  always  follow  the  CJV.  The 
CJV  should  always  be  used  at  address  zero  in  the  timing  generator  in  the 
single-pipeline  mode. 

The  branch  instruction  has  a  PIP  multiplier  field  (PM^g)  in  bits  32  to 
39.  This  serves  as  a  parameter  when  a  CJSB  2910  instruction  is  used.  The 
target  subroutine  can  output  the  field  to  the  PIP  multiplier  port  using  one 
of  the  special  clocks. 

In  both  microinstruction  formats,  eight  special  clocks  (SPCK^_8>  are  used 
for  various  operations,  such  as  one  I/O  buffer  sequential  read/write 
operation  for  PIP  I/O. 


Double-Pipeline  Branching 

Double-pipeline  mode  is  normally  only  used  when  instructions  are  to  be 
executed  consecutively.  Branching  is  usually  done  in  single-pipeline  mode, 
but  can  be  done  in  double-pipeline  mode  also.  If  a  branch  to  address  X  is 
executed  at  location  N,  the  instruction  in  location  N  +  1  is  also  executed 
due  to  the  double  pipeline.  The  instruction  at  N  +  1  must  not  be  a  branch. 


Timing  Generator  Memory  Map 

The  organizing  of  the  timing  generator  memory  ia  to  have  groups  of  clock 
sequences  broken  into  functions.  The  functions  will  be  linked  via 
subroutine  calls  in  macro  sequences.  The  macro  sequences  will  be  linked  by 
address  in  the  macro  address  register  put  out  at  a  high  speed  by  the  SBC 
86/12A  (see  Figures  2  and  3).  A  CJV  instruction  in  location  0  (branch  to 
zero  executed  on  reset  from  the  command  register)  polls  the  macro  register 
until  it  is  loaded  by  the  86/12A.  When  it  is  loaded,  a  branch  to  that 
address  occurs.  A  sequence  of  subroutine  calls  (JSB)  are  executed  to  form  a 
macro  sequence.  The  macro  sequence  terminates  with  a  jump-to-zero  command. 
The  subroutines  are  clock-state  sequences,  which  generate  the  clock  waveform 
patterns  for  a  specific  function.  These  sequences  end  with  a  return  from 
subroutine.  The  clocks  are  output  at  high  speed  in  the  double-pipeline 
mode,  while  all  branching  is  done  in  the  single-pipeline  mode. 

Condition  Codes  (CC  ^  Q) 

There  are  eight  possible  condition  codes,  but  only  two  are  currently 
implemented.  They  are  described  in  Table  3. 

I/O  Buffer  Board 

The  I/O  buffer  is  shown  in  Figure  5.  The  board  stores  input/output  data 
to/from  the  PIP  and  does  digital/analog  and  analog/digital  conversions. 
Buffer  1  is  loaded  via  the  multibus  with  data  going  to  the  pIP  (up  to  4K 
words  of  8  bits  each)  .  Buffer  2  is  loaded  with  data  from  the  PIP  and 
interrogated  by  the  multibus.  Both  buffer  memories  can  be  read  by/written 
to  the  multibus  for  diagnostic  purposes. 
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TABLE  3.  CONDITION  CODES 


ccl-0 

Description 

00 

True: 

the  macro  address  register  has  been  loaded  by 
the  multibus,  but  its  contents  have  never  been 
branched  to  (i.e.,  has  not  been  "unloaded") 

False: 

macro  address  register  is  either  "empty"  or 
has  been  branched  to  ("unloaded") 

01 

True: 

always  true  (unconditional  branch) 

False: 

never 

10 

(Not  assigned) 

11 

(Not  assigned) 

As  Car  as  Che  multibus  is  concerned,  the  two  buffer  memories  appear  as  a 
single  address  8K  x  8  memory,  with  the  lower  41  x  8  as  buffer  1  and  the 
upper  4K  x  8  as  buffer  2. 

I/O  Buffer  Board  Multibus  Commands 

The  following  is  a  brief  summary  of  the  multibus  commands  that  can  be  issued 
to  the  I/O  buffer  board: 


I/O  Address 

R/W 

Command 

20  (HEX) 

W 

Load  Memory  Address  (LSB) 

This  command  loads  the  least  significant  byte  of  the  13-bit  memory  address 
via  the  multibus  data  bus. 


CLOCKS  FROM  INTEL  MULTIBUS 

FOCAL  PLANE 


The  format  of  the  data  byte  is: 


7  6  5  4  3  2  1  0 

Multibus 
Bit 

(Bit*  8-15  unused) 


8  LSB  of  Memory  Address 


This  address  is  loaded  in  a  register  that  can  be  modified  if  external 
control  is  enabled  in  the  command  register.  Otherwise  it  does  not  change, 
except  by  this  command. 

The  command  below  loads  the  five  most  significant  bits  of  the  13-bit  memory 
address  from  the  multibus  data  buss 


I/O  Address 

R/W 

Command 

21(HEX) 

w  • 

Load  Memory  Address  MSB 

The  format  of  the  data  byte  is: 


Multibus 

Bit 

Bits  8-15 
Unused) 


7  6  5  4  3  2  1  0 


Not 

5  MSB  of  13-bit 

Used 

Memory  Address 

where  BIT  4*0  selects  buffer  memory  1  and  BIT  4  *  1  selects  buffer  memory  2 


The  command  below  writes  the  least  significant  byte  of  the  multibus  data  bus 
to  the  i/o  buffer  memory  at  the  previously  specified  address: 


I/O  Address 

R/W 

Command 

22  (HEX) 

W 

Write  Memory  Data  Byte 

The  data  byte  format  is: 

7  6  5  4  3  2  1  0  Bit 


Multibus 

Bit 

(Bits  8-15 
Unused) 


Data  Byte  Written 
To  I/O  Buffer  Memory 


The  command  below  reads  onto  the  least  significant  byte  of  the  multibus  data 
byte  the  contents  of  the  I/O  buffer  memory  at  the  previously  specified  address 


I/O  Address 

R/W 

Command 

22  (HEX) 

R 

Read  Memory  Data  Byte 

The  resulting  data  byte  format  is: 


7  6  5  4  3  2  1  0  Bit 


Multibus 

Bit 

Byte 


Data  Byte  Read  from  I/O 
Buffer  Memory 


The  command  below  writes  the  least  significant  byte  of  the  multibus  data  bus 
into  the  I/O  buffer  command  register: 
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R/W 


Command 


I/O  Address 


23  (HEX) 


W 


Write  Command  Register 


The  data  format  is: 

7  6  5  4  3  2  1  0  Bit 


Multibus 

Not 

Not  Specified 

E 

Bit 

Used 

N 

(Bits  8-15 

E 

Unused) 

X 

c 

wnere 

ENEXC  ■  0  external  control  disabled  (All  control  signals  to  the  I/O 

buffer  board  are  disabled  except  those  from  the  multibus) 


■  1  enable  external  control 


[External  control  signals  to  transfer  PIP 
data  (from  timing  generator)  are  enabled. 

Only  the  multibus  command  register  command 

should  be  used.  Other  multibus  commands  can  provide 

erroneous  results.] 

The  external  controls  are  generated  as  clocks  on  the  timing  generator  board 
and  by  the  PIP  chip. 
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TASK  2 — INITIAL  CCD  CHARACTERIZATION 

The  initial  device  characterization  work  was  performed  using  Honeywell's  pip 
chip  #2214,  which  had  been  fabricated  earlier  in  1980  before  the  start  of  the 
contract.  Only  one  run  of  wafer  processing  was  made  at  that  time,  resulting 
in  a  very  limited  sample  of  devices.  However,  a  thorough  enough  analysis  of 
the  2214  chip  was  made  that  detailed  changes  were  seen  to  be  necessary. 

The  following  subsection  describes  the  array  and  cell  design  of  the  2214 
chip.  Measurements  from  these  chips  are  then  described  to  indicate  the  need 
for  the  changes  that  were  made.  Finally,  the  changes  to  the  array  and  cell 
are  described. 

Original  PIP  Array  and  Cell  Design 

A  prototype  cell  was  used  to  generate  the  whole  array.  Access  to  the  array  is 
through  a  f ill-and-spill  serial  input  at  the  top  left,  a  four-phase  serial  I/O 
buried-channel  shift  register  (with  90°  corner  gate)  along  the  left  and 
bottom  edges  of  the  array,  and,  finally,  a  floating  diffusion  sense  amplifier 
with  reset  at  bottom  right. 

A  polysilicon  clock-bus  set  of  lines  runs  down  the  right  side  of  the  array  to 
connect  the  repeated  metal  clock  lines  which  run  left-right  through  the 
array.  These  are  visible  in  the  photograph  of  the  chip  shown  in  Figure  6. 

A  copy  of  the  Calma  plot  of  the  cell  is  shown  in  Figure  7,  together  with  a 
photograph  of  the  cell  on  the  chip.  The  cell  contains  three  separately 
controlled  storage  areas  labeled  STl,  ST2 ,  and  ST3  in  Figure  7.  ST2  and  ST 3 
have  one  common  additional  gate,  XS,  which  is  used  to  control  access  to  these 
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Figure  6.  Photograph  of  Original  #2214  CCD  chip 


stores.  Also  contained  in  the  cell  is  a  surface-channel  f ill-and-spill  input 
section  labeled  Regenerator.  The  inputs  to  the  regenerator  are  obtained  from 
two  floating-gate  detectors  labeled  Y  and  Z.  Each  floating  gate  has  a  reset 
transistor  source  control--RSI  and  RSNI--but  the  gates  of  the  transistors  are 
tied  together  and  controlled  by  FGR.  The  Z  detector  is  loaded  by  the 
multiplier  capacitor  whose  control  gate  is  labeled  MG.  The  multiplier  in  this 
design  is  equivalent  to  a  surface-channel  gated-diode  structure  in  which  the 
diode  is  connected  to  the  floating-gate  detector.  By  biasing  the  MG  gate 
towards  surface  inversion,  the  total  capacitive  load  on  the  detector  is 
increased.  The  sensitivity  of  the  detector  is  correspondingly  decreased. 
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Figure  7.  Calm*  Plot  of  the  2214  Cell,  Together  with  a  Photograph 
of  the  Actual  Cell  (Poly  I  is  shaded) 
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Finally,  the  cell  contains  a  sink  diode,  ARV,  for  erasing  data  from  the 
array  in  parallel. 

The  periphery  of  each  cell  is  occupied  by  two  three-phase  buried-channel 
shift  registers  with  four  distinct  phases  in  groups  (A,B,C)  and  (B,C,D) . 

The  top  right  corner  is  occupied  by  the  matrix  corner-turning  phase,  labeled 
X. 


The  remaining  area  of  the  cell  is  used  to  obtain  clock  phase  interconnection 

with  two  polysilicon  levels  and  one  metal  level.  Approximately  30%  of  the 

cell  area  excluding  shift  registers  is  used  for  this  purpose  in  the  present 

2 

conservative  design.  The  2214  cell  dimensions  are  265  x  265 ym  (10.43 

2  2 
x  10.43  mil  )  and  the  chip  overall  dimensions  are  230  x  200  mil  . 

Operation  of  the  regenerator  is  shown  schematically  in  Figures  8 (a)  and 
8(b).  The  f ill-and-spill  diode  is  connected  to  the  C  clock,  which  is  the 
phase  into  which  the  generated  charge  flows  directly.  This  scheme  has  been 
demonstrated  at  RSRE,  Malvern,  England  (Reference  2) ,  to  give  a  linearity 
equally  as  good  as  the  more  conventional  method  while  eliminating  the 
separate  dio  3  clock  pulse.  The  SCCD  to  BCCD  transition  occurs  under  the 
extension  of  the  C-phase  gate.  This  provides  a  convenient  potential  barrier 
to  prevent  charge  sloshing  back  into  the  regenerator  when  the  shift  register 
is  used  independently. 

It  should  be  noted  that  the  edges  of  the  array  are  terminated  only  on  the 
left  and  bottom  by  shift  registers.  At  the  top  and  right  there  is  nowhere 
for  charge  to  be  shifted  upwards  or  to  the  right.  The  significance  of  this 
was  not  apparent  until  the  devices  were  operated.  This  is  described  in  the 
following  subsection,  on  measurements. 
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Figure  8.  Regenerator  Operation 


Tbe  modifications  to  the  cell,  described  later,  were  aimed  at  eliminating 
problems  with  performance  degradation  rather  than  changing  tbe  functionality 
of  the  cell  or  the  array. 

Initial  Measurements 


The  clock  timing  for  operating  PIP  is  evidently  quite  intricate.  The  cell 
shift  registers  are  used  bidirectionally  to  maneuver  charge  packets  into  and 
out  of  the  ports  of  the  cell.  An  example  of  a  maneuver  to  cause 
regeneration  of  the  contents  of  ST1  is  shown  in  Figure  9.  Since  tbe  number 
of  clocks  for  control  is  so  large,  it  is  necessary  to  use  a  simple  model  of 
the  cell  to  observe  in  real  time  the  action  of  a  given  logical  state. 

The  example  primitive  shown  in  Figure  9  illustrates  the  machine  code 
constraints  on  hierarchical  software  structure.  It  is  necessary  for  the 
final  state  to  be  identical  to  the  initial  state  so  that  primitives  can  be 
concatenated.  This  is  not  necessarily  the  most  efficient  way  to  use  the 
array  since  a  "next  destination"  may  already  be  available  in  the  front  end 
of  a  subsequent  primitive.  These  questions  are  addressed  in  more  detail  in 
Task  5  on  Software. 

The  basic  operation  of  the  array  can  be  observed  with  simple  input 
conditions.  For  example,  if  the  serial  input  is  off  the  only  charge  that 
can  be  read  out  of  the  array  is  regenerator  offset  and  leakage.  Figure  10 
shows  the  output  from  such  an  operating  condition  with  the  regenerator 
controls,  RSI  and  R'SI,  adjusted  to  give  no  offset,  with  a  frame  period  of 
23  msec.  Notice  that  the  peaks  are  at  tbe  leading  end  of  each  line,  and 
there  are  16  peaks.  This  charge  should  ideally  have  been  shifted  off  tbe 
end  of  the  array  during  the  I/O  sequence,  but  instead  has  summed  at  the 
rightband  cell  in  any  convenient  "on"  gate.  Tu«  charge  merely  spills  to  tbe 
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Figure  9.  Clock  Waveform  Timing  Diagram  for  One  Primitive  (Phases 
not  shown  are  at  rest  in  initial  state.) 


Figure  10.  Leakage  Residual  Buildup  at  Frame  Period  of 
23  msec,  25°C  (2214) 

lowest  potential  if  there  is  no  gate  to  enter  normally.  This  defect  can  be 
eliminated  by  providing  extra  cell  gates  at  the  top  and-  righthand  edge  of 
the  array  plus  sink  diodes  to  remove  the  charge. 

The  linearity  and  sensitivity  of  the  serial  input  I/O  four-phase  register 
and  sense  amplifier  are  shown  in  Figure  11  with  a  triangular  input. 

These  basic  functions  do  not  involve  arithmetic  operations  of  the  cell.  A 
rather  complex  operation  using  all  the  features  of  the  array  is  shown  in 
Figures  12  and  13.  The  array  is  loaded  with  a  6  x  6  uniform  pattern  and 
then  the  following  sequence  is  performed: 

1.  ERASE  STl 

2.  ERASE  ST2 

3.  ERASE  ST 3 

4.  LOAD  STl 


} 


(input  data  retained  in  x) 


Figure  IX.  Input  Fill  and  Spill  to  Sense  7-amp  performance  Using 
Serial  I/O  Registers  Only 


Figure  12.  XY  Display  Using  the  Routine: 

El  E2  E3  Ll  Rl  L2  R1  L2  U2 
SU  SU  SU  SL  SL  SL  L2  M2 
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Figure  13.  video  Waveform  at  Top  Line  of  6  x  6  Test  pattern 
(Zero  reference  at  center  of  screen.) 

5.  regenerate  STl  (copy  of  ST1  retained  in  x) 

6.  LOAD  ST2 

7.  REGENERATE  STl  (copy  of  STl  retained  in  X  again) 

8.  LOAD  ST2 

At  this  point  register  addition  bas  occurred. 

9 .  UNLOAD  ST2 

10.  SHIFT  UP  3  TIMES 

11.  SHIFT  LEFT  3  TIMES 

At  this  point  the  double  copy  of  the  original  pattern  is  now  shifted  to  the 
upper  left  corner  of  its  original  position.  The  original  data  remains  in  STl 
and  in  the  original  position. 

12.  LOAD  ST 2 

13.  SUBTRACT  ST2  FROM  STl 


This  last  primitive  causes  the  difference  between  the  charge  in  ST2  and  STl  to 
be  generated  and  placed  in  X.  The  result  is  then  observed  by  executing  an  I/O 
parallel-to-serial  readout  sequence. 

The  result  of  this  sequence  of  primitives  should  be  a  6  x  6  square  with  one 
3x3  corner  reduced  by  subtraction.  All  other  output  locations  should  be 
zero.  Notice,  however,  that  the  shift  operations  in  both  axes  have  caused  the 
leakage  charge  to  appear  both  at  the  bottom  of  the  display  and  the  right 
side.  By  cooling  the  array  it  is  possible  to  eliminate  these  spurious  outputs 
almost  entirely.  The  major  importance  of  this  result  is  that  the  array  is 
able  to  perform  the  primary  functions  in  a  manner  very  close  to  that  desired. 

Figure  14  shows  a  typical  output  video  from  the  original  tester  buffer 
circuits  directly  after  the  CCD.  The  circuits  are  simply  impedance  and  offset 
buffers  to  provide  drive  for  the  S/H  amplifier  driving  the  data  logging 
system.  No  filtering  or  integration  has  been  applied.  The  low  gain  of  the 
2214's  regenerator  is  evident  when  compared  to  an  unregenerated  output  signal 
of  over  IV.  Figure  14(a)  shows  an  entire  frame  (256  pixels)  and  Figure  14(b) 
shows  several  lines  (16  pixels  each)  in  which  data  is  nonzero.  Nonuniformity 
is  clearly  visible. 

Leakage  current  and  transfer  inefficiency  data  have  been  analyzed,  together 
with  various  uniformity  performance  data.  Typical  data  are  shown  in  Figures 
15  to  31,  which  have  been  logged,  plotted,  and  analyzed  using  an  on-line  HP-IB 
desktop  computer  system.  Since  the  quantity  of  the  data  is  so  large  for 
two-dimensional  arrays  (even  a  small  16  x  16) ,  only  the  major  points  will  be 
covered. 
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a.  Whole  frame  video  (256  pixels)  of  SUBTRACTOR 


b.  Details  of  three  lines  (16  pixels  each) 


Figure  14.  Typical  Video  output  After  Buffer 


Figure  15  shows  the  result  o£  shifting  zero  data  to  the  right  using  the 
instruction  set  ,  E2,  E3»  I/O  {load)  I/O  (unload).  We  expect  to  see 
an  accumulation  of  leakage  charge  in  the  16th  column  as  a  result  of  the  lack 
of  terminations  to  the  array.  However,  although  there  is  evidently  about  60 
mV  output  at  the  appropriate  positions  for  the  16th  column  (pixel  #1  is  bottom 
right,  pixel  #16  is  bottom  left,  etc.,  pixel  #256  is  top  left),  the  output  is 
drastically  smaller  than  that  shown  in  Pigure  16.  For  Figure  16  the 
instruction  set  was  changed  to  E^,  E2»  E^  SU(16)  I/O  (unload).  This 
shows  the  expected  800  mV  output  corresponding  to  the  top  row  lack  of 
termination. 

The  anomaly  was  resolved  by  the  presence  of  an  observed  layout  error  in  the 
16th  column  shift  register.  There  is  a  phase  sequence  connection  fault  on  the 
16th  column  giving  an  (X,  X,  B,  C,  B,  D,  ...)  instead  of  the  intended 
(X,  D,  C,  B,  D, ...).  This  simple  isolated  error  was  corrected  on  the  #3022 
chip  mask  set. 

Figures  17  through  27  show  the  leakage  charge  accumulated  in  each  cell  store 
when  the  routines 

E2  E3  U3  I/O  (unload), 

E^  E3  U3  I/O  (unload),  and 

E^  E2  U2  I/O  (unload) 

are  performed  both  at  room  temperature  C  25°C)  and  cooled  f  °C) .  The 
integration  time  in  each  case  is  about  23  msec. 

-6  2 

The  gate  area  of  ST1  ■  6.9  x  10  cm  , 

■■6  2 

of  ST 2  >  7.2  x  10  cm  ,  and 
of  ST3  »  7.74  x  10  ®  cm^. 
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PIXEL  # 

Figure  15.  Zero  Data  Leakage  Measurement :  Shift  Right  (2214) 


Figure  16.  Zero  Data  Leakage  Measurement:  Shift  up  (2214) 


1/23/^81 


Figure  17.  Integration  Mode  Leakage  of  ST1:  Room  Temperature 


TORE  #1  Leakage  Cooled  2/5/81 
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Figure  19.  Integration  Mode  Leakage  of  STl:  Cooled  0°C 


STORE  *2  Leakage  1/28/81 


Figure  21.  integration  Mode  Leakage  of  ST2 :  Room  Temperature 


A0N3D03dJ  3AIlhTGd  lN3DM3d 
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Figure  22.  Histogram  of  Data  in  Figure  21 


TORE  #2  Leakage  Cooled  2/6/81 


Figure  23.  Integration  Mode  Leakage  of  ST2:  Cooled"  0°C 


STORE!  *3  Leakage  2/28/81 


Integration  Mode  Leakage  of  ST3:  Room  Temperature 


STORE  #3  Leakage  1/28/81 


Figure  25.  Histogram  oC  Data  in  Figure  24 


TORE  #3  Leakage  Cooled  2/5/81 


Figure  26.  Integration  Mode  Leakage  of  ST3:  Room  Temperature 


One  volt  output  corresponds  to  about  1  pC  charge  packet 


Evidently,  STl  uniquely  suffers  from  a  nonuniform  leakage  current.  The 

magnitude  of  the  leakage  current  at  25  C  is  about  2.5  y  A  cm“z  and 

-2 

greatly  exceeds  the  expected  value  of  10  to  100  nA  cm  .  ST2  and  ST3  are 

_2 

close  to  the  expected  leakage  at  about  300  nA  cm  and  are  uniformly 
distributed  (gaussian). 

The  only  physical  correlation  with  the  leakage  behavior  of  STl  that  could 
account  for  a  nonuniform  signature  is  a  nonuniform  proximity  of  the  diffusion 
feeding  the  c  clock  to  the  regenerator  f ill-and-spill  diode.  He  have  searched 
for  such  a  variation  without  much  success  using  an  SEM  on  actual  chips.  The 
diode-poly  spacing  was,  however,  designed  to  be  4ym,  which  is  too  close  for 
comfort  with  snort  channel  effects  evident  in  FETs  at  that  gate  length. 

The  3022  chip  has  double  that  spacing,  which  should  eliminate  the  effect 

entirely.  It  should  be  noted  that  in  normal  operation  of  the  array,  we  erase 

all  stores  immediately  after  loading  the  array  using  I/O  (load);  The  leakage 
should  not,  therefore,  be  a  problem  even  at  its  evident  magnitude  at  25°C. 

Histogram  plots  of  each  set  of  measurements  are  shown  in  sequence  with  the 
pixel  distribution. 

Figure  28  shows  a  plot  of  the  2214' s  response  to  simple  loading  and  unloading 
of  the  perimeter  I/O  registers  with  a  6  x  6  block  of  unit  data. 

Transfer  inefficiency  effects  will  be  the  only  contribution  to  distortion 
here.  Evidently  nothing  very  significant  has  occurred  with  the  2214's  (16  x  8 
and  16  x  8  +  4  »  260)  transfers.  Figure  29  is  a  plot  of  the  same  data  as 
Figure  28,  to  show  the  fat  zero  level  of  “  8%  and  the  trailing  inefficiency 

product  of  about  25  mV.  This  corresponds  to  an  inefficiency,  E ,  of  about  8 

x  10-5  tot  ■  0.99992)  at  a  clock  rate  of  about  100  kHz — the  expected 
performance  of  a  BCCD. 
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ransfer  Inef f 1 c 1 ency  2/16/8 


Figure  28.  Transfer  Inefficiency  Effects  in  I/O  Registers 


Figure  30  shows  the  same  test  as  Figure  2 d,  but  this  time  the  array  is  loaded 
in  audition  to  shitting  through  the  I/O  registers.  Now,  we  observe  the 
leakage  accumulation  outputs  (16th  column  missing)  and  again  a  good e  o£ 
about  the  same  value.  The  number  of  transfers  in  this  case  varies  depending 
on  the  location  of  the  data,  but  is  roughly  equivalent  to  the  previous  case. 
Evidently,  corner-turning  maneuvers  are  no  problem  with  the  2214. 

There  are  no  unexplained  phenomena  in  the  testing  of  the  2214  concerning 
leakage,  transfer  inefficiency,  and  uniformity.  These  results  and  the 
analysis  of  the  defective  properties  were  taken  into  account  in  the  redesign 
of  the  cell  and  array  termination,  described  next. 

TASK  3 --DESIGN  ITERATION  AND  FABRICATION  OF  16  x  16  ARRAY 

Design  Modif ications 

Several  major  modifications  were  made  to  the  cell  and  one  was  made  to  the 
array  termination. 


Cell— 

1.  Regenerator  gate  structure  was  changed  to  reduce  loading  effects  on 
floating  gates  z  and  Y. 

2.  Charge-charge  sensitivity  of  the  regenerator  was  intended  to  be 
increased  from  presently  observed”  0.2  to”  1.0. 

3.  Input  gates  to  the  regenerator  were  screened  from  c-clock 
breakthrough. 


4.  Spacing  between  STl  and  c-clock  diode  was  increased  to  reduce 
"leakage"  into  STl. 


Figure  30.  Transfer  Inefficiency  Effects  in  I/O  and  Array  Registers 


Trans-far  Ihpff  I  c  i  ancM  3/16/81 
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Figure  31.  Magnified  Scale  of  Data  in  Figure  30 


5.  The  multiplier  was  changed  to  produce  a  linear  transfer  function  in 
the  regenerator.  The  present  design  gives  "square  law"  distortion 
products.  This  new  Honeywell  design  has  been  described  in  another 
RADC  Contract,  the  CCD  Applications  Study  (Reference  3). 

Dummy  cell  gates  (A,B)  and  new  phase  T  were  added  to  the  right  edge 
termination,  and  a  diode  was  connected  to  ARV.  Similar  dummy  gates 
(B,C) ,  new  phase  T,  and  an  ARV  diode  were  added  to  the  upper  edge 
termination.  This  will  eliminate  the  major  effects  of  thermal 
leakage  charge  buildup  observed  in  measurements. 

The  new  generator  structure  is  shown  in  cross-section  in  Figure  32. 
Introduction  of  screen  gate  connected  to  RSI  gives  a  slightly  different 
arithmetic  operation.  The  output  consists  of  three  terms: 

l-  <W  ■  K1  °Z  -  KA  *  «A 

■  K1«Z  -  K2«y 

2.  and  K2  are  the  floating-gate  charge-to-charge  sensitivities. 

3.  is  an  extra  charge  input  due  to  quasi-DC  screen  gate  RSI. 

Figures  33  to  36  show  the  new  chip  cell  layout  #3022}  and,  as  a 
comparison,  the  old  cell  layout  (#2214)  .  The  shaded  regions  are  Poly  1 
in  both  Figures  33  and  34,  and  the  metal  interconnect  is  shown  separately 
in  Figures  35  and  36.  Figure  36  indicates  the  location  of  the  two  extra 
control  lines  introduced  to  provide  the  improved  linear  multiplier  (MD  + 

MG2)  and  also  the  exchange  of  memory  control  (Xm) .  The  shift  registers 
are  not  altered  from  the  original  design  save  for  the  different  locations 
of  the  contact  windows  and  Poly  1,  Poly  2  interconnect.  The  new 
regenerator  is  seen  to  have  appropriate  screen  gates  as  described  above. 
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Figure  32.  Modified  Regenerator  Gate  Structure  Showing  Reduced  Loading 
Increased  Sensitivity,  and  C-Clock  Screening 


Figure  37  shows  the  modified  (approximately  to  scale)  40-pin  DIP 
assignments.  Also  noticeable  is  a  more  uniform  pad  placement,  which 
allows  rectangular  placement  of  the  die  in  the  DIP  cavity.  Test  FETs  are 
retained  to  provide  VT  uniformity  data  for  comparison  with  the  array 
regenerator  uniformity  measurements,  but  these  are  not  bonded  out  to  the 
package. 

Although  not  shown  in  the  cell  plots,  the  array  is  terminated  by  dummy 
phase  gates  and  a  sink  diode  at  each  column/row  end.  This  ensures  that 
attempts  to  shift  leakage/data  charge  off  the  array  will  result  in  zero 
charge  remaining  in  those  peripheral  rows/columns. 

The  terminations  use  separate  diodes  and  gates  in  the  same  manner  as  the 
SINK  node  in  the  cell.  An  extra  phase,  T,  was  introduced  to  perform  this 

The  sequence  of  clock  states  for  a  Shift  Right  (SR)  is  shown 
diagramatically  in  Figure  38  and  for  shift  up  (SU)  in  Figure  39. 

improved  Functionality 


The  *3022  cell  architecture  is  essentially  the  same  as  the  #2214. 

However,  there  are  two  major  differences  in  functionality.  The  first  is 
the  new  linear  regenerator,  which  was  described  above.  The  second  is  the 
extra  control  clock  XM,  which  connects  ST 2  and  ST3 .  This  allows  two  new 
primitive  command  types:  rotate  the  contents  of  X,  ST2,  and  ST3 
(bidirectional);  swap  the  contents  of  ST2  and  ST3.  "Move"  primitives 
(MV)  are  also  introduced  to  provide  the  capability  of  performing  memory 
swaps  during  computation  of  a  modulus.  The  updated  list  of  commands  is 
shown  in  Table  4.  This  is  not  an  exhaustive  list  of  the  Level  1 
(primitive)  commands  for  chip  #3022.  Controller  memory  allocation  has 
been  defined  and  performed  under  Task  1  as  part  of  the  set-up  of  the 
Exerciser  unit. 
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Figure  37.  Revised  40-Pin  DIP  Pin  Assignments  for  #3022  Chip 
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OAT  A  IN  I 


III  X  -OFF 


OATA  IN  C 
FLOOD  IN  T 


OATA  IN  X 


DATA  IN  X.B 


I  X  *  ON 


II  B-ON 


OATA  IN  B 


FLOOD  IN  T,  B 
OF  TO  X  -  OFF 


OATA  IN  B.  C 


FLOOD  IN  T,B 
OFTOX  -  OFF 


IV  t  -ON 


IV  C-ON 


TERMINATION 


OATA  IN  C 


IV  B-OFF 


VII  T  -OFF 


Figure  39.  Shift  up  (SU)  Clock  Sequence  with  Termination  Phase  T 
as  Now  implemented 
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TABLE  4.  CCD  MATRIX  PROCESSOR  PRIMITIVES  (CHIP  #3022):  LEVEL  1  COMMANDS 


Function 

Mnemonic  Number  of 

Initial 

Conditions 

States 

Start 

End 

Load/unload  1  stage 

IB,  OB 

17,  17 

Yes 

Yes 

(2  bits  4$ ,  serial  register) 
Load  entire  array 

1A 

4688 

Yes 

Yes 

Unload  entire  array 

OA 

4688 

Yes 

Yes 

Array  shift  right  (»i) 

SR 

21 

Yes 

Yes 

Array  shift  left  (wi) 

SL 

21 

Yes 

Yes 

Array  shift  up  (wi) 

SU 

21 

Yes 

Yes 

Array  shift  down  (Ml) 

— All  shifts:  X  goes  to  X 

SD 

21 

Yes 

Yes 

once  per  pass,  t\  passes 
called 

— SL,  SO:  limited  ri;  SR, 

SU:  unlimited  r) 

Load  store  *1  from  X 

Ll 

Id 

Yes 

Yes 

Load  store  #2  from  X 

L2 

9 

Yes 

Yes 

Load  store  #3  from  X 

L3 

9 

Yes 

Yes 

Unload  store  11  to  X 

U1 

14 

Yes 

Yes 

Unload  store  #2  to  X 

U2 

7 

Yes 

Yes 

Unload  store  #3  to  X 
--"Load"  performs  register 
addition  in  named  store 
--"Unload"  performs  register 

U3 

8 

Yes 

Yes 

addition  with  X 

Erase  store  #1 

El 

7 

Yes 

Yes 

Erase  Store  #2 

E2 

19 

Yes 

Yes 

Erase  store  #3 

E3 

55 

Yes 

Yes 

Rotate  X23 

RF 

14 

Yes 

Yes 

Rotate  X32 

RB 

17 

Yes 

Yes 

Swap  32 

SW 

31 

Yes 

Yes 

Hove  store  #1  to  store  12 

MV 

14 

Yes 

Yes 

Move  Y  to  X 

MX 

17 

HO 

Yes 

Regenerate  store  #1  to  x 

Rl 

23 

Yes 

Yes 

Regenerate  store  *2  to  x 
--original  store  preserved 

R2 

29 

Yes 

Yes 

Difference  ST#1  -  ST#3 

D3 

29 

Yes 

NO 

Difference  ST#1  -  X 

D4 

30 

Yes 

No 

Difference  STll  -  ST# 2 
--D3 ,  D4  :  Result  goes  to 

ST#3>  M2:  result  goes  to  x 

— D3,  D4:  Destroys  ST#1»  M2: 

M2 

68 

Yes 

Yes 

STtl,  ST# 2  conserved 
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The  #3022  architecture,  although  deficient  in  storage,  can  perform  a  simple 
edge  enhancement  algorithm:  the  so-called  Roberts  Cross  operator.  This  is 
a  simple  modulus  sum  on  a  2  x  2  window  as  defined  by: 


R  *  I  a  -  dl  +  I  b  -  d  (negative  results  become  zero) 

The  set  of  primitives  necessary  to  perform  this  is  included  in  Table  4. 

The  sequence  of  major  stages  is  illustrated  in  Figure  40,  and  the  program  is 
shown  below. 

START  Assume  I/O  to  load  array  has  just  occurred,  and  El,  E2,  e3  cleared 
the  stores,  leaving  X  intact. 


I. 

LI 

II. 

Rl,  L3,  Rl,  SU(1),  SL(1),  MV1-2,  Ll 

III. 

03,  generates  sum  SL(1),  ST1  -  ST3 j 

ST1:«  NULL 

IV. 

MV  Y-X,  Ll 

V. 

R2 

VI. 

SU(1),  SL(1) ,  generates  "d"  element 

VII. 

D4,  generates  sum  ST1-X:  ST1:*  NULL 

VIII. 

MV  Y-X 

IX. 

Ll,  El 

clears  "d"  from 

array 

X. 

R2,  SL (1) ,  Ll,  R2,  SU (1) ,  generates 

•b"  and  »c" 

XI. 

D4 

\ 

XII. 

MV  Y-X,  Ll,  R2 

XIII. 

SL(1)  ' 

^  regenerates  "b"  for 

last  difference 

XIV. 

04,  MV  Y-X,  LlT 

El,  generates  R 

This  routine  has  yet  to  be  executed  since  only  fully  functional  #3022  chips 
can  perform  the  new  commands.  XM  is  used  to  provide  temporary  access  to  ST3 
while  charge  is  shifted  along  the  registers. 


THE  WINOOW  IS2x  2  USING  THE  SET 
0.  b.  c.  d}  WITH  a.b  AS  THE  TOP  ROW 
OF  THE  MATRIX. 


DEFINED  QUANTITIES: 

^1  *  (Ka  •  Kd) 

-^2*(Kd-  Ka) 
Sl.Ai+A2-K(1»d1) 
A3  •  (Kb  -  Kc) 
A4-(Ke-Kb) 

R  -A1 +A2  +  A3  +  A4 


Figure  40.  Sequence  of  Primitives  for  43022  to  Generate  Roberts  Cross 
R,  and  Save  Original  Image,  a,  for  Further  Processing 


When  using  the  #3022  chip  with  the  correct  termination,  T,  the  extra  clock 
drive  must  be  coded  as  appropriate.  This  nas  been  performed  also  as  a  part  of 
Task  1.  A  photograph  of  the  entire  chip  is  shown  in  Figure  41. 


Processing  runs  (3022-1,  3022-2,  3022-3)  were  executed  through  November  1981 
and  resulted  in  discovery  of  some  mask  errors  missed  on  the  redesign.  This 
resulted  in  failure  of  the  first  and  second  wafer  runs  (12  wafers  each) .  The 
third  process  run  was  successfully  completed  and  detailed  measurements  made  at 
the  probe  card  level  to  select  parts  for  packaging. 

Since  the  first  two  wafer  runs  were  lost,  two  further  wafer  runs  were  started 
(3022-4,  3022-5)  at  the  end  of  November  1981.  Run  4  was  completed  and 
received  for  testing  with  some  parts  selected  for  packaging.  At  the  time  of 
writing  Run  5  remains  to  be  tested. 

Packaged  parts  from  both  Run  3  and  Run  4  have  been  delivered  to  RADC  with 
further  parts  to  be  sent  as  they  become  available.  It  is  expected  that 
quality  and  performance  of  later  parts  will  be  better  than  the  first  samples. 
Honeywell  is  anxious  to  provide  as  good  a  selection  of  parts  as  possible. 

The  results  of  Run  3  parts  measurements  are  described  in  the  next  subsection. 

TASK  4— FINAL  16  x  16  CCD  CHARACTERIZATION 

Chip  Tests  and  Wafer  Tests  of  #3022  CCDs 

Wafer  testing  has  been  performed  on  four  process  runs  with  good  parts  obtained 
from  the  third  and  fourth  runs.  Further  process  runs  will  be  characterized 
later  on  Honeywell's  internal  programs. 


The  results  from  run  3  wafers,  some  parts  of  which  have  been  packaged  and 
delivered  to  RADC,  are  shown  in  Figures  42  through  58.  The  elementary 
functions  of  shifting  data  through  the  external  four-phase  shift  register  only 
and  then  also  through  the  three-phase  matrix  shift  registers  were  used  as  the 
"go/no  go"  test  criteria. 

yields  of  > 90%  for  the  external  four-phase  registers  were  observed  on  all 
wafers  (these  registers  have  the  least  complex  interconnect).  The  internal 
matrix  registers  showed  about  40%  yield  with  the  major  failure  mode  apparently 
due  to  open-circuit  clock  phases  with  a  random  fault  distribution.  Of  those 
devices  that  passed  this  second  level  of  tests,  there  were  almost  100%  passes 
to  the  next  level  of  testing.  This  next  level  operates  not  only  the  external 
registers  and  matrix  shift  registers,  but  it  also  loads  and  unloads  each  of 
the  stores. 

Figures  42  and  43  show  array  shifting  with  zero  data  for  two  cases:  (1)  using 
"Erase"  primitives  and  "Swap  2«-^3*  primitive;  (2)  using  only  input  array 
(IA)  and  output  array  (OA)  subroutines.  The  primitives  used  are  shown  on  the 
TTY  print-out  in  each  figure.  These  measurements  were  executed  using  the  RADC 
Exerciser  unit  adjusted  manually  to  give  a  slow  frame  time  of  about  47  msec  in 
order  to  show  leakage  charge  effects  on  these  packaged  parts. 

Satisfactory  performance  is  evident  in  these  photographs — the  leakage  charge 
ramps  slightly  from  right  to  left  per  line  of  16  bits  due  to  build-up  during 
the  IA  subroutine  (expected  and  normal) ,  and  ramps  very  slightly  from  left  to 
right,  per  whole  frame  of  256  bits  during  OA  subroutine  (expected  and 
normal) .  It  should  be  noted  that  a  "corner  turn"  is  executed  by  IA,  OA 
subroutines  by  choice — a  different  choice  with  no  corner  turn  is  available  if 
required,  but  the  subroutines  have  not  been  included  in  the  EPROM  set. 
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42.  Leakage  Current  with  Figure  43.  Leakage  Current  with  Zero 

Zero  Data  (25°C,  Data,  without  Erasing  stores 

47  msec) 


The  same  program  was  executed  but  with  a  6  x  6  test  pattern  in  Figures  44  and 
45.  Again,  satisfactory  performance  has  been  obtained.  Note  also  that 
transfer  efficiency  on  this  device  (3022-3-8  chip  #64)  is  good.  An  estimated 
ffe  "  0.1  is  evident,  where 

N~  1/2(16  x  8)  +  1/2(16  x  10)  +  1/2(16  x  10) 

+  1/2(16  x  8)  +  4  -  292  transfers 

-4 

This  suggests  a  mean  transfer  inefficiency  of  3  x  10  (x  »  0.9997),  which 
is  acceptable  performance  for  a  BCCD.  A  slightly  better  transfer  efficiency 
is  expected  from  the  best  devices  which  should  be  as  good  as  the  original 
"best  performance"  #2214  parts  used  earlier. 

Loading  and  unloading  the  three  stores  is  shown  in  Figures  46  through  48.  The 
"swap  memory"  routine  SW  was  tested  in  Figure  49,  with  partial  success.  The 
cause  for  the  degradation  in  the  case  where  L2  is  executed  followed  by  SW  and 
U3  is  as  yet  unexplained.  A  software  error  may  be  present:  it  is  suspected 
because  the  other  choice  (L3,  SW,  U2)  is  quite  satisfactory.  There  is  a 
significant  difference  in  charge  flow  with  the  L2,  SW,  U3  choice  requiring  all 
controls  to  be  correct,  whereas  L3,  SW,  U2  is  a  simpler  use.  This  is  still 
under  investigation  with  more  testing  of  run  4  parts. 

Leakage  charge  collected  by  the  stores  during  IA  +  OA  frame  time  is  unloaded 
and  read  out  for  the  photographs  in  Figures  50  through  52.  There  is  no  major 
difference  between  the  three  stores  and  an  acceptably  low  leakage  is 
attained:  about  10%  full  well  in  the  47  msec  IA  +  OA  frame  time.  Evidently 
the  signature  of  ST1  on  #2214  chips  has  been  successfully  eliminated  on  the 
#3022  chip. 
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ENTER  PRIMITIVES 


ENTER  PRIMITIVES 


Figure  44.  Performance  with  Data 
(70%  full  well,  0% 
fat  zero)  Showing 
Transfer  Efficiency 


Figure  45.  Performance  with  Data, 
Without  Erasing  Stores 
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Figure  46.  Loading  and  Unloading  STl 


Figure  47.  Loading  and  Unloading  ST2 
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ENTER  PRIMITIVES 


Figure  49.  Swapping  Data  (ST2 ,  ST3)j  Software  Error  Suspected 
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The  most  crucial  test  is  that  of  using  the  FGA  as  a  regenerator.  The  output 
from  the  array  after  Rl  has  been  executed  is  shown  in  Figure  53.  The 
Exerciser  was  set  to  operate  the  CCD  at  a  faster  rate  to  reduce  leakage 
charge  effects.  As  a  calibration,  the  subroutine  was  extended  to  store  this 
regenerated  output  in  ST2  and  then  unload  STl.  This  then  shows  the  charge 
that  was  loaded  into  the  FGA  input  Z.  Comparing  the  two  outputs  in  detail 
(see  Figures  53  and  54),  it  can  be  seen  that: 

1.  Relatively  good  uniformity  is  obtained  apart  from  a  baseline  offset 
nonuniformity. 

2.  The  charge-charge  unity  gain  is  approximately  attained,  as  desired. 

3.  There  is  no  discernible  1/f  noise-smearing  of  the  output  pulses. 

This  is  a  rather  pleasing  result  that  validates  the  design  modifications 
introduced  to  the  3022  chip.  A  further  test  on  the  regenerator,  using  the 

Rl  primitive  again,  is  shown  in  Figure  55.  The  6x6  test  pattern  has  been 

set  (by  TTY  command  *E  and  data  location/value)  to  give  a  linear  “staircase" 

in  place  of  the  simple  "box*  of  unit  amplitude  used  previously.  The  data  is 

read  out  in  the  same  manner  as  just  described. 

Notice  that 

4.  A  linear  response  has  been  obtained  up  to  a  limiting 
amplitude — offset  effects  should  be  considered  separately. 

This  validates  the  other  design  changes  to  the  FGA  and  is  as  computed 
from  the  model  given  in  Appendix  A. 

The  multiplier  structure  is  the  remaining  feature  to  be  tested.  The 
operation  of  the  multiplier  (see  Figure  56)  is  briefly  a:  follows.  The 
FGA  control,  RSNZ,  is  assumed  to  have  been  at  OFF  and  then  switched  ON 
ready  to  "float"  the  Z  and  Y  gates  by  switching  off  FGR.  At  this  point, 

MG 2  must  be  set  to  give  the  desired  loading  factor  on  the  FGA  input,  Z. 
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Figure  53.  Regenerator  Performance  Figure  54.  Regenerator  Performance 
--Output  Charge  --Input  Charge 
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Figure  55.  Regenerator  Performance  Showing  Linearity 
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Figure  56.  Physical  Basis  for  Multiplier  Gain  Control  (Patent  applied 
fors  Honeywell  #a4109674-US) 
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Having  set  MG2  to  an  analog  bias,  MD  is  then  switched  from  its  default 
OFF  (OV)  to  its  active  ON  (+15V)  state.  This  allows  excess  channel 
charge  from  under  MG2  to  flow  out  to  the  MD  diode  sink.  The  surface 
potential  under  the  multiplier  gate  (the  large  Poly  I  gate  visible  in 
Figure  34)  is  thereby  set  to  the  desired  value  and  the  z  loading  factor 
is  fixed.  FGR  can  be  switched  OFF  next  and  then,  in  a  subsequent  clock 
state,  signal  charge  can  be  loaded  into  the  Z  floating  gate  as  usual. 

We  have  tested  the  Rl  primitive  with  MD  inactive  at  default  •  +3.5V  and 
MG2  adjusted  to  give  a  "low"  gain  and  then  a  "high"  gain  manually  to 
observe  operation. 

5.  The  gain  control  is  completely  effective. 

Figures  57  and  58  show  the  output  from  an  Ll  Rl  subroutine  with  the  DC 

setting  of  MG2  and  MD  shown  as  the  lower  two  oscilloscope  traces.  MG2  is 

set  high  in  Figure  57  and  low  in  Figure  58. 

Further  testing  of  this  primary  function  is  in  progress  to  observe 
dynamic  gain  setting.  No  data  is  available  at  the  time  of  writing.  When 
data  is  available,  however,  the  results  will  be  sent  to  RADC  through  the 
present  contract  monitor  (Mr.  Paul  Pellegrini). 

The  measurements  show  that  all  the  basic  components  of  the  3022  cell 
function  as  designed.  It  remains  to  be  seen,  however,  where  software, 
drive  waveform,  and  fabrication  defects  limit  the  performance  of  the  3022 
chips  in  executing  a  complete  threshold  and  edge  extract  algorithm.  SEM 
photographs  of  the  3022  cell  and  a  couple  of  the  contact  areas  are  shown 
in  Figures  59  and  60.  The  diagonal  edge  of  the  Z  FGA  input  gate  is 
visible  at  bottom  right  in  Figure  59.  Figure  61  shows  the  connection  of 

RSNI  reset  transistor  source  diffusion  to  the  Poly  I  gate,  which  is  the 

multiplier-to-FGA  noninverting  input  gate  connection.  Figure  60  shows 
the  metal  track  for  the  ST3  bus  directly  above  the  RSNI  contact  on  the 
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Figure  57.  Multiplier  Performance  Figure  58.  Multiplier  performance 
(MG2  High)  (MG2  Low> 
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Figure  59.  SEM  Photograph  of  3022  Cell 

third  track  above  it.  Notice  that  step  coverage  and  alignment  are 
excellent.  A  "spike  inhibit"  undoped  poly  layer  is  applied  after  contact 
windows  are  etched  and  before  metal  is  applied  (Cu-doped  Al) .  It  is 
possible  that  the  low  temperature  (450°C)  sinter  for  making  the  ohmic 
contacts  to  the  gates  is  not  completely  effective  and  that  this  may  be 
the  cause  of  the  low  yield  of  functional  chips.  This  is  under 
investigation. 

System  Testa  with  Original  >2214  CCDs 

During  the  course  of  the  RADC  contract  execution  a  related  contract  from 
Eglin  AFB,  *F08635-80-C-0246  to  Honeywell's  Avionics  Division,  was 
provided  with  samples  of  Honeywell's  #2214  CCDs  together  with  operation 
software  for  all  primitives  plus  simple  subroutines.  The  subroutines 
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were  developed  partly  on  Honeywell  funding  and  partly  on  this  present 
RADC  contract.  At  the  direction  of  RADC  (Dr.  Freeman  Shepherd)  the 
necessary  information  to  operate  the  devices  was  supplied  to  the  Eglin 
contract  to  coordinate  the  two  programs  and  enhance  the  yield  of 
results.  There  was  no  implication  that  fully  functional  parts  were 
available,  however.  It  was  already  known  that  significant  deficiencies, 
which  were  being  corrected  in  the  present  development  program,  existed  in 
the  #2214  chip  design. 

A  somewhat  impressive  image  processing  performance  was  nonetheless 
attained  by  approximately  "tuning"  the  2214  drive  waveforms.  The 
limiting  factors  on  the  chip  were  as  already  described  in  detail  in  Task 
2.  By  using  the  Eglin  DCIU  controller,  however,  it  became  possible  to 
mask  out  border  degradation  from  the  resulting  image  and  also  to  limit 
the  dynamic  range  of  the  final  display  to  mask  the  nonuniformities 
generated  by  the  2214  chip  defects.  The  result  on  the  simple  threshold 
and  partial  edge  algorithms  is  shown  using  FLIR  imagery  (source: 
NVfrEOL/Honeywell-IRL  frame  #01926)  in  Figures  62  and  63.  A  complete  set 
of  photographs  has  been  supplied  to  RADC  showing  various  test  conditions. 

In  Figure  62  the  Tracking  Gate  mode  has  been  invoked  and  a  selected 
target  located  manually  through  the  TTY.  In  Figure  63  the  upper  left 
quadrant  (#0)  displays  a  "thresholded”  version  of  the  target  image  using 
the  subroutine: 

IA  El  E2  S3  LI  M2  SL  Ll  M2  SL  Ll  Ul 

This  enhanced  contrast  image  (stored  in  RAM)  is  then  further  operated  on  by 
the  2214  chip  to  extract  the  bottom/right  edges  shown  in  the  upper  right 
quadrant  (#1)  using  the  subroutine: 

IA  El  E2  E3  Ll  Rl  L3  R1  L3  Ul  SU  SL  Ll  Rl  L2  Rl  L2  US  El  Ll  M2 
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Figure  62.  FLIR  Image  To  Be  Processed  by  2214  Chip 


Figure  63.  Processed  Images  using  Threshold  and  Subtract  Subroutines 


Similarly,  the  upper/left  edges  are  extracted  and  displayed  in  the  bottom  left 
quadrant  (#2)  using  the  subroutine: 

IA  El  E2  E3  LI  Rl  L3  Rl  L3  U1  SD  SR  Ll  Rl  L2  Rl  L2  U3  El  Ll  M2 


Finally,  the  #1  and  #2  quadrant  images  are  added  by  the  microprocessor  in  the 
controller  to  produce  the  result  shown  in  #3. 

The  results  are  remarkable  considering  the  low  performance  of  the  2214  chip. 
They  suggest  that  in  practice  significant  capability  can  be  obtained  without 
the  requirement  for  "perfect"  CCD  performance.  The  algorithms  provide  an 
image  processing  gain  of  S/N,  which  can  in  some  cases  tolerate  considerable 
non idealities. 

The  performance  of  a  functioning  3022  chip  is  expected  to  greatly  exceed  that 
of  the  2214,  however.  Not  only  should  the  analog  arithmetic  be  dramatically 
better,  but  the  3022  will  provide  in  one  subroutine  the  complete  threshold  and 
Roberts  cross  full  edge  extractor  processed  images.  The  3022  will  provide  as 
output  options: 

a.  Original  image  only 

b.  Thresholded  image  only 

c.  Edge  image  (2x2  Roberts)  only 

d.  Sum  of  (a)  and  (c) :  edge-enhanced  original  image 

e.  Sum  of  (b)  and  (c) :  edge-enhanced  thresholded  image 

And  the  3022  will  do  these  entirely  without  the  need  to  use  the  microprocessor 
of  the  controller  in  the  RADC  Exerciser  to  do  partial  sums. 


w  -v 


L 


The  subroutine  to  perform  (c)  is: 


IA 

El 

E2 

SW 

E2 

Ll 

Rl 

L3 

Rl 

SU 

SL 

MV 

D3 

MX  LI 
R2 

SU  SL 

D4 

MX  SD  SD 
Ll  El 

R2  SL  Ll  R2  Su 
D4 

R2  SL 

D4  MX  Ll  El 

Roberts  Cross,  R,  placed  in  position  "a"  in  ST3. 

To  continue  to  produce  (d)  uses: 

all  of  the  above,  plus 
U3  U1 


This  is  a  full  real-world  algorithm,  and  is  the  maximum  window  size  on  which 
the  3022  can  perform  edge  extraction.  Notice  that  this  is  a  full  set  of  two 
"modulus"  operations.  The  3x3  Sobel  Edge  is  a  "better"  algorithm  but  must 
wait  for  the  second-generation  chip  architecture  before  a  demonstration  is 
possible  since  this  architecture  requires  a  fourth  storage  site. 

TASK  5 — SYSTEM  ARCHITECTURE  AND  SOFTWARE  DEFINITION 

In  this  subsection  we  will  assess  the  architectural  strengths  and  weaknesses 
of  PIP  in  performing  image  processing  algorithms.  An  appropriate  starting 
point  is  the  discussion  of  the  required  classes  of  operations  in  target 
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acquisition  and  tracking.  (The  advantages  of  a  parallel  architecture  will 
become  readily  apparent.)  We  will  then  discuss  implementation  of  each  of  the 
algorithm  classes  above  in  the  PIP  by  using  specific  examples.  We  felt  that 
analysis  of  a  system  example  using  a  series  of  algorithms  would  be  an 
appropriate  way  to  provide  the  inputs  needed  for  the  new  cell  design  in  Task 
6.  Toward  that  end  we  picked  a  chain  of  algorithms  very  similar  to  those  used 
on  an  existing  Honeywell  seeker  program  and  obtained  a  set  of  specific 
improvements  incorporated  in  the  new  cell.  Additional  input  for  the  new 
design  came  from  modeling  of  the  existing  PIP  mnemonic  set;  that  will  also  be 
described  in  this  subsection. 


In  this  subsection  we  also  outline  the  basic  form  of  PIP  software  and  discuss 
some  strategies  possible  in  PIP  control. 

Target  Acquisition  and  Tracking  Algorithms 

Acquisition  Subsystems — Acquisition  encompasses  a  number  of  approaches, 
ranging  from  simple  matched  filters  to  advanced  target  screening/cueing 
subsystems.  The  simple  algorithms  closely  parallel  the  tracking  algorithms 
themselves  (like  the  matched  filters),  and  are  effective  only  in  special 
circumstances  that  still  require  man-in-the-loop  interaction  and  are 
ineffective  in  high  clutter  environments.  These  are  typical  of  acquisition 
schemes  found  in  current  trackers.  They  are  generally  designed  to  reacquire  a 
target  following  loss  of  track,  or  to  lock  on  a  target  in  the  window 
positioned  approximately  in  the  target  trajectory  by  an  operator.  Automatic 
target  screening/cueing  approaches,  on  the  other  hand,  are  more  effective  for 
target  acquisition,  although  they  were  originally  developed  for  operator 


cueing.  Their  performances  in  bigb-clutter,  low  signal/noise  environments 
are  superior  to  the  simpler  target  acquisition  schemes  found  in  today's 
trackers.  Some  of  these  schemes  can  also  recognize  targets;  this  is  an 
asset  because  it  enables  the  prioritization  of  targets,  an  essential 
function  in  autonomous  fire  control. 

Automatic  target  acquisition  and  recognition  architectures  and  hardware  have 
been  developed  by  several  contractors,  including  Honeywell  (Prototype 
Automatic  Target  Screener  built  for  the  Army  NV6E0L  and  the  imaging  Sensor 
Autoprocessor  being  designed  for  AFAL) ,  Rockwell  (the  ISA  above) ,  Hughes 
(the  Autocuer  for  the  ATAC  FLIR) ,  Texas  instruments,  and  westingbouse. 
Although  not  all  of  these  are  in  hardware  form,  the  functional  approac  i 
are  similar.  All  employ  a  sequence  of  successive  subfunctions  shown  i  ...e 
block  diagram  of  Figure  64. 

The  video  image  is  first  segmented  into  smaller  regions  based  on  intensity 
thresholds,  edges,  or  texture.  These  regions  are  then  associated  into 
objects  or  object  components.  More  measurements  (features)  are  computed  on 
these  extracted  objects.  The  features  are  then  used  to  classify  the  objects 
as  targets/nontargets  and  into  a  target  type.  Target  recognition  has  the 
added  benefit  of  being  able  to  prioritize  the  object  to  be  tracked. 

Numerous  algorithmic  variations  in  the  subfunctions  are  possible,  of 
course.  The  algorithms  do,  however,  fall  into  seven  classes: 

1.  Smoothing  Algorithms — Linear  boxcar  filters,  3x3  and  5x5 
median  filters 

2.  Edge  Algorithms— Sobel,  Laplacian,  gradient,  Roberts  Cross, 
directional  masks,  etc. 

3.  Thresholding  Schemes— Auto  threshold.  Superslice, 
two-dimensional  histograms 


Figure  64.  Functional  Block  Diagram  of  Generic  Automatic  Target 
Screening/Cueing  Approaches 


4.  Association  of  Classified  Pixels— Object  growing,  connected 
components,  edge/perimeter  match 

5.  Edge  Thinning— Thickening,  nonmaximum  suppression,  perimeter 
extraction,  binary  noise  removal 

6.  Feature  Extraction  Functions— Moment  computation,  perimeter 
extraction  area,  texture,  average  and  peak  intensity 

7.  Target  Classification  Algorithms— Linear  discrimination,  nearest 
neighbor  classifiers,  etc. 

In  later  subsections  we  will  examine  specific  implementation  in  PIP  of 
representative  algorithms  from  each  class. 

Tracking  Algorithms — Commonly  used  tracking  algorithms  can  be  categorized 
as: 

o  Centroid 

— Centroid  of  brightness 
— Centroid  of  area 
—Area  balance 

o  Edge  and  Generalized  Edge 
— Adaptive  gate 

— Leading  and  trailing  edge  averages 

o  correlation 

— Mean  absolute  difference 
— product  correlator 
—Gradient  correlators 

o  Feature  Tracker 

— Feature  matching  (Honeywell  concept) 

o  Matched  Filter 
— Fixed  filter 
— Adaptive  filter 
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All  operations  such  as  convolution,  correlation,  minimum  absolute  distance 
correlation,  matched  filter,  edge  operators,  etc.,  fall  into  the  arithmetic 
category.  The  basic  arithmetic  functions  are  similar  to  the  point 
operations  but  involve  more  variables  and  include  add,  subtract,  multiply, 
and  absolute  value. 

Sliding  window  operations  involving  binary  logic  include  operations  on 
binary  images — usually  in  acquisition  algorithms  for  connectivity  analysis, 
noise  point  removal,  feature  extraction,  and  perimeter  extraction.  In  these 
functions,  every  resultant  image  point  is  a  fixed  boolean  function  of  the 
original  image  points  in  a  fixed  neighborhood. 

Sliding  window  operations  requiring  compare  (or  sort)  operations  include  the 
median  filter,  local  peak  detector  to  find  the  maximum  of  the  correlation 
surface,  and  nonmaximum  suppressor  operators  for  edge  thinning. 

Thus  far  we  have  talked  only  in  general  terms  about  the  types  of  operations 
required  for  acquisition  and  tracking.  A  system  solution,  of  course, 
requires  implementation  of  the  algorithms.  That  is  the  subject  of  the  next 
subsection. 

Processor  Architectures  for  Signal  and  image  processing 

There  are  two  generic  approaches  to  real-time  image  processing  hardware: 
serial  and  parallel.  Present  approaches  are  overwhelmingly  serial  for  two 
reasons.  First,  sensor  output  formats  have  been  generally  serial;  despite 
the  two-dimensional  nature  of  image  data,  the  general  approach  has  been 
serial  image  scanning.  Even  current  staring  arrays  provide  a  multiplexed, 
serial  video  data  stream.  Second,  the  vast  majority  of  processing 
algorithms  have  been  developed  and  simulated  on  serial  computers.  Even 
parallel  operators  like  M  x  M  sliding  windows  and  convolutions  have 

2 

generally  been  implemented  in  serial  fashion  with  M  -  1  FIFOs  and  (M  -  1) 
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single-pixel  delays.  The  weighted  sub  of  all  M  pixels  in  the  window  must 
be  computed  at  video  rates.  This  is  often  a  difficult,  power-consuaing  task. 

Serial  architectures  have  several  drawbacks;  they  are  long  on  hardware  and 
short  on  speed  and  flexibility.  Even  hardwired  serial  modules  often  cannot 
operate  at  875-line  video  rates  and  when  algorithms  change,  either  the 
entire  module  structure  must  be  changed  or  a  new  nodule  developed. 

The  desired  architecture  for  image  processing  must  be  flexible  and  yet 
handle  the  high  throughputs  needed  for  image  processing.  Parallel 
processing,  however,  can  handle  high  throughputs  even  when  using  low-power, 
slower  technology,  a  dramatic  example  is  the  human  visual  system,  which 
does  parallel  processing  orders  of  magnitude  faster  than  the  fastest  serial 
computer  despite  a  basic  processing  element  (the  cortex  neuron)  much  slower 
than  standard  TTL. 

In  addition,  many  of  the  image  operations  we  have  discussed  are  local  and 
parallel  in  nature;  they  can  be  decomposed  into  simple  sequences  of  fast 
parallel  primitives,  as  we  shall  shortly  see.  Unfortunately,  digital 
implementations  of  fully  parallel  approaches  have  not  yet  been  successful. 

The  pip  architecture,  however,  is  both  flexible  and  fast,  can  potentially  be 
mated  directly  with  analog  sensors,  is  compact,  and  consumes  very  little 
power.  The  features  of  the  PIP  chip  were  described  earlier,  but  we  will  now 
discuss  them  from  an  architectural  perspective. 


An  interesting  point  to  note  is  that  parallelism  has  long  been  associated 
with  only  the  "high-throughput"  functions  in  Figure  64— data  conditioning 
and  segmentation.  We  can,  however,  also  cast  other  functions  in  parallel 
form. 

For  example,  let  us  examine  the  processing  sequence  for  a  perimeter 
extraction.  We  have  an  object  consisting  of  binary  "ones"  and  wish  to 
eliminate  superfluous  "ones"  to  form  just  a  boundary. 
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A  simple  operation  for  doing  this  is  merely: 


F  - 

F  •  (B  • 

C  •  D  •  G  •  K  •  J  •  H  •  E) 

parallel  processor 

does  this  as  follows: 

Operation 

Motes 

1. 

SD  (1) 

Bring  C 

over 

the 

center  (m,n) 

2. 

A  A 

A  «  1 

3. 

A  *  I  -*•  A 

A  -  C 

4. 

SL  (1) 

Bring  D 

over 

the 

center 

5. 

A  •  I  -*■  A 

A  »  C  • 

D 

• 

17. 

A  •  I  -*■ 

A 

A  -  B  •  C»  D •  G» 

K  •  J  •  B  •  E 

18. 

A  A 

A  ■  B  •  C*  D  *  G* 

K  •  J  •  H  •  E 

19. 

SR  (1) 

Bring  F  over  (m,n) 

20. 

A  •  I  ♦ 

A 
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We  require  only  20  operations 


Matrix  Processor  Control 


The  rationale  for  parallel  processing  and  the  type  of  instructions  our  CCD 
matrix  processor  can  execute  have  been  summarized.  Use  of  the  device, 
however,  requires  that  we  discuss  control,  in  this  subsection  we  will  outline 
the  control  philosophy  of  the  processor  and  provide  details  on  the  specific 
set  of  cossaands  that  have  been  used  throughout  this  program. 

The  CCD  matrix  processor  operates  from  a  number  of  control  lines,  as  do  all 
CCDs.  These  control  lines  provide  either  DC  levels  or  clocking  waveforms  that 
generally  assume  a  high  or  low  state.  (A  few  lines  could  require  a  continuum 
of  signals;  these  will  be  treated  shortly.)  In  the  case  of  the  2214  we  must 
provide  21  clocked  waveforms.  One  set  of  possible  values  of  these  clocks  can 
be  considered  a  clock  state;  we  can  consider  each  state  as  a  digital  word  in 
which  each  bit  is  the  high  or  low  state  of  a  particular  clock  line.  Control 
of  the  matrix  processor  is  then  effected  by  a  series  of  digital  words  that 
specify  system  states;  that  is,  a  finite-state  machine.  (The  digital  words 
must  be  shifted  and  buffered  for  compatibility  with  the  CCD,  but  that  is  of  no 
consequence  from  a  control  standpoint.)  Specifying  each  control  state  is  a 
time-consuming  task,  since  detailed  manipulation  of  the  charge  packets 
requires  device  knowledge. 

We  can,  however,  concatenate  these  digital  words  into  short  routines  that  are 
used  repeatedly  and  are  themselves  concatenated  to  form  algorithms.  This 
first  level  of  abstraction  will  be  called  a  "primitive.”  (As  an  example,  we 
now  can  request  a  "store  1"  command,  or  ST1,  which  takes  data  stored  under  the 
X  node  and  transfers  that  charge  packet  for  storage  under  STl.)  We  have  now 
freed  the  software  developer  from  one  responsibility.  He  need  only  know  the 
arithmetic  capability  of  the  cell,  not  how  the  cell  does  the  desired 
operation,  coding  of  the  processor  can  be  considerably  more  rapid.  A  higher 
level  of  abstraction  would  refer  to  blocks  of  matrix  processor  primitives  as 
macros;  an  example  of  this  is  given  later  in  this  subsection. 


in  choosing  a  primitive  set  for  the  2214  chip,  we  considered  a  number  of 
factors.  The  basic  tradeoff  is  between  large  memory  requirements,  a  large 
number  of  primitives,  and  high  speed  versus  smaller  storage  requirements,  a 
core  subset  of  primitives,  and  lower  speed.  We  opted  for  the  smaller  set,  but 
this  tradeoff  should  be  made  in  a  more  systematic  fashion  when  the  new  cell  is 
fabricated.  Table  5  lists  all  the  original  matrix  processor  mnemonics  and 
their  functions;  several  have  been  and  are  being  added  to  this  list. 

Appendix  C  provides  a  complete  listing  of  all  clock  states  for  each  of  these 
mnemonics. 

Before  proceeding  to  an  example  in  the  next  subsection,  we  will  discuss  some 
additional  control  issues  that  are  not  yet  resolved.  All  of  the  current 
primitives  require  that  the  start  and  end  of  each  primitive  correspond  to  the 
same  state;  that  is,  there  is  no  restriction  in  concatenation  of  primitives. 
While  this  simplifies  programming,  it  may  sacrifice  speed,  because  states  that 
are  useless  in  a  particular  sequence  have  been  added  for  compatibility  among 
primitives.  Alternatives  to  the  current  approach  might  be  explored. 

A  similar  issue  arises  in  control  of  the  racetrack  storage  loop  (part  of  the 
new  cell  that  will  be  described  later).  The  simplest  control  structure  would, 
in  effect,  label  each  storage  site.  (At  the  beginning  and  end  of  any 
instruction  that  accessed  the  loop  we  would  have  charge  packets  in  the  same 
relative  locations.)  This  may  lose  speed,  however,  since  we  might  again  make 
meaningless  shifts  just  to  get  to  initial  and  final  conditions  prior  to  using 
the  same  storage  locations  again.  Worse  yet,  we  would  not  be  able  to  use  the 
loop  in  a  stack  architecture  mode  by  shifting  the  stack  one  location  at  a  time 
either  forward  or  backward.  The  penalty  for  such  versatility  is  more  complex 
control  logic  that  must  keep  track  of  the  storage  loop  position  before  each 
instruction.  This  issue  requires  closer  examination. 


TABLE  5.  2214  MATRIX  PROCESSOR  MNEMONICS  FOR  FUNCTION  IMPLEMENTATION 


Mnemonic 

Function 

Execution  Time 
(clock  periods) 

I/O 

READ  IN,  SHIFT,  READ  OUT 

17 

SR 

SHIFT  RIGHT  X  -*■  X 

21 

SL 

SHIFT  LEFT  X  *  X 

21 

su 

SHIFT  UP  X+  X 

21 

so 

SHIFT  DN  X*  X 

21 

LI 

LOAD  MEMORY  1  X  -*•  1 

18 

L2 

LOAD  MEMORY  2  X-*-  2 

9 

L3 

LOAD  MEMORY  3  X  ♦  3 

9 

R1 

REGENERATE  (Sj^  •  KjJCOPY  1  to  X 

22 

R2 

REGENERATE  (S2  •  KjJCOPY  2  to  X 

26 

U1 

UNLOAD  MEMORY  1  -*•  X 

14 

U2 

UNLOAD  MEMORY  2 -*•  X 

7 

U3 

UNLOAD  MEMORY  3  -*>  X 

8 

El 

ERASE  MEMORY  1 

7 

E2 

ERASE  MEMORY  2 

19 

E3 

ERASE  MEMORY  3 

55 

M2 

SUBTRACT  2  FROM  1 

67 

MOD 

MODULUS  1  2  -  11  *  (2  -  1)  +  (1  -  2)+  3 

138 

A  final  control  issue  is  again  more  relevant  to  tbe  new  cell.  In  tbe  2214 
many  clock  lines  serve  multiple  functions.  This  reduces  tbe  number  of  control 
lines  but  also  practically  eliminates  tbe  possibility  of  overlapping 
instructions;  that  is,  starting  a  second  instruction  before  its  predecessor 
is  complete,  in  the  new  matrix  processor  cell  design  tbe  control  lines  are 
more  independent;  tbis  affords  tbe  possibility  of  a  control  architecture 
where  two  instructions  are  executed  slightly  faster  than  in  a  conventional 
single-stack-pointer  architecture.  Again,  more  work  is  needed  here. 


PIP  Architecture 


Tbe  2214  chip  is  a  fully  parallel  processor.  Each  arithmetic  unit  cell,  under 
external  control,  performs  exactly  the  same  sequence  of  operations  as  all 
other  unit  cells.  Therefore,  we  can  examine  tbe  sequence  of  operations  at  a 
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single  location  and  realize  that  256  such  sequences  are  simultaneously 
occurring.  We  will  refer  to  quantities  like  A(i,j)  and  imply  a  16  x  16 
array.  The  cell  contains  three  dedicated  storage  locations  (STl(i,j), 

ST2(i, j),  and  ST3(i,j))  which  can  also  serve  as  accumulators.  We  can  perform 
shift  operations  from  cell  (i,j)  to  cell  (i+m,  j+n) j  if  we  define  the  corner 
location  as  X(i,j),  this  implies  that  X(i+m,  j+n) +  X(i,j).  This  notation 
implies  that  X(i,j)  is  the  destination  register  for  a  PIP  result.  The  heart 
of  the  cell  is  the  regenerator  whose  output  is  given  by: 

KlZ (i»  j)  -  KlY ( i ,  j )  -*■  R(i,  j) 

(Charge  can  be  moved  from  any  site  to  Z  or  y.) 

Numerous  operations  (including  subtraction,  initialization,  multiplication, 
and  absolute  value)  can  be  derived  from  this  useful  function.  The  arithmetic 
capabilities  can  be  summarized  as  follows  (STX(i,j)  represents  any  of  the 
three  stores) : 


o 

Add 

STX ( i , j )  +  I (i, j)  ♦ 

STX ( i , j ) 

o 

"Subtract" 

f  K, (Z(i, j)  -  Y(i, j) ) 

if  Z  -  Y> 

0-*"  R(i,j) 

lo 

if  Z  -  Y< 

0-*-  R(i,j) 

o 

Multiply  by  a  constant 

K^Z  (i, j)  -►  R(i,j) 

o 

Clear 

0  -*•  STX  ( i ,  j ) 

o 

Initialization 

a  *►  R(i, j) 

o 

Absolute  difference 

X(i, j)  -  I  (i, j)  ♦  X ( i, j ) 

o 

Store 

I  (i,  j)  STX ( i ,  j ) 

The  modified  subtraction  performed  by  PIP  is  very  useful  in  practice;  the 
algorithm  example  presented  later  will  explain  why.  In  addition  to  the 
arithmetic  operations  above,  certain  binary  operations  can  be  performed  in 
the  cell  by  changing  the  usual  operation  of  the  regenerator.  Binarization 
of  an  image,  for  example,  is  done  as  follows: 
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o 


Define  binarization  as  f  1  if  ((i«j)  >  a  -►  X(i,j) 

l  0  if  l(i,j)  <  a  ♦  X(i,j) 


o  Generate  the  local  threshold  a  as  shown  in  initialization  above 

o  perform  the  single-sided  subtraction  shown  above 

o  Perform  analog  addition  between  the  result  of  the  subtraction  and 
itself  N  times  to  increase  the  gain. 

o  Use  charge  skimming  on  that  result  using  an  appropriate  MOS  gate 
within  the  PIP  cell  as  a  reference  surface  potential  and  an 
adjacent  MOS  gate  as  the  store  of  the  input  charge.  The  excess 
charge  spills  over  the  reference  barrier  and  is  either  stored  in  a 
third  gate  or  sent  directly  to  a  sink  diode.  An  external  D/A  can 
be  used  to  control  the  reference  level  of  the  chosen  clock  phase. 
(The  D  phase,  for  example,  could  be  used  for  this.) 

This  function  greatly  enhances  the  ability  of  PIP  to  perform  several  of  the 
classes  of  image  processing  operations  mentioned  earlier.  Certain  other 
operations  require  AND,  OR,  and  NOT  logical  operations  and/or  comparison 
operations. 


Logical  operations  are  somewhat  difficult  to  implement  in  direct  PIP  form, 
although  following  binarization  of  two  arrays  (as  indicated  above)  we  could 
use  the  regenerator  to  perform  certain  logical  functions.  For  example,  the 
XOR  function  is  achieved  by  matching  the  gain  on  the  inverting  and 
noninverting  inputs  and  performing  the  absolute  difference  operations.  The 
AND  function  is  mechanized  by  using  a  store  as  an  accumulator  and  binarizing 
the  result.  Comparison  operations  of  the  form: 

*1*  12  l£  I2>  11 

can  be  considered  a  special  case  of  the  threshold  operation  and  are  readily 
done  in  pip.  Another  class  of  comparison  operations  (needed  for  median 
filtering)  takes  the  form: 
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This  operation  could  be  designed  into  a  new  cell  but  is  not  readily 
implemented  in  PIP. 

Table  6  lists  a  number  of  image  processing  algorithms  and  summarizes  the  types 
of  operations  required  to  implement  each. 


A  Matrix  Processor  Example 


To  appreciate  the  actual  implementation  of  algorithms  in  PIP  we  have  chosen  a 
system  example  that  performs  representative  algorithms  from  Table  5.  The 
application  is  a  missile  seeker,  and  PIP  performs  both  target  detection  and 
computation  of  centroid  information  for  tracking. 


The  sequence  of  operations  is  shown  in  Table  7.  The  background  filter 
provides  5x5  smoothing  window  with  unity  weighting.  The  output  response 
A(i,j)  with  an  input  of  I(i,j)  iss 
2  2 

A{i, j)  -  r  #  r  j+i> 

k--2  1--2 

The  subtract  and  threshold  operation  that  follows  this  uses  the  background 
estimate  above  to  generate  a  binary  map  as  follows: 

B(i, j)  -jl  if  I  I (i, j)  -  B(i,  j)l  >  a 
|_0  if  I  I  (i,  j)  -  B(i,  j)l  <  a 


If  we  look  at  the  operations  that  follow  this,  however,  we  can  justify  a 
modified  procedure.  The  PIP  provides  a  single-sided  subtraction  that  may 
provide  equivalent  (or  even  better)  results  than  either  straight  binarization 
or  pure  gray-level  processing.  Conventional  image  processing  systems  have 
retained  two  image  maps:  one  binary  and  one  with  gray  scale.  PIP  can  easily 
provide  the  following  function: 


Ttir  «ahold/Bi  nar  i  *• 


TABLE  7.  PROCESSING  SEQUENCE/MISSILE  SEEKER 


Background  average  (5  x  S) 

Subtract  and  threshold  (5  x  5) 
Template  matching 
X-centroid 
Y-centroid 

Sum  of  absolute  differences  tracking 


i.j)  -j'llU. 


j)  -  B(i,  j)l  -a  if  I  I  (i*  j)  -  B{i,  j)l  >  a 
0  if  I  I(i,j)  -  B(i,j)l  <  a 


Note  the  following  ("target"  implies  here  only  a  "region  of  interest"): 

1.  Binarization  presumably  separates  "targets"  from  "nontargets," 
but  in  the  steps  to  follow  we  would  get  sharper  peaks  in  the 
filtering  operations  with  gray-level  information. 

2.  Gray-level  processing  does  not  separate  targets  from  nontargets 

3.  PIP's  subtraction  operation  both* isolates  targets  and  retains 
gray-level  processing  capability  (albeit  with  a  constant 
subtracted  from  the  result)  for  targets. 

Me  would  recommend  simulation  of  this  concept  to  verify  it,  but  a  binary 
thresholder  may  well  be  unnecessary  in  this  application. 


The  next  algorithmic  step  is  template  matching.  The  image  M(i,j)  generated  in 
the  previous  step  is  passed  through  a  two-dimensional  FIR  filter;  that  is. 


C  (i,j) 


■t  t2 


a  M(i+k,  j+1) 


This  is  followed  by  two  operations  that  feed  a  tracking  processor.  If  one  or 
more  points  exceed  an  external  threshold  for  a  target,  we  use  the  M  array  to 
do  centroiding  in  the  X  and  X  directions  as  follows: 


i  M(i» j) 


j 


,  2  i 

±t:l 


1—2 


j  M(i,j) 


When  implemented  in  PIP  we  perform  each  computation  in  fully  parallel  fashion 
and  search  the  external  data  stream  for  a  maximum  in  each  case.  The  final 
tracking  step  computes  minimum  absolute  difference  correlation. 


Let  lQ  (i, j)  be  the  old  frame  and  IN  (i,j)  be  the  new  frame.  We  must  find 
D ( k , 1)  such  that: 


V 


M,n  (D, 


m,  n 


) 


and  D  -EE  IIn(i,j)  -  I„(i-m,j-n)l 
m,n  i  j  0  N 

In  other  words,  we  shift  the  two  images  successively,  subtract,  sum  over  the 
array,  and  find  the  position  (k,l)  that  yields  the  minimum  value  of  this  sum 
of  differences.  We  are  now  ready  to  tackle  the  PIP  implementation  of  this 
algorithm  sequence. 
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A  variety  of  ways  exist  to  implement  a  two-dimensional  filter  in  PIP.  The 
five  basic  classes  are  illustrated  in  Figure  65  along  with  the  numbers  of 
shifts  required  for  each  method.  The  method  to  use  is  dependent  on  M,  the 
number  of  pixels  in  the  window,  but  for  all  practical  cases  methods  (b)  and 
(c)  are  superior.  Note  that  for  separable  filters  the  partial  sums  can  be 
used  in  filter  calculations  at  more  than  one  pixel;  this  is  why  the 
expressions  seem  to  give  lower  numbers  than  the  illustrations  would  indicate. 

To  compute  the  background  average  (a  separable  filter)  we  use  technique  (e) 
from  Figure  65.  The  PIP  mnemonic  sequence  would  go  as  follows: 


PIP  Mnemonic 

Function 

IA 

Load  Array 

L2 

t 

X  (image  In  x)+  REG2 

E3 

• 

9 

0  -*■  REG3 

REGENERATE  (S2,l) 

r 

Standard  Regenerate 

SUM 

• 

• 

SUM  Macro  (See  Below) 

REPEAT  4  TIMES 

; 

LOOP 

SR 

i 

Shift  Right 

SUM 

END  REPEAT 

9 

END  OF  LOOP 

U3 

Horiz.  Sum  in  REG3+ 

SUM 

§ 

Compute  Total  Sum 

REPEAT  4  TIMES 

SD 

9 

Shift  Down 

SUM 

END  REPEAT 

$ 

Sum  in  REG3 

U3 

ST3  (averaged  image  !)-► 

SU 

X 

up  one  position 

SU 

X 

up  two  positions 

SL 

X 

left  one  position 

X  left  two  positions 
X+  STl 
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SL 

LI 


Separable,  propagation  of  e.  separable,  propagation  of 

partial  sums  in  one  dimension  partial  sums  in  both  dimensions 


SUM  MACRO 


Mnemonic  Function 

START  SUM  j  X  -*■  REGl 

Ll 

DOUBLE  REGENERATE  (S1#  .2) 

L3  ;  X  +  REG3  +  REG3 

Ul  ;  REG1  +  X 

END  SUM 

We  next  perform  the  subtract  and  threshold  operation.  This  can  be  done  very 
simply  in  PIP  as  a  variation  on  the  absolute  difference  instruction  described 
earlier.  We  need  only  apply  a  gate  voltage  to  RSI  that  gives  the  desired 
offset  from  zero.  The  regenerator  function  will  then  clip  off  the  input  to 
the  Z  gate  detector  and  leave  the  result.  We  have  designated  Rl  as 
■regenerate  the  contents  of  STl  and  store  the  result  in  x."  We  can  now  define 
■Rl  -  a"  as  the  clipping  operation  described  above,  similarly,  M31 
indicates  (ST3  -  STl),  so  we  will  define: 

M31  -  a  ■  (ST3  -  STl)  -a 

(Note  that  M31  is  equivalent  to  the  existing  M2  instruction  in  function.)  The 
sequence  for  implementation  of  the  subtract  and  threshold  algorithm  is  given 
below. 


Mnemonic 

Function 

MOD 

1  I  -  Al  -*■  ST3 

El 

0-*-  STl 

GEN1  a 

a  ♦  STl 

M31 

1  i  -  Al  -  a  ♦  x 

El 

0  ■*  STl 

Ll 

(  i  -  fi  -  a  )  ♦  STl 
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As  a  starting  point  for  implementation  of  the  remaining  algorithms  we  will 
review  the  contents  of  the  cell  storage  locations,  in  all  cases  the  subscript 
indicates  frame  1. 

«  (I  l^  -  A^l  -a)  ■*>  ST1,  where  i^  »  "raw*  image,  A  -  averaged  image 
I  -*■  ST 2 

1  ll  ~  Al'  *  ST3 
0  ♦  X 

We  consider  each  of  the  first  two  algorithms  in  this  sequence  as  separate 
entities  for  implementation.  To  minimize  the  number  of  required  storage 
locations,  however,  we  should  consider  the  remaining  three  algorithms  as  a 
group. 

The  template-matching  filter  can  use  the  B  array  to  perform  filtering.  We 
will  assume  that  the  window  size  for  the  computation  of  the  template-matching 
filter  is  5  x  5  and  that  the  weighting  coefficients  are  arbitrary,  we  will 
also  assume  a  similar  size  for  the  centroid  computation.  Note  that  5x5  may 
not  be  adequate;  the  actual  window  size  must  be  larger  than  the  largest 
target  expected.  The  template-matching  operation  will  compute  256 
two-dimensional  filtering  functions  in  parallel  and  ideally  will  indicate  the 
presence  or  absence  of  a  target  that  will  then  trigger  centroiding  and 
tracking.  The  template-matching  filter  need  only  provide  an  output  when  one 
or  more  cells  indicate  a  computation  larger  than  some  threshold.  In  the 
existing  chip  such  a  determination  could  be  made  only  after  reading  out  the 
entire  array.  A  random-access  feature  is  incorporated  in  the  new  cell  design. 

If  we  unload  the  entire  array  and  find  one  template  filter  point  above  the 
target  threshold,  then  we  must  proceed  with  centroiding.  However,  since  we 
need  the  array  in  ST1  to  do  centroiding  as  well,  we  may  be  better  off 
combining  some  or  all  of  the  steps  in  the  two  operations. 


The  template-matching  operation  will  be  nearly  identical  to  the  background 
averaging  operation.  We  assume  that  the  array  to  be  used  is  *  0  - 

A^l  -  a)  stored  in  STl.  The  operation  to  be  performed  is: 

5  5 

C<i'j>  "?1  1^1  Bk'1  B<i"k'  j_1> 


As  noted  earlier,  we  cannot  assume  that  the  8^,1  are  separable  and  will 
not  do  so  here.  Table  8  presents  the  nonseparable  implementation  in  pip  of 
the  above  function.  We  will  assume  that  the  (3^,1  are  scaled  properly  to 
preclude  an  overflow  condition  with  any  of  the  C(i,j).  Note  that  the  SUM 
macro  also  includes  computation  of  the  centroid  and  stores  the  running  3ums 
for  i  and  j  in  ST4  and  STS,  respectively.  We  require  four  extra  storage 
registers  in  this  implementation,  which  minimizes  the  number  of  shift 
operations.  In  addition,  we  could  perform  the  minimum  absolute  difference 
operation  simultaneously;  we  have  included  that  in  Table  8  as  well,  we  have 
assumed  here  that  the  output  stream  of  the  template-matching  operation  is 
thresholded  off-chip  (a  trivial  operation) .  As  noted,  the  results  of  window 
operations  will  not  be  centered  but  instead  will  be  stored  in  the  lower  right 
pixel  of  the  window.  Since  all  window  results  are  to  be  read  out  anyway,  we 
do  not  shift  them  up  and  left  prior  to  readout. 

Note  that  minimum  absolute  difference  correlation  requires  data  from  a 
previous  frame.  We  shall  assume  that  the  image  used  for  this  operation  is  the 
image  B^.  We  must  assume  that  no  desired  data  is  destroyed;  one  of  the 
four  extra  storage  registers  required  must  be  used  for  this. 

Although  the  sequence  shown  in  Table  8  is  efficient  with  regard  to  minimizing 
shifting  operations,  it  is  worth  exploring  an  alternative  which  requires  fewer 
storage  registers.  Table  9  gives  the  PIP  instruction  sequence  for  the  same 
seeker  algorithms,  but  each  of  the  four  summing  computations  required 
(template  match,  X-centroid,  Y-centroid,  and  tracking)  is  done  separately.  We 
now  require  only  four  storage  registers  rather  than  seven. 
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TABLE  8.  NONSEP  ARABLE  SEEKER  ALGORITHMS,  TEMP  LATE/CENTROID/TRACKING 
IMPLEMENTATION 


PIP  Mnemonic 

Function 

1.  R1  (0.2) 

0.2  B  (i-2, j-2)  -*•  X 

2.  E4 

Clear  x-centroid  register 

3.  L4 

0.2  B(i-2, j-2)  -*•  ST4 

4.  Rl  (0.2 

0.2  B ( i-2, j-2)  X 

5.  E5 

Clear  Y-centroid  register 

6.  L5 

0.2  B (i-2,  j-2)  +  ST5 

7.  U6 

Bm-l(i'5>  ^  ^ 

8.  L2 

*  ST2 

9.  MOD 

lBm_l(i, j)  -  Bm(i, j) 1  +  ST 3 

10.  U3 

1 B  .  (i,  j)  -  B_(i*  j)  1  ♦  x 

m-JL  m 

11.  E7 

0  ST7 

12.  L7 

,Bm-l(i'j)  ‘  Bm(i'j) 1  *  ST7 

13.  U2 

*  X 

14.  L6 

Bm-l(i,;i)  ST6 

l5‘ 

16.  E3 

0  -*•  ST  3 

17.  L3 

0i-2,j-2  ^  s« 

18.  Rl  (1) 

B(i-2, j-2)  -►  X 

19.  SHIFT 

See  SHIFT  macro  below 

20.  OUTPUT  ST 3 

1  At  this  point  the  output  template  filter  results  are  externally  thresbolded  1 

I  to  determine  whether 

further  processing  steps  (centroid  determination  and 

|  tracking)  are  warranted,  if  so,  then  the  following  steps  are  performed:  ] 

21.  OUTPUT  ST4 

i  (X-centroid) 

22.  OUTPUT  STS 

j  (Y-centroid) 

23.  OUTPUT  ST6 

Tracking  results 

24.  Ul 

Bro-1  ♦  x 

25.  L6 

Bm-1  -  ST6 

L6 


TABLE  8.  NONSEP ARABLE  SEEKER  ALGORITHMS,  TEMPLATE/CENTROID/TRACKING 
IMPLEMENTATION  (continued) 


PIP 

Mnemonic 

Function 

SHIFT  Macro 

1. 

REPEAT  5  times 

j 

2. 

SR 

Shift  right 

3. 

SUM 

See  SUM  macro  below 

4. 

END  REPEAT 

5. 

SD 

Shift  down 

6. 

SUM 

7. 

REPEAT  4  times 

8. 

SL 

Shift  left 

9. 

SUM 

10. 

END  REPEAT 

11. 

SD 

• 

12. 

SUM 

13. 

REPEAT  4  times 

14. 

SR 

15. 

SUM 

16. 

END  REPEAT 

17. 

SD 

18. 

SUM 

19. 

REPEAT  4  times 

20. 

LS 

21. 

SUM 

22. 

END  REPEAT 

23. 

SD 

24. 

SUM 

25. 

REPEAT  4  times 

26. 

SR 

27. 

SUM 

28. 

END  REPEAT 
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TABLE  8.  NONSEP ARABLE  SEEKER  ALGORITHMS,  TEMPLATE/CENTROID/TRACKING 
IMPLEMENTATION  (concluded) 


In  a  conventional  digital  processor  most  instructions  require  the  same 
execution  time  in  a  given  processor.  (There  are  obvious  exceptions,  of 
course:  the  double-length  and  single-length  instructions  in  many 
microprocessors  form  two  classes.)  In  PIP,  however,  execution  of  each 
sequence  for  each  mnemonic  requires  a  different  number  of  clock  states.  This 
implies  that  not  only  the  number  of  instructions  but  also  the  types  are 
important  in  determining  total  throughput.  Table  10  summarizes  the  actual 
execution  times  required,  grouped  by  instruction  type.  These  are  expressed 
both  in  numbers  of  clock  states  and  in  percentage  of  the  total  execution 
time.  (Even  with  no  changes,  however,  PIP  can  perform  all  calculations  in  59 
msec  at  a  1  MHz  clock  rate.)  In  addition  to  possible  major  cell  changes,  we 
wanted  to  identify  missing  mnemonics.  A  summary  of  comments  regarding  the 
adequacy  of  the  existing  instruction  set  follows. 
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TABLE  9.  SEPARABLE  SEEKER  ALGORITHMS ,  TEMPLATE/ 
CENTROID/TRACKING  IMPLEMENTATION 


PIP 

Mnemonic 

Function 

Template  Match 

1. 

R1 

B.  _  .  ,  B(i-2,j-2)  X 

1-a, 3-2 

2. 

E3 

0  ♦  ST  3 

3. 

L3 

Si-2,j-2  ST3 

4. 

R1  (1) 

B(i-2,j-2)  ♦  X 

5. 

SHIFT 

SHIFT  macro  in  Table  6 

6. 

OUTPUT  ST3 

If 

threshold  results  justify  further  processing,  then  we  perform  the  1 

following  steps: 

X-Centroid 

7. 

R1  (0.2) 

0.2  B(i-2,j-2)  ♦  X 

8. 

E3 

Clear  X-centroid  register 

9. 

L3 

0.2  B ( i-2, j-2)  -»•  ST 3 

10. 

R1  (1) 

B(i-2,j-2)  -*■  X 

11. 

SHIFT+ 

12. 

OUTPUT  ST3 

Y-Centroid 

13. 

R1  (0.2) 

0.2B(i-2, j-2)  ♦  X 

14. 

E3 

Clear  Y-centroid  register 

15. 

L3 

0.2  B ( i-2, j-2)  ♦  ST3 

16. 

R1  (1) 

* 

Bm  (i-2,  j-2)  ♦  X 

11. 

SHIFT 

12. 

OUTPUT  ST 3 

Replace  SUM  macro  with  SUMl 
^Replace  SUM  macro  with  SUM2 
^Replace  SUM  macro  with  SUM3 
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TABLE  9.  SEPARABLE  SEEKER  ALGORITHMS,  TEMPLATE/ 

C ENTR0I O/TRACKING  IMPLEMENTATION  (continued) 


19. 

U4 

B»-l*  X 

20. 

E2 

0  ST2 

21. 

L2 

Vi  *  ST2 

22. 

MOO 

IVl{i'j)  -  Bm{i'j)l  *  ST3 

23. 

24. 

Rl  U) 

* 

SHIFT 

B  (i,j)  -*•  X 
m 

25. 

OUTPUT  ST4 

26. 

Ul 

Bm(i,j>  -  X 

27. 

E4 

0  -*■  ST4 

28. 

L4 

Bm(i,j)  -*•  ST4 

SUM1 

Macro 

1. 

E2 

0  -►  ST2 

2. 

L2 

B  (i-k, j-1)  >  ST2 
m 

3. 

Si-k.j-i  *  x 

4. 

L3 

Template  sum  ST 3 

5. 

U2 

Bm(i-k,  j-1)  -  X 
m 

SUM2 

Macro 

B 

E2 

0  >  ST2 

m 

L2 

+  ST2 
ro 

R2  (i/5) 

(i/5)  Bm(i-k,j-l)  ♦  X 

i  : 

L3 

X -centroid  sum  ■*>  ST 3 

5. 

U2 

Bm(i-k,j-l)  ♦  X 

Replace  SUM  aacro  with  SUM4 


TABLE  9.  SEPARABLE  SEEKER  ALGORITHMS ,  TEMPLATE/ 

CENTROID/TRACKING  IMPLEMENTATION  (concluded) 


PIP 

Mnemonic 

Function 

SUH3 

Macro 

1. 

B2 

0  +  ST2 

2. 

1*2 

Bm  ( i-k ,  j-1 )  -►  ST2 

3. 

R2  (3/5) 

( j/S)  B  (i-k, j-l)  x 

rn 

4. 

L3 

y-centroid  sum  ♦  ST 3 

5. 

U2 

Bm(i-k,  j-l)  *  x 

SUM4 

Macro 

1. 

E2 

0  ♦  St2 

2. 

L2 

^(i-k.j-i)  +  ST2 

3. 

MOD  (2,4,3) 

IST2  -  ST4I  -*•  ST 3 

°f  iBm-l(i,j)  "  Vi_k'j"l,l  +  ST3 

4. 

U2 

Bm(i-k, j-i)  ♦  x 

TABLE  10.  PIP  EXECUTION  TIMES  FOR  MISSILE  SEEKER  ALGORITHMS  (2214) 


Instruction  Type 

Clock  States 
Required 

percent 
of  Time 

I/O  (Load  Array,  unload  Array) 

48,730 

83.0 

Shift  (SL,  SU,  SD) 

777 

1.3 

Erase  (El,  E2,  E3) 

640 

1.1 

Load/unload  (Ll,  L2,  L3,  ul,  U2,  U3) 

2,565 

4.4 

Regenerate  (Rl,  R2) 

2,275 

3.9 

Modulus 

3,726 

6.3 

58,713 

100.0 
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Background  Filter — A  fourth  storage  cell  is  necessary  for  rapid 
implementation.  Multiple  passes  through  the  regenerator  are 
necessary  to  achieve  some  scaling  factors. 


o  Subtract  and  Threshold — The  modulus  instruction  can  be  modified  (as 
mentioned  earlier)  to  directly  obtain  the  desired  result. 

o  Template  Match — Random-access  output  is  desirable. 

o  centroid — No  problem,  except  that  even  sore  storage  registers  (up 

to  six)  are  needed. 

o  Minimum  Absolute  Difference  Correlation— No  additional  functions 
are  needed. 

The  final  data  required  before  we  can  discuss  architectural  improvements  in 
the  new  cell  comes  from  cell  modeling.  We  describe  that  in  the  next 
subsection. 

PIP  Modeling 

Modeling  of  the  PIP  cell  can  be  done  in  many  ways.  We  initially  generated  a 
register  transfer  model  of  PIP  (see  Figure  66)  which  treats  the  PIP  cell  as  a 
distributed  processor  architecture.  Although  useful  in  some  contexts,  this 
approach  has  a  weakness;  it  is  not  a  physical  model,  so  some  operations  shown 
occur  simultaneously  and  throughput  calculations  are  impossible.  We  chose 
then  to  perform  physical  modeling. 

The  modeling  of  PIP  is  considerably  more  complex  than  modeling  of  CCD  delay 
lines,  filters,  or  imagers  because  the  destination  of  spurious  dark-current 
charge  and  signal  charge  separated  from  the  main  packet  by  charge  transfer 
inefficiency  cannot  be  readily  ascertained.  This  occurs  for  two  reasons. 

Both  the  temporal  and  spatial  characteristics  of  such  charge  can  be  modeled  in 
a  more  conventional  CCD.  in  PIP,  however,  the  charge  left  behind  because  of 
charge-transfer  inefficiency  must  wait  a  variable  amount  of  time  for  a 
trailing  packet  to  come  along.  (This  provides  more  time,  in  general,  for 


trapped  charge  to  rejoin  the  main  packet  and,  except  in  the  case  of  the 
intercell  shift  registers,  would  probably  provide  a  hiqher  effective  transfer 
efficiency  within  the  cell.)  Also  PIP  gates  can  be  turned  on  and  off  with  no 
signal  charge  present  and  the  dark  current  generated  during  the  ON  cycle  has 
no  clear  destination.  An  example  of  this  is  the  set  of  B  and  C  gates  that  are 
not  adjacent  to  any  of  the  corner  nodes  X.  When  a  "shift-right"  operation  is 
performed  the  signal  charge  packet  moves  from  one  X  node  to  the  next.  The  B 
and  C  gates  between  the  nodes,  however,  are  on  during  part  of  the  transfer  and 
then  off  when  charge  reaches  the  X  node.  The  dark  current  generated  while 
they  are  on  must  go  somewhere  when  the  gates  are  turned  off;  it  redistributes 
itself  in  areas  of  lowest  potential.  The  exact  redistribution  could  not  be 
determined  without  solving  a  complex  two-dimensional  potential  problem  for 
every  possible  PIP  instruction  sequence — a  tedious  exercise  of  questionable 
value.  An  additional  PIP  nuance  is  the  varying  gate  size,  which  leads  to 
varying  numbers  of  dark-current  carriers  being  generated  per  unit  time  under 
the  various  electrodes.  Nevertheless,  we  have  generated  a  first-order  model 
that  can  predict  the  performance  degradation  within  a  cell  resulting  from 
dark-current  generation  and  charge  transfer  inefficiency.  It  can  be  extended 
easily  to  include  the  I/O  instructions,  shift  right,  shift  up,  etc.  To 
generate  a  first-order  model,  several  simplifying  assumptions  were  necessary 
as  follows: 

1.  A  generation  rate  of  10  n^/cm^  was  assumed  throughout  the 
active  channel  regions. 

2.  When  dark-current  charge  was  redistributed  among  various 
potential  wells,  it  was  assumed  to  stay  entirely  within  the  cell 
being  modeled.  The  charge  was  assumed  to  distribute  itself 
among  charge  packets  of  interest  (STl,  ST 2 ,  ST3,  X)  in  a  way 
proportional  to  both  their  areas  and  the  amount  of  time  each 
gate  was  on  (i.e.,  serving  as  a  potential  well). 

3.  Dark  current  was  assumed  conserved  within  the  cell,  (in 
practice,  this  means  that  no  dark  current  was  lost  through  sink 
ARV.) 
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Despite  the  varying  gate  sizes  of  the  PIP  cell,  we  assumed  an 
average  transfer  inef f iciency e  per  transfer. 

5.  The  regenerator  circuit  was  modeled  using  the  equation  derived 
in  an  earlier  quarterly  RADC  report,  namely: 

Regenerator  output  *  A  -  K2B,  where 
A  -  charge  in  Z 
B  =  charge  in  V 

2 

6.  Second-order  effects  fe  terms  and  dark-current  transfer 
inefficiency)  are  ignored. 

The  model  equations  for  each  instruction  within  the  cell  are  shown  in 
Table  11.  Note  that  x  is  the  time  for  a  single  clock  state  in 
microseconds.  Everything  is  reduced  to  numbers  of  electrons.  The  model 
divides  the  behavior  into  ideal  and  non-ideal  terms  to  aid  in  predicting 
ideal  performance  for  long  instruction  sequences.  Note  that  noise  terms 
are  omitted,  since  we  would  expect  dark  current  (particularly  for  long 
instruction  sequences  including  input/output)  to  dominate.  Nonlinearity 
in  the  regenerator  performance  can  be  readily  incorporated.  (Primed 
quantities  indicate  results  after  an  instruction;  unprimed  quantities 
are  the  original  charge  packets.)  This  model  plus  the  data  from  the 
seeker  example  give  us  the  inputs  necessary  to  discuss  our  "wish  list” 
for  design  of  a  new  cell. 

Prepared  Architectural  Improvements 

We  can  divide  the  architectural  discussion  into  three  parts  as  follows: 
o  Input/output 

o  Storage 

o  Regeneration/scaling 
These  are  discussed  below. 
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TABLE  11.  MODEL  OF  THE  PIP  INTRACELL  INSTRUCTION  SET 


Mnemonic 

Register 

Ideal  Result 

Degradation  (Ada  to  Ideal  Result) 

El 

mmmm 

0 

+0.4e  *  STl  +  7 . 4t 

ST  2 

+  1.0e  *  STl  +  18. 8t 

ST3 

+1 . Oe  »  STl  +  18. 3t 

X' 

X 

+1.5e  *  STl  +  26. 1T 

E2 

ST1  ’  = 

ST1 

+2 . 5e  *  ST2  +  58. 8t 

ST2 '  = 

0 

+0.6e  *  ST2  +  15. 2T 

ST3 '  = 

ST3 

+2 . Oe  *  ST2  +  48. 2t 

X1 

X 

+2.9e  *  ST2  +  67. 7t 

E3 

ST1 '  = 

ST1 

+  3.9e  *  ST 3  +  5 . 9e  *  X  +  180. 7T 

ST2 '  = 

ST2 

+3.1e  *  ST3  +  4.6c  *  X  +  142. 6t 

ST  3'  = 

0 

+0.6e  *  ST3  +  1.0c  *  X  +  29. 6t 

X' 

X 

+4. 3c  *  ST 3  -  11.5c  X  +  199. 5t 

R1 

ST1  '  = 

ST1 

+(-4.2e  t  1.8K  E  )  *  STl  +  69. 7t 

ST2 '  = 

ST  2 

+  (1.2e  +  1.2  ICe  )  *  STl  +  45. 6T 

ST3  ‘  = 

ST3 

+(1.2c  +  1.2  K*e  )  *  STl  +  45. 6T 

X‘ 

K1*ST1 

+(1.8c  -  7.2  )  *  STl  +  69. 1T 

R2 

ST1 '  = 

ST1 

+(3.1c  +  1.3  K  e  )  *  ST2  +  59. 4T 

ST2  ’  = 

ST2 

+(-10. 0c  +  1.71K  e  )  *  ST2  +  77. 6t 

ST  3 '  = 

ST3 

+(2.9e  +  1.2  K  E1)  •  ST2  +  56. 3t 

X' 

K1*ST2 

+ (3. Be  -  11.4  KjC  )  *  ST2  +  76. 7X 

LI 

ST1 '  = 

X 

+ (-3 . 5c  )  *  X  +  100. 9t 

ST2 '  = 

ST2 

+1.7e  *  X  +  37. 8X 

ST3 '  * 

ST3 

+1.7c  «  X  +  37. 8t 

X' 

0 

+0.2e  *  X  +  3.5t 

111 

ST1'  = 

c 

+0. lc  +  STl  +  2.5t 

ST2 '  = 

ST2 

+2 . Oe  *  STl  +  35. 0T 

ST3 '  = 

ST  3 

+2 . Oe  *  STl  +  35. 0T 

X' 

STl 

+  (-4 . le  )  *  STl  +  67. 5t 

U2 

ST1 '  » 

ST  I 

+0.7e  *  ST2  +  12. 7t 

ST2 '  = 

0 

to.le  *  ST2  +  1.8T 

ST3 '  ■ 

ST3 

+0.9e  *  ST2  +  15. 8X 

X* 

ST2 

+(-1.7c  )  +  ST2  +  39. 7t 

U3 

ST1 '  = 

STl 

+1 . 4c  *  ST 3  +  27. 3X 

ST2 '  = 

0 

+0.8c  *  ST3  +  15. 1T 

ST3 '  = 

ST  3 

-0.2e  *  ST3  +  3 . 4t 

X' 

ST2 

+ (-2 . 3c  )  *  ST3  +  39. 7X 

M2 

ST1  ’  * 

STl 

+(1.6e  +  1.6c  K  )  *  STl  +  177. 6X 
+  (7 .  Oe  -  1 . 6e  lO  •  ST2  +  2.4e  *  ST3 

ST2 '  = 

ST2 

+  (1.9e  +  1 . 9c  Kp  *  STl  +  207. 3t 
+  ( 5 . 6e  -  1.9e  lO  *  ST2  +  2.8e  *  ST 3 

ST3 '  = 

0 

+0.1c  *  ST 2  +  2.0t 

X' 

K1*STl-K2»ST2 

♦ (-16.5c  K  )  +  STl  +  6 . 5e  K  *  ST2 
+  273. 1T  2 
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Input/Output — This  function  clearly  dominates  the  current  statistics  on 
execution  times.  We  will  discuss  input  and  output  separately. 

Input  speed  can  obviously  be  increased  tremendously  by  using  fully  parallel 
array  loading.  If  a  focal  plane  array  is  mated  with  the  PIP  array  using  bump 
interconnection,  then  we  can  readily  solve  the  input  problem  if  we  include  in 
each  PIP  cell  a  node  for  inputting  signal  charge  and  add  a  gate  for 
controlling  that  input  within  the  focal  plane  unit  cell.  We  need  only  include 
in  the  PIP  cell  a  region  around  the  node  large  enough  for  bump 
interconnection.  The  ARV  region  is  a  fine  candidate  region  for  direct  sensor 
input  data. 

Outputting  from  the  array  in  serial  fashion  cannot  be  improved,  but  in  general 
we  do  not  really  need  serial  output  data  for  the  entire  array.  Instead, 
locations  of  selected  subarrays  and  the  data  therein  are  more  interesting.  A 
random-access  feature  in  PIP  would  be  ideal  for  this.  The  rationale  for 
including  a  random-access  feature  in  PIP  has  been  explored  in  a  contract  with 
Eglin  AFB.  Three  random-access  features  are  desirable  as  follows: 

o  Indicating  that  a  threshold  excession  has  occurred  somewhere  in  the 
array. 

o  Providing  the  controller  with  an  address  (addresses)  of  threshold 
excessions. 

o  Allowing  the  controller  to  randomly  address  the  array  to  obtain 
analog  data. 

Although  in  many  applications  only  the  first  level  is  required  (decision 
point  for  further  processing),  all  three  types  are  desirable.  For  example, 
the  missile  seeker  summary  in  Table  10  is  dominated  by  the  time  required  to 
repeatedly  unload  the  entire  array  when  only  a  portion  of  it  is  needed. 
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Storage — Tables  10  and  11  indicate  several  interesting  tbings  about  the 
existing  storage  architecture.  First,  since  storage  locations  each  must  be 
accessed  differently,  the  degradations  resulting  from  erasing  each  of  three 
registers  in  turn  are  not  identical.  (This  was  also  clear  in  the  second 
quarterly  report  when  leakage  values  from  the  three  stores  were  presented.) 
Next,  although  not  readily  apparent,  portions  of  the  time  required  for 
regeneration  and  modulus  functions  are  necessary  just  to  clear  a  path  for 
desired  results  by  temporarily  storing  the  value  at  location  X.  Finally, 
data  loading/unloading  probably  consumes  more  time  than  any  other  operation 
except  I/O. 

Consideration  of  all  these  factors  has  led  to  proposal  of  a  ring  store 
architecture  that  consists  of  a  racetrack  CCD  (see  Figure  67).  This 
approach  offers  numerous  advantages  over  the  present  storage  architecture. 
Some  of  these  include  the  following: 


MBENIMTOR 


Figure  67.  Ring  Store  Architecture 
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o  Mnemonics  can  be  readily  modified  to  have  any  source  and 
destination  registers.  The  only  software  change  it. 
minor — involving  control  of  a  few  extra  shifts  around  the  racetrack. 

o  The  input/output  node  of  the  racetrack  can  be  placed  between  the 

regenerator  output  and  node  X.  This  eliminates  the  time  wasted  to 
temporarily  buffer  the  contents  of  X  to  go  "around  the  corner"  to 
stores  1  and  2. 

o  The  storage  loop  can  be  operated  as  a  FIFO  or  LIFO,  which  allows, 
if  desired,  a  "stack"  architecture  in  addition  to  the  independent 
mode  used  now. 

o  Total  throughput  is  greatly  increased. 

o  Additional  storage  locations  require  no  additional  control  gates. 

Regeneration/ Scaling — The  regenerator  was  significantly  improved  on  the  3022 
chip.  Nevertheless,  we  feel  that  even  more  can  be  done  to  help  regeneration 
and  scaling.  The  existing  regenerator  has  only  one  multiplier;  it  is 
associated  with  the  noninverting  input.  Scaling  of  the  inverting  input  is 
not  possible  directly. 

If  a  second  multiplier  is  added  (associated  with  the  inverting  input) ,  then 
this  problem  can  be  corrected.  In  addition,  the  seeker  example  illustrated 
that  the  limited  dynamic  range  of  the  regenerator's  multiplying  constant 
leads  to  repeatedly  regenerating  a  result  to  obtain  proper  scaling.  This 
can  be  eliminated  by  adding  a  coarse  scaling  function  (based  on  charge 
splitting  by  area)  that  can  quickly  scale  to  a  value  close  to  that  desired. 
Only  one  regeneration  is  needed  for  any  scaling  operation. 

To  illustrate  the  throughput  improvements,  we  computed  the  number  of  clock 
states  required  for  two  of  the  PIP  mnemonic  sets.  For  M2,  the  subtraction 
instruction,  the  new  architecture  (assuming  five  stores)  requires  only  36 


129 


states  versus  65  with  the  3022.  Erase  instructions  with  the  new 
architecture  require  only  9  to  13  states;  the  old  architecture  requires  8 
to  53.  The  overall  throughput  improvement  is  conservatively  estimated  at 
40%  within  a  cell. 

When  the  entire  array  is  considered  (including  I/O)  the  improvement  is  even 
more  dramatic.  We  performed  an  analysis  of  the  missile  seeker  algorithm 
sequence  as  implemented  with  a  PIP  cell  that  includes  all  of  these 
features.  The  results  are  listed  in  Table  12.  In  Section  ill  we  will 
discuss  the  design  of  new  matrix  processor  cells. 

TASK  6 — INITIAL  DESIGN  OP  32  x  32  ARRAY  CELL 

The  layout  of  a  full  capability  cell  is  shown  in  Figures  68,  69,  and  70. 

The  elements  contained  in  the  cell  are  as  required  from  the  analysis  of  the 

TABLE  12.  EXECUTION  TIMES  FOR  SEEKER  ALGORITHMS  (New  Cell  Functions) 


Instruction  Type 

Clock  periods 

Percent  of  Total 

I/O 

90 

1.4 

Shift 

777 

12.1 

Erase 

480 

7.5 

Load/unload 

1,881 

29.4 

Regenerate 

1,820 

28.5 

Modulus 

1,350 

21.1 

6,398 

100.0 

Improvement  in  execution  time  vs  old 

cell:  91% 
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Figure  70.  Conposite  CALMA  Plot  of  Second-Generation  Matrix  Processor 


algorithms  discussed  in  Task  5  on  Software  and  System  Architecture.  Bach 
component  and  its  function  is  listed  below. 

1.  Ring  Store 

This  is  a  re-entrant  three-phase  CCD  delay  line  that  is  entirely 
within  the  cell  and  does  not  interfere  with  the  cell-cell  shift 
registers  except  for  sharing  the  C  and  x  phases  for  erase  and  FGA 
input  functions. 

There  is  communication  to  the  horisontal  cell-cell  shift  register 
through  the  gate  labeled  "XS"  that  is  similar  to  the  XS  function  on 
the  3022/2214  chip.  The  racetrack  has  five  distinct  phase  gates 
(three  plus  two  extra  for  loading/unloading  the  FGA  ports):  E,  F ,  G, 
EE,  and  GG. 

There  is  access  to  a  SINK/ INPUT  diffusion  at  the  left-middle  part  of 
Figure  69  controlled  by  Phase  C.  The  access  to  the  FGA  z  and  y  ports 
is  through  phases  EE  and  GG,  respectively,  and  can  be  seen  in  the 
center  of  Figure  69.  Entry  to  Z  and  X  is  controlled  by  phase  X  as 
before  in  the  3022/2214  chips. 

2.  FGA 

This  is  the  same  basic  design  as  used  on  the  3022  except  that  each 
input  z  and  Y  have  identically  loaded  multipliers  with  independent  MG 
controls:  MGZ,  MGY.  It  is  visible  near  the  top  of  the  cell  in 
Figure  69. 


0.2X  Coarse  Splitter/Scaler 

This  is  a  new  component  located  to  the  left  of  the  FGA,  and  it 

comprises  a  simple  three-gate  potential  equilibration  charge 

divider.  This  has  been  demonstrated  to  work  extremely  well  on  a 

Honeywell  CCD/HCT  IR  imager  chip.  Details  are  available  for  that 

chip  if  needed.  The  operation  is  very  simple.  Charge  is  loaded  into 

the  larger  gate  region  (S3)  from  the  cell-cell  shift  register  using 

gate  A.  Both  S2  and  SI  gates  are  at  ON  potentials  at  this  time,  so 

charge  from  A  equilibrates  in  channel  potential  following  A  Switch 

2 

OFF.  The  time  constant  is  diffusion-limited  (L  /D) ,  where  the 
longest  path  is  perhaps  two  shift  gates.  Phase  S2  is  then  switched 
OFF  and  divides  the  charge  into  two  regions  under  S3  and  SI,  The 
areas  are  chosen  to  give  an  accurate  0.2  ratio  of  charge  in  SI  to 
charge  in  S3.  Having  completed  that  division,  the  x  control  is 
switched  ON,  opening  access  to  the  SINK/INPUT  diffusion  as  used  for 
the  Ring  Store.  S3  can  then  empty  the  excess  charge  into  the  sink 
followed  by  X  Switch  OFF.  The  scaled  charge  can  then  be  returned  to 
the  A  phase  and  continue  in  further  processing.  The  net  effect  of 
this  is  that  virtually  any  sensitivity  of  charge-charge  can  be 
programmed  for  the  regenerate  primitive.  The  MGZ  and  MGY  controls 
can  be  set  to  any  appropriate  analog  voltage  and  then  multiple  splits 
of  0.2  each  applied  to  that  gain.  It  is  usual  in  algorithms  to 
require  a  simple  set  of  values  for  threshold  and  smoothing,  so  that 
the  controller  need  only  store  a  few  coefficients  in  the  primitive 
routine.  A  D/A  is  then  called  for  on  the  controller  VLSI  chip. 

Random  Access  Analog  Output  Transistor 

This  is  an  FET  located  at  the  mid-righthand  side  of  the  cell.  The 
gate  of  this  FET  is  a  small  extension  of  the  Y  input  gate  to  the  FGA 
and  links  the  FGA  to  the  MGY  multiplier  structure.  The  source  and 
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drain  of  this  RAAO-FET  ars  then  simultaneously  accessible  in  the 
X-direction  and  x-direction  at  the  edges  of  the  array. 

The  structure  therefore  resembles  a  conventional  DRAM  except  that  the 
data  bit  is  not  influenced  by  the  word  and  bit-line  capacitance  save 
for  a  small  gate-diode  overlap  capacitance.  A  simple  MOS  comparator 
must  then  be  supplied  with  a  common  reference  voltage.  Depending  on 
the  gain  of  the  simple  comparator — programmable  if  desired— Either  a 
digital  or  an  analog  output  can  be  obtained  in  one  or  two  clock 
states  (e.g.,  0.5  y  sec  total)  for  an  entire  array  or  any  selection 
of  lines  of  the  array. 

The  selection  of  lines  is  controlled  by  the  other  axis.  Each  of  tne 
word-lines  is  to  be  driven  either  ON  or  OFF  with  a  simple  drive. 

This  can  be  some  small  overhead  to  the  matrix  processor  or  N  pins  can 
be  added  to  the  list  of  the  total  pin-outs  for  the  chip  to  allow  the 
controller  chip  (or  Exerciser  unit)  to  drive. 

The  significance  of  this  component  is  to  be  considered  in  terms  of 
the  higher  level  algorithms  for  adaptive  sequences.  The  output  is 
not  only  very  fast  from  the  entire  array,  but  it  is  nondestructive. 
This  means  that  conditional  jumps  can  be  executed  by  the 
controller/sequencer  depending  on  the  presence  or  absence  of  events 
in  the  data  matrix. 

A  truly  smart  sensor  can  therefore  be  implemented  that  has  virtually 
no  constraints  on  programmability. 

Sink/Input  Diffusion 

This  is  an  input  node  to  the  cell  of  a  DC  sink  node,  depending  on  the 
application  of  the  array.  It  in  intended  that  a  sensor  chip  (e.g.,  a 
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silicon  imager)  can  be  indium-bumped  to  give  a  direct  parallel  input 
to  the  cell.  This  sensor  must  provide  the  mode  select  switches  for 
program  control.  A  simple  pair  of  FETs  are  required  on  that  chip  (or 
other  device;  e.g.,  tactile  input)  to  select  between  "Read  image" 
(charge)  or  "Sink  data"  to  reset  voltage. 

These  functions  are  all  provided  without  in  any  way  interrupting  the 
bit-serial  I/O  presently  available  on  the  3022/2214  chips.  The 

design  has  been  laid  out  to  CALMA  using  VLSI  design  rules  for  both  a 

2  2  2 
9x9  mil  cell  and  a  3  x  3  mil  cell.  The  9x9  mil  cell 

uses  the  same  rules  as  have  been  used  on  the  VLSI  PtSi  imager  chip 

2 

#3318  presently  under  development  at  Honeywell.  The  3x3  mil 
cell  uses  a  state-of-the-art  fabrication  process  that  is  not  yet 
readily  available.  A  simple  scaling  factor  can  be  applied  to  the 
design  to  shrink  the  cell  to  the  appropriate  cell  pitch  required  by 
the  system  application.  The  smallest  cell  achievable  on  a  given 
process  is  determined  by  the  second  metal  (including  via)  pitch. 

There  are  20  horizontal  second-metal  lines  per  cell,  as  seen  in 
Figure  68. 

TASK  7 — DESIGN  OPTIMIZATION  REVIEW  WITH  RADC 

Three  major  reviews  were  held  with  RADC;  one  was  early  in  the  program  at  SRC 
when  the  design  corrections  and  modifications  to  the  original  chip  were 
discussed  and  agreed  upon,  the  second  was  when  the  initial  32  x  32  cell  design 
was  presented  to  RADC,  and  the  final  review  for  customer  acceptance  was  at  the 
end  of  the  contract.  At  each  of  these  reviews  care  was  taken  to  consider  in 
detail  not  only  the  actual  chip  design  but  also  the  system  implications  of 
that  design.  RADC  requested  that  the  CCD  be  capable  of  operating  on  as  broad 
a  range  of  applications  as  feasible,  including  non-image  SIMD  operations. 

All  the  requirements  of  RADC  and  Honeywell  have  been  met  in  concept-.  It 
remains  only  to  verify  the  actual  chip  performance  in  a  real  application. 
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The  Exerciser  unit  delivered  to  RADC  is  sufficient  to  operate  all  tbe 
presently  conceived  generations  of  the  architecture.  Sufficient  expansion  has 
been  provided  so  that  even  extra  block  phases  for,  say.  Parallel  RAAO  can  be 
accommodated  with  only  software  development  effort. 
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SECTION  III 


PROJECTIONS  OF  FUTURE  DEVELOPMENTS  OF  MATRIX 
PROCESSOR  CCDs  WITH  PARALLEL  I/O 


In  this  section  we  present  an  estimate  of  the  follow-on  development  options 
that  seem  both  desirable  and  feasible  in  the  near  future  and  perhaps  out  to 
the  1987-88  time  frame. 

Current  CCD  and  sensor  technology,  and  near-term  advanced  technology 
development,  are  considered  for  a  feasibility  demonstration  of  a  complete 
subsystem,  possibly  by  late  1985. 

Controller  and  software/memory  optimization  will  be  described  in  the  context 
of  a  realizable  three-  or  four-chip  assembly:  sensor,  processor,  controller, 
and  optional  host  CPU. 

SENSOR/PROCESSOR 

Initial  interest  is  for  a  minimum  offset  sensor  that  provides  an  existing 
capability  of  integrated  FET  switches  for  multiplexing,  or  CCD  line-parallel 
(least  desirable)  addressing  to  mate  with  a  second-generation  matrix  processor 
array.  An  attractive  candidate  is  the  PtSi  monolithic  CCD  imager  developed  by 
RADC  and  currently  under  advanced  development  both  by  Honeywell  and  other 
manufacturers  in  the  US. 

Other  candidates  that  may  possibly  be  usable  include  inSb,  HCT,  and 
pyro-electric  (on  a  silicon  chip),  which  are  maturing  rapidly,  and  in  some 
cases  have  demonstrated  monolithic  CCD  structures  (in  the  technical 
literature)  for  self-scanned  arrays.  The  worst  nonuniformities  of  spectral 
responsivity  with  PV  operation  of,  say,  HCT  would  be  the  most  difficult  to 
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deal  with  in  a  matrix  processor,  although  we  think  not  impossible.  An 
exception  may  be  automatic  compensation  methods  from  off-chip  references  and 
shutterless  compensation  algorithms  techniques  executable  by  the  processor, 
as  mentioned  in  Task  6. 

Thermal  mismatch  for  cooled  sensor  hybrids  must  be  considered  in  detail  when 
using  non-silicon  sensor  arrays.  Considerable  experience  has  already  been 
gained  by  Honeywell  (SRC/EOO)  with  indium  bump-connected  HOT  and  LIT  PV 
arrays.  This  technology  has  been  under  development  using  IRSD  funding  and 
also  on  several  contracts  dating  from  about  1976,  and  is  now  at  a 
sufficiently  advanced  stage  for  use  on  matrix  processor  development. 

A  special  PtSi  imager  chip  would  be  designed  to  allow  either  line-parallel 

readout  to  line-parallel  input  circuitry  on  the  new  matrix  processor  or, 

more  ambitiously,  a  full  two-dimensional  parallel  accessed  array.  Even  with 

2 

limited  design  rules,  an  8  x  8  mil  processor  cell  seems  practical  with 
current  technology.  This  implies  a  32  x  32  array  on  a  320  x  300  mil  chip 
with  a  full  set  of  peripheral  circuitry  for  RAAO  readout.  The  operation  of 
such  an  array  would  initially  be  using  the  Exerciser  unit  already  developed 
since  this  gives  greatest  flexibility  and  serves  also  as  an  algorithm 
development  system  for  parallel  operators.  However,  early  design  of  at 
least  part  of  the  final  chip  for  control  of  the  matrix  processor  will  be 
undertaken  on  Honeywell's  IfUD  funding  to  finalize  the  processor-controller 
,  partitioning  of  functions. 

SOFTWARE  AND  CONTROL 


Maximum  benefit  from  the  matrix  processor  concept  will  only  accrue  when  we 
miniaturize  the  controller  as  well  as  the  matrix  processor/sensor  array. 
Although  an  Exerciser  unit  serves  valuable  functions  in  both  validating  the 
matrix  processor  concept  and  testing  various  combinations  of  primitives,  the 
unit  is  too  large  (and,  frankly,  unnecessarily  powerful)  for  a  system 
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application.  We  £eel  that  design  o£  a  controller  chip  is  an  important  step 
toward  development  of  a  very  compact  sensor/processor/controller  module 
which  could  revolutionize  the  approach  to  applications  requiring  moderate 
resolution  and  high  throughput. 

The  control  architecture  being  used  in  the  Exerciser  box  is  conventional — a 
2910  microprogram  sequencer  steps  through  successive  words  of  ROM.  No  new 
primitive  may  begin  until  its  predecessor  is  complete.  We  foresee 
improvements  in  two  areas  of  the  controller  that  would  facilitate  its 
integration  on  a  single  chip:  architecture  and  storage. 

Honeywell  has  developed  (as  part  of  the  DoD  VHSIC  program)  an  innovative 
approach  to  microprogram  sequencer  design  that  incorporates  two  stacks — one 
address  stack  and  one  counter  stack.  The  counter  stack  concept  allows 
considerable  flexibility  in  vesting  of  unfinished  loops,  such  a  feature 
would  be  a  great  asset  in  upgrading  the  controller  to  respond  to  commands 
like  "SR5, "  which  would  require  repetition  of  one  primitive  and  also 
repetition  of  a  subprimitive  within  it.  In  other  words,  "SR*  includes  a 
partial  sequence  that  controls  the  A,  B,  and  C  clocks;  that  sequence  is 
repeated  three  times  within  the  primitive.  We  need  only  store  that  sequence 
once  and  loop  through  it  three  times.  We  can  also  use  the  same  sequence  for 
"SL"  but  decrement  addresses  rather  than  increment. 

Another  architectural  improvement  would  use  the  greater  flexibility  of  the 
new  cell  design  developed  in  this  program.  In  Section  II,  Task  5,  we 
outlined  the  possibility  of  overlapping  instructions  to  enhance  throughput. 
This  may  now  be  possible  because  operations  in  different  parts  of  the  cell 
can  be  carried  out  simultaneously  with  little  or  no  interference.  (In  the 
2214/3022  design,  each  clock  often  performed  functions  at  several  locations 
within  the  cell) . 
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Storage  of  primitives  in  the  present  Exerciser  unit  is  in  48-bit  words, 
where  up  to  32  bits  may  be  allocated  for  matrix  processor  clock  states. 

This  is  desirable  from  an  experimentalist's  point  of  view,  since  changes  in 
primitives  can  be  readily  inserted.  In  a  system  application,  however,  we 
need  not  waste  so  much  storage  when  it  is  highly  likely  that  only  a  few  bits 
change  each  time.  By  either  coding  only  differences  in  states  from  cycle  to 
cycle  (similar  to  image  bandwidth  compression)  or  using  ROM  compression 
techniques,  we  could  substantially  reduce  the  memory  size  required  to 
implement  a  set  of  primitives  in  a  system  application.  Such  a  reduction 
obviously  makes  a  single-chip  implementation  simpler. 

Honeywell  is  using  1982  internal  funding  to  investigate  controller 
architectural  and  storage  improvements.  This  work  will  culminate  in  a 
design  and  (if  contract  funding  can  be  identified)  fabrication  of  a 
single-chip  controller  that  would  take  a  major  step  toward  system 
applications  for  the  matrix  processor  concept. 
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SECTION  IV 


SUMMARY  AND  CONCLUSIONS 


An  architecturally  novel  Parallel  CCD  Matrix  Processor  has  been  developed. 

The  processor  is  fully  programmable  and  in  its  most  elaborate  form  can 

exceed  the  throughput  of  virtually  any  available  image  processor  presently 

known.  At  the  heart  of  the  device  is  a  unique*  linear,  programmable 

floating  gate  amplifier  (PGA),  which  is  small  enough  to  fit  into  a  3  x  3 
2 

mil  pixel  together  with  six  storage  sites,  a  fixed-ratio  scaler,  and 
SINK/INPUT  node.  One  axis  of  the  two-input  FGA  is  used  additionally  as  a 
control  for  periphery-access  in  a  random  access  analog  output  readout  mode. 
This  cell  with  a  full  complement  of  features  represents  the  state-of-the-art 
of  packing  density  for  a  general-purpose  signal  and  image  processing 
algorithm  implementing  IC. 

Although  not  completely  debugged,  the  two  chips  fabricated  and  tested  have 
shown  all  the  basic  design  concepts  and  parameters  to  function  as  intended. 
It  remains  to  exercise  a  "perfect"  chip  and  demonstrate  a  single-frame, 
high-order  image  processing  algorithm  for  contrast  enhancement  using 
threshold  and  edge-extraction  subroutines.  The  first  demonstration  of 
partial  operations  has  been  achieved,  albeit  with  assistance  from  the 
host/controller  CPU.  There  are  no  known  design  concept  barriers  to 
implementing  the  first  fully  programmable  smart  sensor  system.  Fabrication 
limits  to  uniformity  of  offset  in  the  FGA  have  as  yet  prevented  a  complete 
demonstration,  but  this  is  expected  to  be  accomplished  in  the  near  future. 

A  flexible  Exerciser  and  Algorithm  Development  System  (microprocessor-based) 
has  been  constructed  and  exercised  using  a  multiplicity  of  signal/image 
sources.  Both  real-time  and  freeze-frame  modes  of  activity  are  accommodated 
by  this  powerful  unit,  which  provides  all  the  necessary  drive  voltages  to 
operate  both  first-  and  second-generation  CCD  matrix  processor  chips. 
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A  cheap,  robust,  multispectral  subsystem  of  a  few  watts  power  requirement  is 
foreseen  as  a  fabrication  goal.  Leading  up  to  that  goal,  construction  of  a 
second-generation  CCD  chip  and  an  appropriate  image  sensor  is  seen  as  a 
necessary  stage  of  development,  included  in  that  development  is  the  design 
and  possible  fabrication  of  a  programmable  microcontroller/clock-driver  VLSI 
CMOS  chip. 

Unique  analog  SiMD  algorithms  are  now  provided  with  hardware  to  speed 
software  development.  Many  unique  arithmetic  features  of  the  CCD  processor 
provide  new  opportunities  for  efficiency  and  speed  of  execution  of 
ultra-high  throughput  parallel  operators.  There  remains  a  great  deal  to  be 
researched  and  tested  with  this  new  processor  but  construction  is 
demonstrated  to  be  feasible  on  a  standard  IC  CCD  process  line  to  yield 
full-capability  arrays  of  perhaps  32  x  32.  An  advanced  CCD  fabrication  line 
would  permit  up  to  128  x  126  arrays  with  only  geometry  scaling  required  from 
the  larger  cell  design. 
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APPENDIX  A 


HP41C  POCKET  COMPUTER  LISTING 


The  following  pages  give  the  HP41C  pocket  computer  listing  for  general 
single-gate  electrode  calculations  of  CCD  and  Floating  Gate  CCD 
one-dimensional  approximations. 
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4i«  XE?  95 
419  ♦ 

429  SCI  i? 

4?i  * 

422  PCI  49 
42?  / 

VE-?  47 
425  PTH 


avfg 

SECOND  PASS? 
STOP 


USE  Xw 
M 

SECOND  PASS 
COMPUTE 

COMPUTE  X'  OR  X 
O  O 


429*1  PL  «5  COMPUTE  qN 
42?  PCL  1 ! 

424  pa  i: 

4?9  * 

410  PTM 


4? '.♦•.PL 

A* 

4?2  PCL 

19 

4?? 

r  *-  ■ 

90 

4?4 

PCI 

?0 

4?5 

F?r 

40 

436 

♦ 

43' 

ytj 

43? 

SCI 

30 

439 

* 

440 

- 

44! 

PCL 

0? 

442 

♦ 

44  7 

PCL 

1? 

4*4 

■E9 

41 

44? 

07 

446  XE* 

00 

447 

xa 

44? 

♦ 

44r- 

PTH 

COMPUTE 

*2  *  2cp  esi  N  *  VPB  -  B(X>  -  V2 

2  “o'si  -  q% 


OR  X0 


45P*L0L  0?  COMPUTE 

451  PCI  12  T  (2Beoeci  ~  <aWr»> 

452  PCL  30 
45?  *£$  8! 

454  VE0  05 

455  - 

456  / 

45?  PTH 


45  -nPL  0? 


455 

F5? 

00 

464 

PCL 

21 

46! 

pr* 

80 

462 

PCL 

22 

463 

PTH 

464*L8L 

85 

465 

•XDEP71 

46t 

w;ek 

467 

95E 

4hi- 

0 

STO 

?4 

X: 

■u  -pvR 

471 

PTH 

4?2*L5l 

10 

473 

VE0 

05 

474 

PCL 

05 

475 

1  E- 

-12 

476 

• 

47? 

i 

478 

PCL 

04 

474 

» 

400 

$Tfl 

26 

441 

Ra 

23 

♦82 

» 

48? 

CHS 

484 

PCL 

09 

485 

♦ 

♦  PC 

ST? 

*7 

C» 

♦87 

PCL 

12 

488 

7 

♦85 

• 

450 

XE9 

«7 

491 

STft 

25 

492 

PTH 

♦4?*L0l 

1! 

454 

CE  88 

495 

E«l 

RETURN  X„  OR  X„ 
O  O 


ON  ERROR 
DISPLAY 
XDEP  <  0 


RESET  X  = 
RESTART1 


COMPUTE  E,  D,  G 


t  (2Beoesi  •  qV 


BREAK  LOOP 
AND  STOP 
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ANALYSIS  OF  BCCD  FLOATING  GATE  SENSING 


Refer  to  Figure  Al.  The  analysis  starts  by  assuming  charge  balance  in 


conditions  (a) 

and  (b) : 

So,  „0 

„0 

0  + 
VFG 

°D1  “  0 

(A-l) 

<4+ 

“  0 

(A-2) 

where  =»  surface  space-charge  in  x^»  X^. 

Conservation  of  charge  yields: 

Q  Q  0  0 

C  V  +0  ■  c  V  +0 

S  FG  yFG  S  FG  UFG 

From  Haken  (Reference  4),  we  know  that 

4  •  -nVro 


(A-3) 


(A-4) 


QFG  *  ~qNDXlAFG 


(A- 5) 


where 


Xl  *  X, 


»•[ 


2  *Si  ‘»i  - 
0  ' 


♦  V>  ' 


1/2 


(A-6) 


and  X^  is  of  the  same  form  as  (A-6)  but  X*  replaces  XQ  and  Xj  -  X(1  replaces  X^ 


.  CS1(2AXJ  C0X  *  q  D* 

"  C0X12*Si  -  ’"d> 

€si  (2A(Xj  -  Xm)  Cox  +  gN0) 
C0X  (2*Si  “  qND> 
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'  ■.?'  , 


Combining  the  above  equations  gives  a  simple  expression  for 


2  v°  qNP  X1AFG  +  qNDXlAFG 
VFG  *  VFG  "  Cg  Cg 

Q  0 

By  defining  the  following  quantities,  a  solution  for  VpQ  -  Vp(.  is  found; 


(Notice  x!-  is  a  function  of  V?_.) 

1  rvs 


D  =*  V 


,0 

FG  '  C„ 


qNDAFG 


.2  2Esi  (a  (X j  -  Xm)  2  ♦  V^) 


F  *  X°  "  ^SiA  -  <JND 


*SiA  -  qN0 


Then  (A-8)  becomes:  v?G  ■  D  ♦  EXo  -  E  [p  ♦  ovro] 


which  solves  for  as: 

ru 


AV  -  -  EXX  +  EXq  +  (E2G/2)  -  \/dE2G  +  (E3GXq  )  +  (eV/4)  +  E2F  (A>10) 


This  is  the  desired  solution  for  constant  load  capacitance  Cg. 


For  the  #3022  chip  the  total  load  capacitance  with  vre  ■  10v  (substrate  » 
-3V)  is  0.154  pF  and  the  floating  gate  active  area  is  1093.3 ym  . 

With  these  parameters  in  the  BASIC  program  "FGBCCD, *  which  computes  A v 

FG 

and  percent  change  in  slope  (not  given  here) ,  the  responses  shown  in  Figures 
Bl  and  B2  (Appendix  B)  were  generated. 


APPENDIX  B 


HP9845C  DESKTOP  COMPUTER  BASIC  PROGRAM 


The  following  pages  give  the  HP9845C  desktop  computer  BASIC  program  for 
several  cases  of  one-dimensional  CCD  and  floating-gate  CCD  approximations. 

The  programs  in  most  instances  produce  hard-copy  graphics  over  user-specified 

ranges  of  input  variables. 


10  '  THE  NAME  OF  THIS  PROGRAM  IS  "SPFCCH".  IT  COMPUTES  THE  CHANCE  IN  THE 

FLOATING  GATE  VOLTAGE,  Dtlvfg,  FOR  A  SURFACE  N-CHANNEL  DEVICE.  THE 
20  '  VARIAL8E  PARAMETERS  ARE  INPUT  FROM  THE  KEYBOARD. 

30  1 . 602E- 1 9  1  C-  *- 

40  PRINT  -  TYPE  IN  THE  VALUES  OF  THE  PARAMETERS  Np,  Tox,  Lfg,  Hfg,  Ci,  Q*tg 
ofg,  Vf  6 “ 

30  PRINT  "  “ 

60  PRINT  “  N*  in  unit*  of  IE15  cm — 3“ 

70  PRINT  •  " 

80  PRINT  “  To x  in  ANGSTROMS  <»  •  MICRONS  *  10,000>" 

90  PRINT  “  " 

100  PRINT  "  Lfg  in  MICRONS  <»  *  ANGSTOMS  *  IE-4  «  #  CM  ♦  1E4'" 

110  PRINT  *  “ 

120  PRINT  “  Wfg  in  MICRONS" 

130  PRINT  “  “ 

140  PRINT  "  C*  in  picoFARABS" 

130  PRINT  ■  " 

160  PRINT  “  Dug  in  pi  coCOULONBS” 

170  PRINT  “  " 

180  PRINT  “  Vofg  in  VOLTS" 

190  PRINT  "  " 

200  PRINT  “  Vfb  m  VOLTS" 

210  Input  N»,To..Lfg.Wfg,C*,Q*ig,Vofg,vfb 

220  N«»N**1E15  1  cni~-3 

230  To**To *  *  1 E-3  i  cm 

240  Lfg«Lfg«lE-4  '  cm 

250  Wfg«Wfg*lE-4  '  cm 

260  C*«C**1E-12  '  F 

270  Afg»Lfg*Nfg*lE8  !  un-2 

280  Q» i  gu*2"Q» i  gxAf  g  !  pC'uP-'B 

290  Qs i g*Q» i g»lE- 1 2  >  C 

300  E1*etron»«-Q»ig< Q  1  *~ 

310  El*ctu«2»-Qsig  i  Afg*0)  '  •-'um-'i 

320  EO-8.86E-14  '  Fxcui 

330  E» i « 1 1 . 7 

340  Eox»3. 9 

330  Cox»Eo*Eo*'Tox  i  F'cb'-j 

360  Afg-Lfg*Ufg  !  c«i~2 

370  N»C»  x  <  Af  g*Cox  ^ 

380  Vo«Q»N**E* i *Eo  -Co x~2 

390  8*t  «*2*Af  g^  C»*<  G*i  g  Afg+Cox»-  Vo  *2*2* Vo*'  Vof g-Vfb>  >'  .  3+Af  g*Cox'-2*Vo  C*> 

400  Gl«'Q»ig  Afg+Cox*SQR<  Vo-2*2*Vo*<Vofg-Vfb.<  1  :  '2 
410  G2»Cox''2»<Vo'2»2«Vo*<  Vofg-Vfb*0»ig'Af  g-  Cox)  > 

420  G»«m>»»<  Af  gxC»'-'2«<Gl-G2> 

430  Di  »c  «SQR  v  B«t  «^2-4*G*in«i») 

440  D*  1  wf  g«< -B*i  »*Di  »c  >-'2 

430  Tox»Tox«lE8  !  Angstrom* 

460  Afg«Afg*lE8  !  M  •■'2 

470  C»«C»*1£!2  '  pF 

480  Q*ig>Qstg*lE12  1  pC 

490  Cox«Cox* 1E4  '  pF'um~2 

300  PRINTER  IS  7,1 

310  PRINT  "Progr««i  SRFCCH" 

320  FLOAT  l 

330  PRINT  “N«  <C«r'-3*"|  SPA<  t  7>  f  N* 

340  FIXED  0 

330  PRINT  "Tox  i Angst ro»» > " i SPA< 1 2> | Tox 
360  FIXED  1 

370  PRINT  "Afg  < M i c non* '2 > " | SPA< 1 2  )  I Af g 
380  FIXED  5 

390  PRINT  "C»  <  P i  c of  «r  »<js  > "  j  SPA<  1 2  > | Cl 

600  FIXED  2 

618  PRINT  "Qstg  < P i coeou 1 o»b» > " r SPA < 8 > j Os t g 
620  FLOAT  3 

630  PRINT  "0»ig  <  p  i  c  oCoul  ombi-'um  '2 ) -  )  SPA<  3  •  I  Q»  t  gui*2 
640  FIXED  1 
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630  PRINT  “Vofg  < Vol  t % > “ j SPA< 13 > | Vof g, ” Vf b  t Vol t s > “ i SPA<2; i vfb 

660  FLOAT  4 

670  PRINT  “Cox  <P i  cof»r»di<'Nicron''2 '  “  i  SPR(  2 '  j  Cox 

600  FIXED  2 

690  PRINT  “N  (loading  fac t or > “ | SPA< 9> | N 
700  FLOAT  3 

710  PRINT  “Electron*” | SPA< 18>|El«ctron» 

720  FIXED  1 

730  PRINT  “E 1  «c  t  ron*/'Mi  cnon'2”  |  SPA<9>  |  El  «c  t  u»2 
740  FIXED  3 

730  PRINT  “lot  4“ | SPA <  23 ) | B« t  a 
.  760  PRINT  “Ga*«a“ t  SPA( 22  > | Gaaaa 

770  PRINT  “Dalvfg  ( Vo  1 1 » > “ J SPA< 1 3 > | D* 1 vf g 
780  PRINTER  IS  16 
790  END 


10  PRINT  "Th*  name  of  this  program  is  'GLSFCL'  (GeneraL  SurFac*  ChanneL).  It 
computes  Delvfg,  the  chtngt  in  th*  * 

20  PRINT  "floating  gat*  voltage  Vfg  for  *  surface  N-chann*l  device.  It  alio  c 
•  I  cul  ates" 

30  PRINT  "th*  partial  derivative  of  Delvfg  with  respect  to  Qsig.  Graphical  ou 
t  put  i s “ 

*0  PRINT  "optional.  An  INPUT  asking  for  '  ( YrN> '  Means  type  Y  for  a  Yes  respon 
s*  and  N  for" 

SO  PRINT  "a  No  response.  Any  response  other  than  Y  or  N  »ay  result  in  erroneo 
us  operation" 

60  PRINT  “of  th*  progra*. “ 

70  PRINT  "  " 

00  OPTION  BOSE  1 

90  DIM  Qs i g< 51 ' , Qs i gpc < 31 >, D* I vfg(31 >, Pd*rvq<31 >, Pvq(51 > 

100  i 
110  I 

120  PRINT  “Th*  first  section  of  th*  program  calculates  th*  Chang*  in  Vfg  as  a  f 
unc t i on  of “ 

130  PRINT  "Qsig  for  N  values  of  Qsig  defined  on  an  interval  Qsig<l)  to  Qsig<N) 
via  a" 

1*0  PRINT  "FOR-NEXT  loop.  0  s i ng 1  *  point  <N«1)  computation  is  allowed.  Press 
CONTINUE  to" 

150  PRINT  "access  th«  first  section  of  the  program." 

160  PAUSE 

170  PRINT  L I N v  4  > 

180  Q« 1 . 602E- 1 9  '  C 

190  PRINT  "TYPE  IN  THE  VALUES  OF  THE  PARAMETERS  Hi,  Tox,  Lfg,  Wfg,  Cs.  Vofg,  Vf 

b" 

200  PRINT  “  " 

210  PRINT  "  Na  in  units  of  1E13  c»"-3" 

220  PRINT  “  " 

230  PRINT  "  To<  in  ANGSTROMS  C«  4  MICRONS  •  10,00O>" 

2*0  PRINT  "  “ 

230  PRINT  “  Lfg  in  MICRONS  c«  *  AnGSTOMS  *  IE-4  «  •  CM  •  1E*>" 

260  PRINT  “  “ 

270  PRINT  "  Wfg  in  MICRONS” 

280  PRINT  “  " 

290  PRINT  *  Cs  in  picoFARADS" 

300  PRINT  “  " 

310  PRINT  "  Vofg  in  VOLTS" 

320  PRINT  “  " 

330  PRINT  "  Vfb  in  VOLTS" 

3*0  INPUT  Na. Tax, Lfg, Mfg.Cs, Vofg, Vfb 

350  Na»Na»lE15  1  cm~-3 

360  Afg»Lf  g*Mf  g  '  um"2 

370  PRINTER  IS  7, 1 

300  PRINT  "Program!  GLSFCL" 

390  FLOAT  1 

400  PRINT  “Na  <Cm"-3)"lSPAcl21fNa 
410  FIXED  0 

*20  PRINT  "To»  < Angstroms )"|SPA<7>|Tox 
430  FIXED  1 

440  PRINT  "Afg  < M i c rons "2 ) " | SPA< 7 ) | Af g 
430  FIXED  5 

460  PRINT  "Cs  (Picofarads)"! SPA( 7>| C* 

470  FIXED  l 

480  PRINT  "Vofg  c Vo  1 1 s > " ! SPAi 10 > | Vof g, * Vf b  ( Vol l * > " | SPA< 2 ) | Vf b 

490  ToxmTox#lE-8  '  cm 

300  LfgmLfgelE-4  1  cm 

310  Hfg«Mfg*lE-4  1  cm 

320  C*>Cs*lE-12  1  F 

330  Eo-O.OOE-14  •  Frcm 

340  Es I *1 1 . 7 

330  Eox-3.9 

360  Cox»Eo*EoxxTox  '  Frcm "2 

370  Afg»Lf g»Hf  g  1  c**2 
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980  Lf-Csr<Afg*Cox> 

990  FLOAT  4 

600  PRINT  “Cox  <F-'ca''2>*|SPA<10>|Cox 
610  FIXED  2 

620  PRINT  “LF  Ooading  fac  t  or  >  "  |  SPA  ( 3 )  i  Lf 
630  PRINT  LINC2) 

640  PRINTER  IS  16 

690  Vo«Q*Na*Esi  »Eo/Cox-'2  •  V 

660  INPUT  “N,Qsig<l>  pC  ,Qsig<N>  pC  ” , N, Qs i g 1 o, Qs i gh i 
670  REDIN  Qs i g<N) , Qs i gpc (N) , De 1 vf g<N> , Pdervqt N) , Puq(N) 

680  Qsiglo*Qstg1o*lE-12  •  C 

690  Qs i ghi »Qsi ghi *1F -12  1  C 

700  PRINT  “Qs i g  <pC ) “ , “ De I  of g  <V>“ 

710  PRINT  “  “ 

720  FOR  I»1  TO  N  STEP  1 

730  IF  N-l  THEN  Qs i g< I) >Qs i g 1 o  <  C 

740  IF  N-l  THEN  760 

790  Qsi g< I )«Qsi g) o+(Qsi ghi -Qsi g 1 o >*< 1-1 >✓ <.N-1 >  !  C 

760  Bet a-2*Af g'Cs»<Qsi g<  I  )  'Rf  g*Cox*<  Vo"2*2»Vo» <  Vofg-Vt'b  >  > -  .  S*Afg*Cox  '2* Vo  Cs  ■ 1  V 
770  Gl-<Qsig<I>rAfgtCox*SQR<Vo'2+2*Vo*<Vofg-Vfb>>>~2  1  C'2  ca>4 

780  G2*Cox-'2*<  Vo'2+2*Vo*CVof  g-Vfb*Qs  i  g(  I  >/Af g/Cox '  >  1  C  2'C»  4 


790  Gaaea-<Afg'Cs)~2*<Gl-G2>  <  V'2 

800  Disc-SQR<Beta'2-4*Gaaina>  •  V 

810  D« t vf g< I ) -<-Bel a»Di sc >'2  1  V 

820  FIXED  3 

830  Qsi gpc < I )»Qsi gil )* 1E12  1  pC 


840  PRINT  Qsigpc<I>,Delufg<I) 

890  NEXT  I 

860  PRINT  “End  of  Output'1 

870  ! 

680  j 

890  PRINT  LIN(2) 

900  PRINT  “Press  CONTINUE  to  access  th»  n»xt  section  of  the  progru  which  cslcu 

Utes  the" 


910  PRINT  “partial  derivative  of  Delvfg  with  respect  to  Qsig  'the  slope  of  the 
Delufg  us“ 

920  PRINT  “Qsig  graph)  using  the  parameters  end  Qsig  arrev  from  the  previous  se 
ct i on.  " 

930  PAUSE 

940  PRINT  LIN(9) 

990  PRINT  "Qsig  pC","Pderuq  (in  units  of  1E+12  V/C>” 

960  FOR  1-1  TO  N  STEP  1 

970  Bet  «-2*flf  grCs*<.Qs  ig<I)''Afg*Cox*'Vo''2  +  2*Vo*<Vofg-Vfb.')'.9*Afg*Co*'2*VoC*'|V 
980  Gl*<Qsig( I ) 'Afg*Cox»SQR< Vo~2+2*Vo*(Vofg-Vfb>  >  >~2  !  C*2  c«  4 

990  G2«Cox~2*<  Vo~2+2*Vo*<  Vof  g-Vf  b+Qs  i  g(  I  >'Af  gxCox  >  )  ■  C'2'ca  4 

1000  Geaaa-<AfgrCs)'2*<Gl-G2>  ■  V--2 

1010  X- 1  rSQRC  Bet  a^2-4*Gaaaa)  '  V-'-l 

1020  Y«<  Af  g/'Cs>“‘2*<2»<Qsig<  I  )rflfg*Cox*<:Vo''2*2*Vo*<Vofg-Vfb  >  >  - .  9  >-2*Vo*Cox)  f  V~2  C 
1030  Pdervq<I)-l-'Afg*<.9*<-ltSeta*X>*<2»Afg-'Cs>-X*Y>  !  V  C 

1040  Puq< I >-Pdervq< I )/lE12  1  VrC 

1090  FIXED  3 

1060  PRINT  Qs i gpc  < I > | SPR< 10) | Pvq( I > 

1070  NEXT  I 

1080  PRINT  “End  of  Output” 

1090  PRINT  LINUO) 

1100  INPUT  “Do  you  desire  a  graphical  output  of  this  data'  (Y  N)”,Ansl# 

1110  IF  flnsl*-“N“  THEN  1760 

1120  INPUT  “Do  you  want  the  graphical  output  on  the  9872R  or  9345C '  >72A  43C',”,A 

ns2* 


1130  IF  Ans2*-“49C“  THEN  1790 

1140  >  9872R  PLOTTER  ROUTINE - 

1190  PRINT  “BE  SURE  THERE  IS  PAPER  ON  THE  PLOTTER  then  press  CONTINUE" 
1160  PAUSE 
1170  PRINT  PAGE 

1180  INPUT  “Do  you  need  to  draw  new  axes'  ( Y-N ' " . Ans3* 

1190  IF  Ans3»-”N“  THEN  1980 
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Vy  "  > 


1200  PLOTTER  IS  7,9,"98?2A“ 

1210  LIMIT  20,260,20,190  •  am 

1220  LOCATE  10,120,10,97  !  GDUs 

1230  INPUT  “ENTER  THE  MAX.  VALUE  TO  BE  PLOTTED  ON  THE  X  AXIS.Y  AXI S“ , Xmax, Yaax 

1240  SCALE  0, -Xmax, 0, Yaax 

1290  PEN  1  '  b  1  ack 

1260  LINE  TYPE  1 

1270  AXES  1-10,110, 0,0, 10, 10, 2 

1280  I  -  LABEL  ROUTINE 

1290  '  -  X  »xis  title  1 

1300  SEC  I 

1310  MOVE  .9, -.2  • 

1320  CSI2E  9,  .5  i 

1330  LORG  6  ■ 

1340  LSIR  0  i 

1390  LABEL  USING  "K'VQsig  <pC)“  ! 

1360  1  -----  x  axis  numbers  » 

1370  CSIZE  3,  .9  1 

1380  LORG  9  • 

1390  FOR  A»0  TO  -Xmax  STEP  .1  • 

1400  MOVE  A,-. 1  ! 

1410  LABEL  USING  “ K “ 1 “ - “ J  A  ! 

1420  NEXT  A  t 

1 430  !  -  Y  ixit  title  1 

1440  MOVE  -.08,1.9  ' 

1490  CS12E  4, .9  i 

1460  LORG  6  i 

1470  LSIR  90  I 

1480  LABEL  USING  “K"|"Oelufg  <V>“  ! 

1490  1  - Y  i>i>  numbers  1 

1900  CSIZE  3, .9  • 

1910  LORG  8  i 

1920  LOIR  0  i 

1930  FOR  A-0  TO  Yaax  STEP  1  • 

1940  MOVE  -.01, A  i 

1990  LABEL  USING  “K"|A  i 

I960  NEXT  A  ! 

1970  '  - - - 

1980  PEN  4  i  red 

1990  Qs i gpc  1 1 >«-Qsigpc<li  1  pC 

1600  PLOT  Qs igpc < 1 > , De 1 uf gt 1 > , -2 

1610  FOR  J-2  TO  N  STEP  1 

1628  Qs  i  gpc  <  J>»-Qt  igpc  <  J  ■>  <  pC 

1630  PLOT  Qs i gpc  tJ),Delvfg(Ji,-l 

1640  NEXT  J 

1690  1  -  t  ox  label 

1660  Tox«Tox»lE8  *  Angstroms 

1670  PEN  1  i  black 

1600  CSI2E  2,. 9 
1690  LORG  2 

1700  MOVE  -Xmax*. 0089, De 1  of g<N> 

1710  LDIR  0 

1720  LABEL  USING  “K " | " t ox«" | Tox | " A" 

1730  '  - 

1740  PENUP 
1790  PEN  0 

1760  INPUT  "Do  vou  wish  to  run  the  program  again  with  a  new  set  of  parameters'" 

Y 'N  >  " , Ant 4S 

1770  IF  Ans4S»"Y“  THEN  180 
1780  STOP 

1790  i  9849C  CRT  GRAPHICS  ROUTINE - - - 

1000  PRINT  PAGE 

1810  INPUT  "Do  you  need  to  draw  new  axel7  < Y/N) " , AnsS* 

1820  IF  Ans9*»“H"  THEN  1940 
1830  PLOTTER  IS  “GRAPHICS" 

1840  •  GRAPHICS 


1890  LIMIT  8,148,0,140  >  a* 

1868  PEN  2  >  rad 

1870  LINE  TYPE  1 
1880  FRAME 

1890  LOCATE  10,97,10,97  >  GDUs 

1900  INPUT  “ENTER  THE  MAX.  VALUE  TO  BE  PLOTTED  ON  THE  X  AXIS.Y  AXIS" , X«a Y*a* 
1910  SCALE  0, -Xmax, 0, Ymax 
1920  PEN  1  •  whit* 

1930  AXES  1/-10,  l'lO.O.O,  18,  10,2 
1940  Qai gpc < 1 gpc < 1 )  !  pC 

1990  PLOT  Qai gpc  <l),D«1ufg(l),-2 
1960  FOR  J«2  TO  N  STEP  1 
1970  0»igpc<T>a-0slgpc<J)  1  pC 

1980  PLOT  Qai gpc < J ) , D« 1 of g< J > , -1 
1990  NEXT  J 
2000  PENUP 
2010  PEN  0 

2020  INPUT  “Do  you  uiah  ■*  o  clear  th*  graphic*  display"'  *.Y  N>”,A  n*6$ 

2030  IF  An»6*»"N“  THEN  2060 
2040  SETGU 
2090  GCLEAR 

2060  INPUT  “Do  you  wish  to  run  i  ht  program  again  i.nth  a  new  j«t,  of  par  am«*  «r  $ ' 
Y/N> “ , Ana7F 

2070  IF  An*7«-“Y”  THEN  180 
2000  STOP 
2090  END 
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10  PRINT 
d  chsnncl 
20  PRINT 
d  Vg.  " 


30 
40 
S0 
60 
?0 
30 
90 
100 
110 
120 
130 
140 
130 
160 
170 
130 
190 
200 
210 
220 
2  30 
240 
230 
260 
270 
280 
290 
300 
310 
320 
330 
340 
330 


FLOAT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 


"The  n*nt  of  thi*  progr»»  is  BURCHL'.  It  computes  the  maximum  burl* 
potent i si,  Vz,  of  s  buried  N-ch*nn*l  device  if  w*  know  Os  end  Vg  OP" 
"it  finds  the  chennel  ehsrge,  Qs,  of  the  seme  device  if  w*  know  Vz  *n 


"TYPE  IN  THE  VALUES  OF  THE  PARAMETERS  Nd,  Ns,  Tox,  Xn,  Vg,  Vfb" 
“donor  c one ent r *t 1  on,  Nd,  in  units  of  E15  '  CENTIMETER'S" 

"Acceptor  concent  ret  1  on,  Ns,  in  units  of  E13  CENT  I  METER'  3  “ 

"Oxide  thickness,  T ox .  in  ANGSTROMS  <■  *  Micron*  *  10,600)“ 

“channel  depth,  Xn,  in  MICRONS  <■  •  Angstroms  *  IE-4  »  *  Cm  *  1E4)’ 

II  H 

"gat*  volttge,  V<j,  in  VOLTS" 


cm 

cm 

cb 

F  'c  m 


F  -cm"2 


flat-band  voltage,  Vfb.  in  VOLTS" 

INPUT  Nd,N», Tox.Xn, Vg, Vfb 
Nd«Nd* 1E1 3 
Na«Na«lE15 
Tox«Tox*lE-8 
Xn»Xn»lE-4 
Q«1 .602E-19 
Eo«8. 86E-14 
Eox«3 . 9 
E*i«ll.7 

Cox«Eo*Eox'Tox 

B«Q*Nd»(Nd«-Na) x.  2*Eo*Esi *N*> 

Vgpr im»Vg-Vfb 
PRINT  “  " 

PRINT  ■  ” 

PRINT  "  “ 

PRINT  “  » 

PRINT  "To  compute  Vz,  t..p*  'Vr  and  press  CONTINUE;  to  compute  0*.  type  0* 
'  *nd  press  CONTINUE;  to  stop,  press  STOP" 

360  INPUT  Sections 

370  IF  Sect ion*»"Vz"  THEN  400 

380  IF  Sect ionS«"Qs"  THEN  760 

390  1  - 

400 
410 
420 
430 

440 
450 
460 
470 
460 
490 
300 
310 
320 
330 
340 
330 
360 
370 
380 
390 
600 
610 
620 
630 


!  Vz<0s,Vgi 
INPUT  “Enter  'he 

Qs»Qs*lE-4  1 

XzsQ*  x  (  Q*Nd 1  1 

XmXn-X*  i 


■slue  of  0*  in  picoCOULOMBS 
C  x  c  m  '  2 
cm 

cm 


MICRON '2",Q* 


Xopr 1 m»Eo*Es i *< 2*8*X»CoxtQ*Nd>r <Cox*<2*B*Eo*E* 1 ~Q*Nd) > 

XI «Xopr i m-SOR  t Xopr im~2-2*Eo*E*i *<  B*X'2-Vgpr l*)/(2*B*Eo*Esi-Q*NJ>) 
PRINT  “X1«",X1 

700 


Angstroms 
Micron* 
pF  'um-'2 
pC/u»»2 


IF  X 1 < 8  THEN 
v*»B*  c  x-x  1  >-'2 
Tox»Tox»,E8  1 
Xn«Xn*lE4  • 

Cox»Cox»JE4  1 
Os»0s«lE4  • 

PRINTER  IS  7, 1 
PRINT  “Progr amt  BUPCML 
PRINT  “Nd  tC#-'-3)",Nd 
PRINT  “Ns  <Cm'-3>“,Na 
PRINT  “Tox  <Angstromt>“,Tox 
PRINT  “Qs  I  P  I  c  ocoul  ombs-' M  i  c  ron*2  )  “  ,  0* 
PRINT  “Xn  > Mi  crons) ", Xn 
PRINT  "Vg  (Volts)", Vg 
PRINT  “Vfb  < Vol t  s  > “ , Vfb 
PRINT  “Cox  (Pi cof sr sds 'Ml cron~2> " , Cox 


Section;  Vz(0s,Vg)" 
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640  PRINT  "V2  <Vol »»>■“, VI 

690  PRINTER  IS  16 

660  Tox»Tox#lE-8  1  c» 

678  Xo«Xn*lE-4  !  cm 
600  Cox»Cox»  IE-4  !  F/'C«^2 

690  GOTO  390 

700  INPUT  “XI  <  0  <not  allowed?  Do  vou  want  t  o  «nt«r  new  data?  <Y'N>",A* 
710  IF  A»»“Y"  THEN  60 

720  PRINT  “I  won't  accept  no  for  an  anawerl  Pre*»  CONTINUE" 

730  PROSE 
740  GOTO  60 

790  1  - 

760  I  8*CVz,Vg> 

770  INPUT  "Entar  the  value  of  Vz  In  VOLTS" , Vz 
780  fl«2»Eo*E»i  ',<2*B»Eo*E»i-Q«Nd> 

798  C»SQR< Vz^B)-Xn 

800  D»A»B 
810  EwRaVgprin 

820  F»Eo*E»i''<Cox*<2*B*Eo*£»i  -Q*Nd)> 

030  G“C+2*8*F*Cox*Xn+F*0*Nd 

040  H>E-D*Xn^2 

890  I  •<  1  -2*B*Cox*F )  '*<Q*Nd> 

860  J«2*B*F  »Cox*Xn+F »0»Nd 

870  K«2*B*F*C0X''<G*Nd' 

880  L«GA2-JA2-H 

890  M»2*K»J+2*G»I-2*D*Xn/(Q*Nd) 

900  N«I-'2-K'2+D/tO*Nd>''2 

910  P"M/N 

920  R«L'N 

930  D1  »c»P"-2-4*R 

940  PRINT  “Di»c«",Dt»c 

990  IF  D1»C>«0  THEN  980 

960  PRINT  "  Dt»c  -  P'2-4*R  it  <  0  » >S0R< Di  *c  >  i*  not  raal.  ENTEP  NEW 
DATA  “ 

970  GOTO  60 

980  0»*<  -P-SOR  <ti tc  > ) '2  t  C-'c»-'2 

990  Qsab»"ABS<Qs> 

1000  IF  Qsabs< ■ 1 . 39000000000E-1 7  THEN  Q»»0 

1010  Tox»Tox»lE8  1  Angatroaa 

1020  Xn«Xn»lE4  <  Hicrons 

1030  Cox»Cox*lE4  t  pFua-2 

1040  Q»-Q»»1E4  •  pC-  ua-'2 

1090  PRINTER  IS  7. 1 

1060  PRINT  "Program  8URCHL  Saction:  Q»<Vz,Vg>" 

1070  PRINT  "Nd  (C*-'-3)",Nd 

1080  PRINT  "Na  <C»~-3)",Na 

1090  PRINT  “Tox  < Ang»t ro*s > " , Tox 

1100  PRINT  "Q»  <P t cocou' 0»b»rMi cron-2> " , Q» 

1110  Q»»Q»*lE-4  '  C'c  «~2 

1120  PRINT  "Q*  <Couloi»b»/c»-2>",Q» 

1130  PRINT  "Xn  < M i cron* > “ , Xn 

1140  PRINT  “Vg  <Volt*>“,Vg 

1190  PRINT  "Vfb  «  Volt  *  1 ■ , Vf b 

1160  PRINT  “Cox  <Picofarad»^Nicron'-2>“,  Cox 

1170  PRINT  "Vz  <Volt»)»",Vz 

1180  PRINT  "Di»c»“,Dt»c 

1 190  Xl«Xopr  i  »-SQR'  Xopr  1 1«  '2-2*Eo*E*i  *<•  B*X  ~2-Vgpr  i  « > / <  2*B*Eo*E*  l  -0*Hd.»  > 

1200  PRINT  «X1«",X1 
1210  PRINTER  IS  16 
1220  IF  X1<0  THEN  1248 
1230  STOP 

1240  INPUT  “XI  <  0i  Do  you  wi»h  to  antar  naw  data?  < V- N > " ,Cho1 ca» 

1290  IF  Choi ca*a“Y"  THEN  60 
1260  STOP 
1270  END 
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"The  name  of  this  program  is  FGBCCD  tFloating  Gate'Buned  ch»nn«l " 
"CCD).  It  computesi" 

"  <l>  the  change,  Delvfg  <«Vqfg  -  Vofg),  in  the  floating  gate" 

■■  voltage,  Vfg,  of  a  buried  N-channel  CCD” 

“  the  partial  derivative  of  Delvfg  with  respect  to  Osignal" 

<3*  the  percent  change  in  the  partial  derivative" 

“  < 4 >  the  maximum  channel  potential,  Vz” 

“The  derivation  of  the  equations  may  be  found  in  the  21  May  1381“ 
"memo  by  James  Joseph  entitled  “Bur ied-Channel  Floating-Gate  Sensing" 
“With  Constant  Capacitive  Load'.  The  required  parameters  are" 
"entered  via  INPUT  statements.” 

"Note:  We  calculate  the  quantity  Delvfg  which  by  physical  argument .  ” 
“must  be  a  negative  number  since  Vofg  is  always  greater  than  Vqfg" 
“for  this  type  of  system.  However,  rather  than  label  the  Delvfg" 
"axis  with  negative  numbers,  we  plot  the  absolute  value  of  Delvfg" 

"and  label  the  vertical  axis  '  -  <  de  1 1  a)  Vf  g  ' .  Also,  Xm  must  be  a" 

"positive  number  since  it  is  simply  the  width  of  that  portion  of  the" 
"depletion  region  which  has  been  charge-neutr al i zed  by  the  signal" 
"charge.  Therefore,  although  it  is  defined  as  Qs-"Q*Nd)  in  the" 
"memo,  we  must  define  it  as  RBS CQs i gx C Rf g*G#Nd ) >  since  our  Qsig" 
“values  are  negative." 

"i press  CONTINUE)" 

PAGE 

“Is  this  progr am  being  run  on  the  HP9835A  or  HP9S43C"  1 A'C  1  ", Comput <r 


268 

IF  Computer*- 

278 

IF  Computer*” 

288 

GOTO  250 

298 

OPTION 

BASE  1 

300 

DIM  Os 

i gt 101  > 

rvq< 

1  0  1  )  ,  P  v 

q< 101  ) , 

318 

1 

328 

I 

330 

PRINT 

"The  fi 

348 

PRINT 

"Os lg  f 

350 

PRINT 

"via  a 

368 

PRINT 

“The  ma 

378 

PRINT 

"rev i sc 

380 

PRINT 

"Press 

390 

PAUSE 

488 

PRINT 

L 1N<  4 ) 

418 

Q-1.602E-19 

420 

Eo”8. 86E-14 

438 

Es  i  ”  1 1 

.  7 

440 

Eox«3. 

9 

458 

PRINT 

"ENTER 

460 

PRINT 

M  M 

470 

PRINT 

* 

488 

PRINT 

M  M 

490 

PRINT 

" 

500 

PRINT 

"  " 

510 

PRINT 

520 

PRINT 

530 

PRINT 

“ 

540 

PRINT 

..  « 

550 

PRINT 

M 

560 

PRINT 

"  " 

570 

PRINT 

" 

580 

PRINT 

..  .. 

590 

PRINT 

" 

688 

PRINT 

"  1. 

618 

PRINT 

628 

INPUT 

MN*,  Nd 

638 

•  Na»4 

640 

Na«Na* 

IE  1 5 

•  C 

1  F  '  c  m 


Na  in  units  of  IE  1 5  cm' -3" 

Nd  in  units  of  1E1S  cm-3" 

To  ir.  ANGSTROMS  <»  «  MICRONS  *  10,000>" 

Xj  in  MICRONS" 

Lfg  in  MICRONS  ■”  #  ANGSTROMS  *  IE-4  ■  *  CM  «  1E4>“ 
Wfg  in  MICRONS" 

Vofg  in  VOLTS” 

Vf b  in  VOLTS” 


650  !  Nd-5e 

660  Nd«Nd*  t  E 1 5 

670  I  Toxang»200 

680  Tox«Toxang*lE-8 

690  !  Xju»».3 

700  Xj-XjumME-4 

710  INPUT  "Lfg,  Wfg,  Vofg,  Vfb“ , Lf gum, UfguM, Vof g, Vf b 

720  •  Lfgua-40 

730  Lfg«Lfgu*»lE-4 

740  !  Ufgua-20 

750  Wf g»Wfgua» IE-4 

760  Afgu«2»Lf gu**MCgu» 

770  Afg»Lfg*Wfg 

780  >  Vofg-10 

790  !  Vfb—.9 

800  Cox»£o»£ox/Tox  •  F-'e»A2 

810  Coxpf»Cox»lE12  1  pF'ca~2 

320  INPUT  "Which  do  you  knout  LF  or  CS  in  pF''  (LF.'CS)  ,  Eni  *r  it*  " 

Value 

830  IF  An»*«“LF“  THEN  860 

840  IF  An*#«“CS“  THEN  910 

850  GOTO  820 

860  Lf  “Value 

870  !  Lf«.l 

880  C**Lf  *Cox«Af  g 

890  C*pf«C**lE12 

900  GOTO  950 

910  C»pf«Value 

920  !  C*pf». 25130 

930  Cs-C»pf*lE-l2 

940  Lf»C»x<Af g*Cox  ■ 

950  IF  Computer*»“A“  THEN  PRINTER  IS  7,1 
960  IF  Comput  er  #»“ C “  THEN  PRINTER  IS  0 
970  PRINT  ’’ Program t  FGBCCD" 

980  FLOAT  1 

990  PRINT  “Ha  <Ca~-3> ” » SPAv 1 2 > ; Na 
1000  FLOAT  1 

1010  PRINT  “Nd  (Cm~-3)"|SPA<12)jNd 
1020  FIXED  0 

1030  PRINT  "To*  ( flngitromt ) " ; SPA< 7 > ; Toxang 
1040  FIXED  1 

1050  PRINT  “Xj  <Micron*>”| SPA< 1 0) j X jum 
1060  FIXED  1 

1070  PRINT  "Afg  <  M  i  c  r  on*-'2 ) "  |  SPA<7 )  J  Af  gua2 
I08O  FIXED  1 

1090  PRINT  "Vofg  <Volt*.>"|SPA<10>|Vofg,,,Vfb  < Vo  1 1 a > " 5 SPA< 2 > ( Vf b 

1100  FLOAT  4 

1110  PRINT  "Cox  <Fxc«-'2)"|SPA<  10)»Cox 
1120  FIXED  2 

1130  PRINT  “LF  (loading  f ac t or >“ j SPAt 3> J Lf 
1140  FIXED  5 

1150  PRINT  “CS  <P 1  cofarad*  > “ 1 SPA<7> 1  C*pf 
1160  PRINTER  IS  16 

1170  INPUT  “Which  do  you  knout  Q»tg  In  pC  OR  number  of  e-"1  1 0  E>“,  Charge* 
1100  IF  Charget»”E“  THEN  1210 
1190  IF  Charge»«"Q“  THEN  1370 
1200  GOTO  1170 

1210  INPUT  “N,  Electron*! 1 >,  E 1 ee t r on* < N> “ , N, E 1 ec % ron* I o, E 1 ec l r onjhi 

1220  !  N* 1 1 

1230  1  E I ec t rontl o»0 

1240  !  Elec t r on»h I *6E6 

1250  0*1 gl opc«-El ec Ire -i* 1o»QelE12 

1260  0*1 ght pc •-£ lec t ronehl *0* IE  12 

1270  IF  Co*putert»“A“  THEN  PRINTER  IS  7,1 

1280  IF  Computer*»"C“  THEN  PRINTER  IS  0 

1290  FIXED  0 
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1  300  PRINT  "N" |  SPfl< 12) | N 
1310  FIXED  1 

1320  PRINT  “ Qs i g<  1 '  pC"  ,  Qs  i  g  1  opc  ,  “Q*  i  g<N>  pC",Qsighipc 
1330  FLOAT  3 

1340  PRINT  "El «c trons< 1 » “ , E 1 ec t r ons 1 o , “Elect rons  CN) ” , El ectronshi 
1330  PRINTER  IS  1 6 
1360  GOTO  1320 

1 370  INPUT  “N,8sig<l>  pC  ,Qsig<N>  pC  " , N, Qs i gl opc , Q» 1 gh i pc 

1330  •  Nall 

1390  1  Qsig1opc«0 

1400  1  Qsigh\pc«-1 

1410  E1ectronslo«-Qsiglopc*l£-127Q 
1420  El ec t ronthi «-Qsi ghi pc # IE- 12rQ 
1430  IF  Comput er  #« “  A“  THEN  PRINTER  IS  7,1 
1440  IF  Computer*- “C"  THEN  PRINTER  IS  0 
1430  FIXED  0 

1460  PRINT  “N“(SPA<.12>|N 

1470  FIXED  1 

1480  PRINT  “ Qs • g> 1 >  pC “  ,  Q» i  g '  opc  ,  "0* i  g< N )  pC.Qstghipc 

1490  FLOAT  3 

1300  PRINT  “Electrons'll*, Elect ron» I  o, "El*ctrons<N>“, Elect  ronshi 
1310  PRINTER  IS  16 

1  520  REDIM  Q  s  i  g  ■  N  > ,  Q  s  i  gpc  <  N>  ,E1  ec  *  rons<N> ,  Megael  ec  t  rons<  N  ■ ,  D*  1  ■  f  g>  N  > ,  Pd*r"qi  N  1 ,  P 
uq<  N> ,  Pder  per  cent  <  N>  ,  Vz'  N  > ,  Xu  ar<N),YoartN> 

1330  Qsi gl o»Qs i g 1 opc *  IE-  12 
1540  Qsi ghi *Q| i gh i pc  *  IE- 1 2 

1330  PRINT  "Qsig  '.pc  "Mi  I  I  ions  of  e- " ,  CHPSi  1  93.’  {  “Vfg  '  V  > '■ 

1360  PRINT  *  “ 

1570  Nn»N  1  so  we  c in  redefine  Nn  without  changing  N 

1580  A«Q*Nd*iNdtNi) ' <2*Eo*Esi *N»> 

1590  B»-17(Afg*Q»Nd>  1  PD  of  Xm  wrt  Qsig 

1600  Y-£o*Esi  (2*A»£o*Esi -Q*Nd> 

1610  C»-Y*2*A  1  PD  of  Xoprime  wrt  X* 

1620  Xo»Y»<2*A*Xj*Coxt-Q»Nd>'Cox 

1630  Xl»Xo-SQR<Xo~2-2eY«<A*Xj"2-Vofg*Vfb>) 

1640  D«Vofg-Q*Nd*Xl*AfgrCs 
1630  E«Q*Nd*AfgrCs 
1660  C-2*Y 

1670  FOR  I«1  TO  Nn  STEP  1 

1680  IF  Nn»l  TH~.N  Qs  i  g<  I  )  »0  =  i  gl  o 
1690  IF  Nn»l  THEN  1710 

1700  Qs i g< I )»Qs i gl oeCQs i ghi -Qs i g I o >*< l -1 ‘  ' C Nn-1 >  1  C 

1710  E I ec trons< I >*-Qs i g» I >rQ  i  e- 

1720  XmtABS vQs i gt I 1  >0 

1730  X  Dpr  ime«Y»'2*Ae<Xj-Xm>*Cox»Q*Nd.>7Cox 

1 740  F«Xopr i m»-2-2«Y*< AeCXj-Xm) -2+Vfb> 

1730  IF  F*D*G+E*G*XoprimetE-2*G~2 '4<0  THEN  4800 

1760  De  1  uf  g<  I '» “-£*< XI  -Xopr  i  mt-E*Gx2*SQR <f* D*G*E*G*Xopr  i  me*E  2*G  '2 ' 4 )  >  1  V 

1770  Vqfg»Vofg-»Del  vfg<  I ) 

1780  Xlpr ime*Xopr  i  me-SQR'  Xopr  i  me~2-2*Y*t  A*(  X  j -Xm )  ^2-Vgf  g+Vf  b  )  > 

1798  IF  Xlpr i me<  0  THEN  4930 
1800  FIXED  3 

1810  Qsi  gpc  1 1 1  «Qs  i  g>.  I  >•  IE  1 2 

1820  Meg*e1  ectronst  I  ’•»£  1  ectronK  I  >*lE-6  '  mil  liens  of  e-.need  *rr  aw  for  plotting 
1830  PRINT  Qs  i  gpc  <■  I ' ,  Meg«e  1  ect  ronK  I  > ,  De  I  vfg<  I  > 

1840  NEXT  I 

1850  PRINT  "End  of  Output " 

1860  Gpeum2"Qs i gpc 'Nn > ^ Af gu»2 
1870  Kelectronsum2«El eetrons<Nn>xAfgu*2*lE-3 
1880  IF  Computer**" A"  THEN  PRINTER  IS  7,1 
1890  IF  Computer*«"C"  THEN  PRINTER  IS  0 

1900  PRINT  " Qs i g<  Nn  >  Af  g  ‘'maximum  allowable  signal  charge  per  gate  area)  pC'um"2 
-|Qpcum2 

1910  PRINT  "Thousand*  of  e-  per  urn "2  c orrespond 1 ng  to  Gsig'Nn)  “ I Ke 1 ec t r onsum 

2 

1920  PRINT  L I N  <  2  >  1  seperates  hardcopy  output 
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1930  PRINTER  IS  16 
1940  ! 

1950  > 

1960  PRINT  “The  »»con«  wet  ion  of  th«  program  calculates:" 

1970  PRINT  "  t 1 •  the  partial  dtrivtt iv«  of  Delvfg  with  rtipt't  to  Osig 

1980  PRINT  “  i2>  the  ptrctnt  change  in  the  partial  dtriuitiv*" 

1990  PRINT  “  <.3>  the  noxiaum  channel  potential,  Vz" 

2000  PRINT  “using  tha  pariMttrt  and  0*ig  array  (>o*  th*  previous  faction. 
2010  PRINT  “Press  CONTINUE  to  Access  the  second  section  of  the  program. " 

2020  PAUSE 
2030  PRINT  LIN<6> 

2040  PRINT  "Os  i g  pC “ , "S 1  ope < xl£  1 1  V'C>“,  slope  change “,“ Vz " 

2030  FOR  1*1  TO  Nn  STEP  1 

2060  X»*ABStQsi g< 1 )*B) 

2070  Xopr  i  »e*Y* < 2*A* < X J -Xm ) #Cox+0*Nd> ' Cox 
2000  F-Xopr  i«e''2-2*Y*<A*<Xj-X»)-2+Vfb> 

2090  H«CD*G+F>e£A2+G*Xoprt!ee*E''3+G~2*E"4/4 
2100  Vqf g*Vof g+De I  of g< I  > 

2110  Xlpr  i  me*Xopr i  me-SQR<Xopr  i  me ^2-2* Yet  A#<X  j-X» »  "2-Vqf  g*Vf  b  >  > 

2120  Pdervqt  I  >*E»B*C- .  5»E'2»B*C.  H-'.  5*<E*G*2*Xopr  l  »e-2*<Xj-Xm  >  > 

2130  Pvq< I >*Pdervq< I >*1E-1 1  I  need  array  for  plotting 

2140  p  derperc  ent  C I  )  •  <  Pdervqt  1 ' -Pder  vq<  1 ) > -Pderoq'  1)*10Q 
2150  Vz<  I  )*A*<Xj-Xm-Xlpr  i  «e)A2 
2160  IF  Vz<  I XVqf g  THEN  3090 
2170  FIXED  3 

2180  PRINT  Qs  i  gpc  >  I  >  .  P»q>.  I  • ,  Pder  percent  c  I  > ,  Vz<  I  ) 

2190  NEXT  I 

2200  PRINT  “End  of  Output" 

2210  PRINT  LINC2) 

2220  < 

2230  • 

2240  INPUT  "Do  you  desire  a  graphical  output  of  this  data'  •  '<  N  > " ,  Ans  1  * 
2230  IF  Bnsl**"Y“  THEN  2280 
2260  IF  Anel*»“N"  THEN  3040 
2270  GOTO  2240 

2200  PRINT  “You  may  choose  which  quantity  you  want  on  the  horizontal  ais: 
2290  PRINT  “  i>  Qsig  in  pC  CODEt  Q" 

2300  PRINT  "  ill  Number  of  e-  CODEt  E" 

2310  PRINT  "Enter  the  code  of  the  variable  for  the  X  axis." 

2320  INPUT  "XCODE't"  ,  Xcode* 

2330  PRINT  PAGE 

2340  PRINT  "You  may  choose  which  quant ity  you  want  on  the  vertical  axis!" 

2330  PRINT  "  i)  Delta  Vfg  CODEt  DV" 

2360  PRINT  “  it)  Partial  Derivative  CODEt  PD" 

2370  PRINT  “  iii)  X  change  tn  PD  CODEt  PC" 

2380  PRINT  "  iv)  Vz  CODEt  VZ“ 

2390  PRINT  "Enter  the  code  of  the  variable  for  the  Y  a. is." 

2400  INPUT  "YCODE?", Ycode* 

2410  PRINT  PAGE 

2420  PRINT  “XCODE  »  “ j Xc ode*, “YCODE  »  "t Ycode* 

2430  FOR  1*1  TO  Nn  STEP  1 

2440  IF  Xcode*»"Q“  THEN  Xuar < I ) »Qs i gpc < I > 

2430  IF  Xcode*«"E"  THEN  Xuar < I >*Megae l ectronsc I ) 

2460  IF  Ycode**“DV“  THEN  Yvart I )»De 1 vf gC I > 

2470  IF  Yeode*«"PD"  THEN  Yvar < I >*Pvq< I > 

2480  IF  Ycode*»“PC"  THEN  Yvar < I >*Pderperc*nt < I > 

2490  IF  Ycode**“VZ“  THEN  Yvar  <  I ) *Vz <  I  > 

2300  Xvarmi n*ABS<Xvar < 1 ) ) 

2510  Xvarmax* BBS <  Xvar  < 1 ) ) 

2320  Yvarmtn»A8S< Yvar< 1 ) ) 

2330  Yvarmax*ABS< Yvar < 1 ) ) 

2340  IF  1*1  THEN  2590 

2330  IF  ABStXvar < I >  >< ABSCXvar 1 1- 1 >  >  THEN  Xvarm i n*ABS t Xvar < I > > 

2360  IF  BBS  <  Xvar  < I ) ) >ABS<Xvar < 1-1 ) )  THEN  Xvarmax-ABSc Xvar < I ) ) 

2570  IF  BBS < Yvar < I >  X BBS < Yvar < I  —  1 >  >  THEN  Yvarm i n*ABS < Yvar < I ' > 

2580  IF  ABS<  Yvar  i  I  '  '  >ABS 1  Yvar  <  I  -  I  )  >  THEN  Yvar  max*ABS  <  Yvar  <  J  " 
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2590  NEXT  l 

2600  PRINT  "Xuarmirt*  n  ;  Xv»rm  n,  “Xuaraax*  “(Xuaraax 
2610  PRINT  “Yvar»in»  “ j Yuaratn, “Yuar»4x«  *|Vv(r»tx 

2620  INPUT  "Do  you  want  lh«  graphical  output  on  th«  9872ft  or  9845C7  <  72ft -45C  > " ,  ft 

n»2* 

2630  IF  An»2*«"?2A"  THEN  2660 
2640  IF  An»2*«"45C"  THEN  3650 
2650  GOTO  2620 

2660  1  9872ft  PLOTTER  ROUTINE - 

2670  PRINT  "BE  SURE  THERE  IS  PAPER  ON  THE  PLOTTER  than  praa*  CONTINUE" 

2680  BEEP 
2690  PAUSE 
2780  PRINT  PAGE 

2710  INPUT  "Do  you  n««d  to  draw  n# w  axes7  (YxN)",An»3* 

2720  IF  ftn*3««“Y“  THEN  2750 
2730  IF  ftn»3*»"N"  THEN  3440 
2740  GOTO  2710 

2750  IF  Co»put«r*»"C“  THEN  PLOTTER  IS  13, “GRAPHICS" 

2760  IF  Co*put«r*«"C“  THEN  PLOTTER  13  IS  OFF 

2770  PLOTTER  IS  7,5,"9072A" 

2780  PLOTTER  7,5  IS  ON 

2790  LIMIT  20,260,20,190  <  mm 

2800  LOCATE  18,97,10,97  »  GDU» 

2810  X» l n* INT» Xv ar m i n • 

2820  Xm*x»  INT(  X"»r*i«x  1 
2830  Yin  i  r>«  INT  <  Yi-arm  i  n  ' 

2840  Ymax  •  I  NT  i  Yi<*r  ma  ■  1  ♦  1 
2850  SCALE  X»in,X»»«,Ynin,Y«uy 
2860  PEN  1  'black 

2870  LINE  TYPE  1 

2880  AXES  1'10, I- lO.Xain, Yarn, 10, 18,2 

2890  I  - - -  AXIS  LABELING  ROUTINE 

2900  '  -  X  axil  till#  i 

2910  DEG  > 

2920  MOVE  X*i  rvt .  5»  <Xi»ax-Xi»i  n ' ,  Y«  i  n-.  067*<  Y»ax~Y«*i  n>  ' 

2930  CSI2E  5, .5  ! 

2940  LORG  6  1 

2950  LOIR  8  ' 

2960  IF  Xcod#*»“Q"  THEN  2980  ! 

2970  IF  Xcod«»*"E“  THEN  3000  ' 

2980  LABEL  "Qaig  <pC""  > 

2990  GOTO  3020  ' 

3000  LABEL  "Electrons  (xlE+6"'"  > 

3010  GOTO  3020  ! 

3020  ! -  X  axi*  tiu»b*r*  ! 

3030  CSI2E  3,  .5  ' 

3040  LORG  6  < 

3050  MOVE  0,  Y»in-.03*<.Y»ax-Yi»in>  > 

3060  LABEL  USING  "K"|“0"  ! 

3070  FOR  A* , 5  TO  X»a*  STEP  .5  ! 

3080  MOVE  A,Yain-.03»<Y»4x-Yi»in>  > 

3090  IF  Xcod«*»"Q"  THEN  3110  ! 

3100  IF  Xcod«»«"E"  THEN  3130  ' 

3110  LABEL  USING  "K"|"-“jA  ' 

3120  GOTO  3150  < 

3130  LABEL  USING  "K‘|A  ! 

3140  GOTO  3150  < 

3150  NEXT  A  ! 

3160  i  -  V  axil  till*  ! 

3170  MOVE  X*it it- .  07* < X*>ax-X» in),Ykint, 5»< V*>ax-Y*it n>  ' 

3180  CSI2E  4, .5  i 

3190  LDIR  90  i 

3200  IF  Yc od€*« "Dv"  THEN  3240  ! 

3210  IF  Ycod#*-"PD"  THEN  3260  ! 

3220  IF  Yco<J«*«“PC"  THEN  3280  ! 

3230  IF  Yeotj«*«"VZ“  THEN  3340  ! 


3240  LABEL  * - " ) CHR* U 98) > " Vf CV>“  i 

3250  GOTO  3350  < 

3260  LABEL  “SLOPE  CxlE+lt  VxC>“  i 

3270  GOTO  3350  i 

3280  IF  Pd*rp*rc*nt<Nnx2>>0  THEN  3300  i 

3290  IF  Pd*rp*rc#m  (Nn  '2><0  THEN  3320  i 

3300  LABEL  “5j  CHANGE  IN  SLOPE”  ! 

3310  GOTO  3350  ! 

3320  LABEL  *-  5;  CHANGE  IN  SLOPE"  ! 

3330  GOTO  3350  > 

3340  LABEL  “Vr  <V>“  i 

3350  ! -  V  axis  nuabtrs  ! 

3360  CSI2E  3,  .5  < 

3370  LORC  8  < 

3380  LOIR  0  > 

3390  FOR  A«0  TO  Ymax  STEP  1  > 

3400  MOVE  Xm i n- . 083*  < Xmax-Xm i  n) ,  A  ! 

3410  LABEL  USING  “K,X"tA  > 

3420  NEXT  A  ! 

3430  !  - 

3440  PEN  4  i  r*d 

3450  PLOT  ABS<Xuar< 1 ) ), ABS< Yusri 1 ) ),-2 
3460  FOR  J«2  TO  Nn  STEP  1 
3470  PLOT  ABS<  Xuar (J>>, ABS ( Yuar <  J>  > , -1 
3480  NEXT  J 

3490  !  -  locus  1  ab*  1 

3500  PEN  1  i  black 

3510  CSIZE  2.5, .5 
3520  LORG  2 

3530  MOVE  ABS< Xuar <Nni)+.01,ABS<Yuar<Nn)> 

3540  !  LABEL  USING  "X, K " ; “ t  ox«“ j  Toxang; CHRt ( 208 ) ; “  LF«" ; Lf 
3550  LABEL  “LF-  "  |  Lf 
3560  I  - 

3570  PRINT  "LETTER  command  actioat«d.  Pr«ss  CONTINUE  to  deacttuat*.  DO  NOT" 
3580  PRINT  “LEAVE  THE  PLOTTER  IN  LETTER  MODE  <p*n  will  dru  out  > ■ ' » 

3590  CSIZE  2.5, .5 
3600  LETTER 
3610  PRINT  PAGE 
3620  PENUP 
3630  PEN  0 
3640  GOTO  4710 

3650  !  9845C  CRT  GRAPHICS  ROUTINE - 

3660  PRINT  PAGE 

3670  INPUT  “Do  you  n*«d  to  draw  new  ( Y.'N  >  “ ,  Ans4» 

3680  IF  Ans4*»"Y"  THEN  3710 
3690  IF  Ans4*«"N“  THEN  4410 
3700  GOTO  3670 
3710  PLOTTER  IS  7,5,“9872A“ 

3720  PEN  0 

3730  PLOTTER  7,5  IS  OFF 
3740  PLOTTER  IS  13, “GRAPHICS" 

3750  IF  Comput*rS-“C“  THEN  GRAPHICS 
3760  LIMIT  0,125,0,125  •  mm 

3770  LOCATE  10,99,10,99  I  GDUs 

3780  Xmin«lNT<Xuarmin> 

3790  XmaxmINT<Xuarmax)*l 
3800  Ymi n«INT< Yuarmi n) 

3810  YmaxmlNT < Yuarmaxl+l 
3820  SCALE  Xmi n, Xmax, Ymi n, Ymax 
3830  PEN  1  I  whit* 

3840  LINE  TYPE  1 

3850  AXES  1X10, 1x10, X*tn, Ymin, 10, 10, 1 

3860  I  -  AXIS  LABELING  ROUTINE 

3870  I -  X  axis  tit  1*  ! 

3880  DEG  i 

3890  MOVE  Xmi n*. 5»< Xmax-Xmi n) , Ymi n-. 067*< Ymax-Ymi n > i 
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3990  CSIZE  5, .5  I 
3910  LORG  6  ! 
3920  LOIR  0  < 
3930  IF  Xcod**"“Q"  THEN  3990  ! 
3940  IF  Xeod4*«“E"  THEM  3970  ! 
3990  LABEL  “Qsig  (pC>“  • 
3960  GOTO  3990  • 


3970  LABEL  "E1«ctrons  <xtE*6>"  ! 

3900  GOTO  3990  ! 

3990  I  -  X  axis  nuubtri  • 

4000  CSIZE  3,  .5  < 

4010  LORG  6  • 

4020  HOVE  0,Yi»in-.93*<Y*i*x-Y«in>  ! 

4030  LABEL  USING  "K"l“0"  ! 

4040  FOR  A-.9  TO  X«4x  STEP  .9  ! 

4090  MOVE  A,  Y»i  n-.  93*<  Y»*x-Y»t  r>)  I 

4060  IF  Xcod«*«"Q”  THEN  4000  ! 

4070  IF  Xcod«*""E"  THEN  4100  ! 

4000  LABEL  USING  “K“|“-“|A  • 

4090  GOTO  4120  • 

4100  LABEL  USING  "K"|A  I 

4110  GOTO  4120  ( 

4120  NEXT  A  ! 

4130!  -  Y  txil  titlt  ! 

4140  HOVE  X*  1  n- .  07*<X«i*x~X«)  n>  ,  Yioi  r>+.  9*<  Y»»x-Yi«i  r>)  ! 
4190  CSI2E  4, .9  ! 

4160  LDIR  90  ! 

4170  IF  Ycod**»”OV”  THEN  4210  ! 

4100  IF  Ycod**»“PD“  THEN  4230  ! 

4190  IF  Ycod«*“"PC"  THEN  4290  ! 

4200  IF  Ycod**«“VZ"  THEN  4310  ! 

4210  LABEL  “ - “ ; CHR* < 1 98 ) i "Vf g  <V>"  • 

4220  GOTO  4320  ! 

4230  LABEL  "SLOPE  <xlE*ll  VxC>"  *  ! 

4240  GOTO  4320  ■ 

4290  IF  Pd*rp*rc*m <Nnx2)>0  THEN  4270  ! 

4260  IF  Pd€rp*rc«fU<Nn-'2><0  THEN  4290  ! 

4270  LABEL  CHANGE  IN  SLOPE"  ! 

4200  GOTO  4320  ! 

4290  LABEL  *  CHANGE  IN  SLOPE"  ! 

4300  GOTO  4320  < 

4310  LABEL  “Vz  <V>"  ! 

4320  !  - -  Y  ixil  numbers  ! 

4330  CSIZE  3, .9  < 

4340  LORG  0  ! 

4390  LOIR  0  ! 


4360  FOR  A-0  TO  Y»*x  STEP  1  t 

4370  MOVE  X»i n- . 003»<  X»4x-Xai n  > , A  • 

4300  LABEL  USING  "K,X"|A  ! 

4390  NEXT  A  ! 

4400  !  - 

4410  PEN  2  !  r#d 


4420  PLOT  ABSCXumr < l ) > , AB$< Yu*r < 1 ) > , -2 
4430  FOR  J-2  TO  Nn  STEP  1 
4440  PLOT  ABS<Xv*p<3)>,ABSCYv«r<J>>,-l 
4490  NEXT  J 

4460  I  - -  locus  1 4b • I 

4470  CSIZE  2. 9,. 9 
4400  LORG  2 

4490  MOVE  ABS(Xu*r<Nn))*.01,ABS<Yu*p<Nn>> 

4900  •  LABEL  USING  “X,K*|*tox«*| Toxangt CHRP < 200 > | "  LF»“ | Lf 
4910  LABEL  "LF»  "|Lf 
4920  I 

4930  PRINT  “LETTER  eo«t«i*nd  *ctw*t*d.  Pros*  CONTINUE  to  d*4ctlv«to. 
4940  CSIZE  2. 9,. 9 
4990  LETTER 
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4560  PRINT  PAGE 
4570  PENUP 
4580  PEN  0 

4590  INPUT  "Do  you  wish  to  get  a  hardcopy  printout’  'Y  N>",An*5Z 

4600  IF  Ans5**“Y”  THEN  4630 
4610  IF  Ans5*-"N“  THEN  4650 
4620  GOTO  4590 

4630  IF  Computer*»“C"  THEN  DUMP  GRAPHICS 
4640  !  DUMP  GRAPHICS  *7,1 

4650  INPUT  “Do  you  wish  to  clear  th«  graphic*  display’  c  Y-- N>  “  ,  Ans6< 

4660  IF  An*6S«"Y“  THEN  4690 
4670  IF  Ans6*«"N“  THEN  4710 
4680  GOTO  4650 
4690  SETGU 
4700  GCLEAR 

4710  INPUT  "Do  you  ui*h  to  run  th*  program  again  with  a  ntu  s* t  of  parameters’  < 
YrN) “ , Ans79 

4720  IF  Ans7*«“Y"  THEN  410 
4730  IF  Ans7*«"N"  THEN  5750 
4740  GOTO  4710 
4750  PRINT  PAGE 

4760  PRINT  Th*  progi  a*  ha*  ended.  To  r**tart  the  program,  press  PUN.  - 


4770 

4780 

4790 

4800 

4810 

4820 

4830 

4840 

4850 

4860 

4870 

4880 

4890 

4900 

4910 

4920 

4930 

4940 

4950 

4960 

4970 

4980 

4990 

5000 

5010 

5020 

5030 

5040 

5050 

5060 

5070 

5080 

5090 

5100 

5110 

5120 

5130 


PRINT  L1NC5) 
STOP 


Nn«I-l 

PRINT  “End  of  Output" 

PRINT  “A  calculation  ha*  baen  mad*  in  the  Kicuimn  of  » he  program" 

PRINT  "resulting  in  a  negative  number.  The  next  statement  request*  the" 
PRINT  “square  root  of  that  negative  number.  This  t»  not  a  valid  operation" 
PRINT  "*o  the  computation*  have  been  stopped  and  the  r**ul*»  computed" 

PRINT  "up  to  thi*  point  have  been  printed  on  the  screen  as  the"  normally" 
PRINT  "uould  be.  The  remainder  of  the  program  may  be  executed  as  if  thi*” 
PRINT  “had  not  occurred  EXCEPT  that  now  the  program  will  only  handle  the" 
PRINT  "results  obtained  so  far  and  not  th*  N  point*  requested  in  t he  INPUT" 
PRINT  "statement.  Th*  request  for  the  square  root  of  a  negati"*  number" 
PRINT  "indicates  that  th*  operator-entered  data  is  i  ncompat  ab  I  e  with  the" 
PRINT  "physical  situation  represented  by  th*  program  equations." 

GOTO  I860 


Nn-I-I 

PRINT  “End  of  Output" 

PRINT  "A  calculation  has  been  made  in  th*  execution  of  the  program" 

PRINT  "resulting  in  a  negative  value  of  Xlprim*.  This  is  not  a  physically" 
PRINT  "realizable  situation  so  th*  computation*  have  been  stopped  and  the" 
PRINT  "results  computed  up  to  thi*  point  have  been  printed  on  th*  screen" 
PRINT  "as  they  normally  would  be.  Th*  remainder  of  th*  program  may  be  “ 
PRINT  "executed  as  if  this  had  not  occurred  EXCEPT  that  now  th*  program" 
PRINT  "will  only  handl*  th*  results  obtained  so  far  and  not  th*  N  points" 
PRINT  "requested  in  the  INPUT  statement.  Xlprim*  will  in  general  be  '  0." 
PRINT  "It  will  be  »  0  when  th*  surface  depletion  region  collapses  at" 

PRINT  "maximum  channel  charge." 

GOTO  1860 

i  - 

PRINT  PAGE 
PRINT  "Vz  <  Vqfg" 

PAUSE 
GOTO  2170 


Program  FGBCCD 
Na  (Ca-'-S) 

Nd  (Ca^-3) 

Tox  (Angstroas) 

Xj  (Microns) 

Afg  (Microns~2> 

Vofg  (Volts) 

Cox  (F/ca-'E) 

LF  (loading  factor) 
CS  (Picofarads) 

N  181 


4.8C+13 
4. 2E+16 
1188 
.38 

1093.3 

13.8 

3. 1413E-08 
.43 

. 13448 


Vfb  (Volts)  -.3 


Qsig(l)  pC 
Electrons(  1 ) 
Qsi  g(Nn)/'flf  g 
Thousands  of 


8.0  Qsig(N)  pC  -.3 

8. 000E+00  Electrons(N)  3.121E+86 

(aaxiaum  allowable  signal  charge  per  gate  area)  pC^ua^ 
e-  per  um2  correspond! ng  to  Qsig(Nn)  2.833 


Figure  Bl.  HP9845C  Printout  of  the  Floating-Gate  voltage  Response 
for  *3022  Matrix  Processor  Chip  (Paraneters  are  listed 
at  the  top.) 
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Program  FGBCCD 
Na  (Cn*-3) 

Nd  <C«-'-3> 

Tox  (fingstrons) 

Xj  (Microns) 
flfg  (Microns'^) 

Vofg  (Volts) 

Cox  <F/'c»*2) 

LF  (loading  factor) 
CS  (Picofarads) 

N  101 

Qsig(l)  pC 
El tctrons( 1 ) 

Qsi g(Nn)/fifg  (naxin 
Thousands  of  a-  par 


4.0E+15 

4.2E+16 

1100 

.56 

1093.3 

13.0 

3. 1413E-08 
.45 

. 15440 

0.0 

0.090E+00 

m  al 1 owabl a  signal 
unA2  corrasponding 


Vfb  (Volts)  -.5 


Qsig(N)  pC  -.5 

Elactrons(N)  3. 121E+86 

charga  par  gata  araa)  pC^um2 
to  Qsig(Nn)  2.855 


Figure  B2.  HP9845C  Printout  of  the  Deviation  from  initial  Sensitivity 

for  the  #3022  Chip  (Parameters  are  listed  at  the  top.) 
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APPENDIX  C 


SOFTWARE  USER’S  GUIDE 

This  appendix  is  a  detailed  user's  guide  for  the  Exerciser  unit  delivered  to 
RADC.  The  following  topics  are  covered: 

o  Procedure 

o  Hierarchy 

o  Memory  allocation 

o  Startup  sequence 

o  Description  of  functions 

The  RADC  Exerciser  is  supplied  with  software  resident  in  ROM  and  is  normally 
to  be  '<sed  with  Honeywell-supplied  primitives.  However,  modifications  to  the 
set  can  be  implemented  by  hookup  to  an  MDS  (8086  compatible) . 

A  complete  listing  of  PLM  operating  software  can  be  obtained,  if  desired, 
direct  from  Honeywell  SRC  (Dr.  P.C.T.  Roberts,  612/378-4992). 
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DISCUSSION 


Procedure  Hierarchy 

The  software  is  composed  of  a  collection  of  procedures.  A  procedure  may  have 
other  procedures  embedded  within  it  or  it  may  use  procedures  that  are  externally 
defined.  To  show  the  overall  relationship  between  procedures  a  hierarchy  is 
defined.  The  hierarchy  scheme  used  here  is  based  on  whether  or  not  a  procedure 
is  used  by  any  other.  If  so,  it  is  placed  on  a  level  belbw  it  in  the  hierarchy. 
If  it  is  used  by  two  different  procedures,  it  is  defined  to  be  a  utility  proce¬ 
dure.  These  are  referenced  by  their  number  only  and  are  included  as  a  separate 
list  below  the  hierarchy  list.  This  should  leave  one  procedure  at  the  outer 
level  of  the  hierarchy;  however  another  type  of  procedure,  namely  the  interrupt 
procedure,  is  also  placed  at  the  outer  level.  Interrupt  procedures  are  not 
referenced  by  any  other  procedure  directly  but  rather  by  the  hardware  mechanism 
that  controls. the  interrupt;  therefore,  they  are  placed  at  the  outer  level. 

With  the  above  definitions  the  following  hierarchy  results: 

PIPOMO  (1,2,4,7,8,9,14,20) 

COMPUTEJG  3RANCH_ADDR 
SET  BAUD  RATE 
iniT  intPrrupt_cont 
SELECT  FROM  MENU 
GREY_SCALE  ~ 

TEST  PATTERN  (1,2,3,6,7,8,9,12) 

EDIT  ELEMENT 

ELEMENT  XY 
UPDATE  SCREEN 
INIT  INPUT  ARRAY 
DIGITIZE  FRAME- (9)  ~ 

TRACK  GATE  (2, 3, 5,6, 7,8,9) 

MOVE_GATE 

RESTORE  VIDEO 
SAVE  VlCEO  AND  DRAW  EDGE 
HOME  GATE  -  -  -  - 

VIDEC_T0_PIP 

XFER_VIDE0  T0_PIP  ARRAY 
XFER  PIP  ARRAY  TO  VIDEO 
THRESHOLD  ~  “ 

DARK  THRESHOLD 
LIGHljHRESHOLD 
BINARY  THRESHOLD 
QUICK  TG  MEMORY  TEST  (8,10711) 

~  "  TG  ADDRESS  FORMATTING 

Print  error  addr 

QUICK  10  MEMORY  TEST  (8,10,11)  “ 

-  10  AOCRESS  formatting 

Print  error  addr 

TG  UTILITIES  (16,17,18,19720)  - 
GET  ADDR 
10  ffUFFER_BYTE 
TG-S INGLE  STEP 


TG  HALT 
7G”SEG_A0DR 
7B_JAM  AOOR 
TG  VERIFY 
TG_SAVE 

HALT  TG  INTERRUPT 
DIGITIZE  AT  INTERRUPT 


Utility  Procedures 


1)  RG8256  INIT 

2)  INTERPRET_PRIMITIVES 

3)  IOB  MEMIO 

4)  BRANCH  AOOR 

5)  EXECUTl_PRIMITIVES 

6)  PROMPT  AND  GET  CHAR 

7)  SIO  OUTCHAR 

8)  SIO  OUT  STRING 

9)  SIO  GET  CHAR 

10)  SIO_OUT_BYTE 

11)  SIO  OUT  WORD 

12)  SIO  ASCll_TO_HEX 

13)  SIO  CRLF 

14)  SIO  MS  DELAY 

15)  SIO'GET  STRING 

16)  TG  RESET 

17)  TG_RUN 

18)  TG  MACRO_SELECT 

19)  TG  MEMIO" 

20)  TG_LOAO 


Memory  Allocation 


The  8086  CPU  address  is  20  bits  long  giving  a  range  of  0  to  FFFFF  hexidecimal 
(0  to  1,048,575  decimal).  The  single  board  computer  (SBC)  has  been  configured 
such  that  RAM  is  located  at  addresses  0  to  7FFF  (0  to  32,767  decimal)  and  EPROM 
is  located  at  addresses  F8000  to  FFFFF  (1,015,808  to  1,048,575  decimal).  The 
following  shows  the  starting  location  for  the  Intel  supplied  monitor  program 
and  scratchpad,  the  PIP  wavefont!  data,  the  program  data,  and  the  program  code 
and  constants. 


FIRST  ADDRESS 

HEX  (DECIMAL)  DESCRIPTION 

0  (0)  Intel  957A  Monitor  RAM 

700  (1,792)  Program  Data 

2,000  (8,192)*  Program  Code 

6,000  (24,576)*  Constants  and  P.IP  Waveforms 

FEOOO  (1,040,384)  Intel  957A  Monitor  Code  and  Constraints 

♦The  Program  Code  and  Constants  when  in  EPROM’s  are  located  starting 
at  address  F8000  (1,015,303).  Constants  and  PIP  Waveforms  when  in 
EPROM's  are  located  starting  at  address  FCOOO  (1,032,192). 
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Startup  Procedure 


From  RAM: 


When  the  program  is  to  be  executed  from  RAM  the  following  steps  are  per- 
formed: 

The  following  should  be  turned  on:  MOS  computer,  SBC,  camera,  monitor 
and  the  TI  terminal  (if  it  is  to  be  used).  The  following  steps  should 
be  performed  in  the  order  shown: 

1)  Put  disk  "PIP  MASTER"  in  drive  1  and  disk  "SYSTEM. OBY"  in  drive  0 

2)  Press  RESET  on  MDS. 

3)  Press  RESET  on  the  SBC. 

4)  Type  SU8MIT  :F1:PIP 

5)  When  the  MDS  displays  :F1:TAL.K,  press  RESET  on  the  SBC. 

The  program  has  now  been  down  loaded  to  the  SBC  RAM.  If  the  MDS  terminal 
is  to  be  used  two  control  U's  must  be  typed  on  the  MDS  keyboard.  Now  the 
command  G200:2  can  be  typed  to  begin  program  execution. 

If  the  TI  terminal  is  to  be  the  I/O  device,  its  cable  must  be  put  on  the 

J2  connector  in  place  of  the  download  cable.  The  SBC  must  be  reset  and 

two  control  U's  typed  on  the  TI  keyboard.  Now  the  command  G200:2  can  be 
typed  to  begin  program  execution. 


From  EPROM: 


When  the  program  is  to  be  executed  from  EPROM  the  following  steps  are  per¬ 
formed: 

Turn  on  the  S8C  chassis,  camera,  monitor  and  the  TI  terminal  or  the  MDS 
computer.  The  cable  from  the  TI  terminal  or  the  MDS  computer  (whichever 
device  is  to  be  used  for  I/O)  is  to  be  connected  to  the  J2  connector  on 
the  SBC  card. 

If  the  MDS  computer  is  to  be  used  the  disk  "PIP  MASTER"  must  be  placed 
in  drive  1  and  disk  SYSTEM. DBY  in  drive  0.  The  following  must  be  typed: 

:F1:BAUD  112 
: FI: TALK 

Now  with  either  terminal  the  monitor  is  Invoked  by  pressing  two  control 
U's.  The  interrupt  vector  pointer  locations  must  now  be  loaded  as  follows: 


type: 


SW  AO,  <er> 
2184, 

F800, 

21DB, 

F800, 

087E, 

F800, 

08DA, 

F800  <cr> 


The  program  is  then  started  by  typing  GF800:2. 


Description  of  Functions 

When  the  program  is  started,  the  following  is  printed: 

LIST  OF  AVAILABLE  FUNCTIONS 

G  *  DISPLAY  GRAY  SCALE 
P  *  OI SPLAY  TEST  PATTERN 
D  =  DIGITIZE  ONE  FRAME 
T  *  TRACK  GATE 
M  -  TG.  MEM.  TEST 
0  *  10.  MEM.  TEST 
U  -  TG  UTILITIES 
I  =  INTERPRET  ANO  EXE 
L  »  LIST  FRAME  MEMORY 

SELECT  FROM  LIST 

Each  of  these  functions  will  be  described  in  the  following  paragraphs.  The 
following  conventions  are  used  in  command  syntax  descriptions  (note  that  "A" 
and  "8"  represent  arbitrary  ASCII  character  strings): 

[A]  indicates  that  "A"  is  optional 

[A  b}  indicates  that  either  ‘A"  or  “B"  must  be  selected 

< A  >  indicates  that  "A"  is  an  ASCII  hexadecimal  constant 

<cr>  indicates  that  a  carriage  return  is  entered 

U  indicates  that  a  blank  or  space  is  entered 

DISPLAY  GRAY  SCALE: 

This  function  displays  on  the  last  16  lines  of  the  monitor  a  scale  of  16  gray 
shades  ranging  from  black  on  ’*',j  left  to  white  on  the  right.  It  is  invoked 
by  pressing  the  "G"  key  on  tt..  keyboard. 

DISPLAY  TEST  PATTERN: 


This  function,  when  invoked  by  pressing  the  "P"  key,  displays  the  test  pattern 
in  the  upper  left  quadrant  of  the  monitor.  The  terminal  then  responds  with 
an  asterisk  prompt.  This  prompt  means  that  the  program  is  still  in  the  "test 
pattern"  function  and  there  are  remaining  options  available.  The  response  to 
the  prompt  can  be  "I",  "E",  "P",  or  "carriage  return",  where: 


I 

E 

P 

<cr  > 


indicates  initialize  test  pattern 
indicates  edit  test  pattern 
indicates  enter  primitives 
indicates  return  to  SELECT  FROM  LIST  level 
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The  "I"  and  <cr>  need  no  additional  input.  The  "E"  option  prints  the  terminal 
response  "ELEMENT/SHAQE?".  The  required  input  is  3  hexidecimal  characters. 

The  first  represents  the  column  of  the  element  that  is  to  be  edited;  the  second 
character  represents  the  row.  The  third  character  is  the  shade  of  gray  which 
the  element  will  be  assigned,  black  being  a  value  of  "0"  and  white  a  value  of 
"F".  After  each  entry,  the  response  is  repeated  and  another  element  can  be 
changed;  if  no  more  changes  are  necessary,  a  carriage  return  is  entered. 

The  "P"  option  responds  with  "ENTER  PRIMITIVES".  The  required  response  is  one 
of  the  two  character  mnemonics  corresponding  to  the  PIP  primitives.  These  are 
listed  in  the  description  of  the  INTERPRET  AND  EXE  function.  When  the  list 
of  primitives  have  been  entered,  a  carriage  return  executes  the  primitives  and 
displays  the  first  256  elements  of  PIP's  write  portion  of  the  I/O  buffer  in 
the  test  pattern  format.  The  first  time  this  option  is  selected,  the  result 
is  placed  in  the  upper  right  quadrant  of  the  monitor,  the  next  time  in  the  lower 
left,  and  then  the  lower  right.  The  fourth  time,  the  result  will  overwrite 
the  upper  right  quadrant  and  the  cycle  repeats.  Each  time  the  program  comes 
back  to  the  prompt  level  and  another  “P"  must  be  entered. 

DIGITIZE  ONE  FRAME: 

This  function  is  invoked  by  pressing  the  "D"  key. 

The  program  responds  with  "SELECT  SOURCE(S)  (0-3)".  The  correct  response  is 
any  combination  of  the  numbers  0,  1,  2  or  3.  The  numbers  selected  correspond 
to  the  four  video  sources  which  are  connected  to  the  video  memory.  When  a  source 
is  selected  it  will  be  digitized  once  every  50  milliseconds  for  two  seconds. 

Then  the  next  video  source  selected  will  be  digitized  in  a  similar  manner,  and 
so  on.  To  freeze  the  frame,  any  key  is  pressed.  The  program  then  returns  to 
the  SELECT  FROM  LIST  level. 

TRACK  GATE: 

This  function  is  invoked  by  pressing  the  "T"  key.  The  program  responds  with 
"GATE  SIZE?  (1-8)".  The  single  number  1  through  8  is  then  entered.  The  number 
entered  determines  th?  size  of  a  side  of  the  gate  which  is  drawn  (1  representing 
16  pixels  and  8  representing  8*16  or  128  pixels).  The  gate  can  now  be  moved 
up,  down,  left  or  right  by  pressing  the  ”U",  "D" ,  "L" ,  or  "R"  keys.  When  the 
gate  is  positioned  in  the  desired  position,  the  "H"  key  is  pressed  to  place 
the  gate  and  its  current  contents  back  to  the  upper  left  of  the  monitor.  The 
rest  of  the  monitor  is  set  to  a  black  level. 

The  asterisk  prompt  is  now  printed  indicating  the  program  is  still  in  the  "track 
gate"  function  and  there  are  remaining  options  available.  The  response  to  the 
prompt  can  be  "P",  "S",  "A",  "E",  "T",  "B",  or  "carriage  return",  where: 

P  indicates  enter  primitives 

S  indicates  subtract  gates 

A  Indicates  average  gates 

E  indicates  edge  detect 

T  indicates  threshold 


MICROCOPY  RESOLUTION  TEST  CHART 
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B  Indicates  background  average 

<cr>  Indicates  return  to  SELECT  FROM  LIST  level 

The  program  prints  "INPUT  GATE(S),  OUTPUT  GATE?"  The  valid  response  is  two 
or  three  numbers  each  between  0  and  3  which  represent  quadrants  of  the  monitor 
where  the  input  is  to  be  taken  and  the  output  is  to  appear. 

The  "P"  option  responds  with  "ENTER  PRIMITIVES".  The  required  response  is  one 
of  the  two  character  mnemonics  corresponding  to  the  PIP  primitives.  These  are 
listed  in  the  description  of  the  INTERPRET  AND  EXE  function.  When  the  list 
of  primitives  have  been  entered,  a  carriage  return  executes  the  primitives  and 
displays  the  first  256  elements  of  PIP 1 s  write  portion  of  the  I/O  buffer  in 
the  test  pattern  format. 

The  "S"  option  subtracts  the  pixel  values  of  the  first  input  gate  from  the  pixel 
values  of  the  second  input  gate  and  places  the  results  in  the  output  gate. 

The  "A"  option  averages  the  pixel  values  of  the  two  input  gates  and  places  the 
results  in  the  output  gate. 

The  "E"  option  responds  with  "THRESHOLD?".  The  valid  response  is  a  single  hexa¬ 
decimal  number  0  (zero)  thru  F  (0  *  dark,  F  =  light).  The  program  then  uses 
the  Robert’s  cross  algorithm  to  perform  an  edge  detect  function  on  the  input 
gate  and  places  the  results  in  the  output  gate. 

The  "T”  option  requests  a  threshold  type  and  a  threshold  value.  This  option 
first  responds  with  "THRESHOLD  TYPE  (D,L,B)?'\  The  valid  responses  are  "D", 

"L",  "B"  or  carriage  return  where: 

D  -  (Dark  Threshold)  will  change  all  pixels  with  a  value  greater  than 
the  threshold  value  to  white 

L  -  (Light  Threshold)  will  change  all  pixels  with  a  value  less  than  the 
threshold  value  to  black 

B  -  (Binary  Threshold)  will  change  all  pixels  with  a  value  less  than  or 
equal  to  the  threshold  value  to  black  and  those  with  a  value  greater 
than  the  threshold  value  to  white 

<cr>  -  (carriage  return)  will  perform  no  action  and  allow  the  user  to  select 
another  TRACK  GATE  option. 

The  option  then  prints  "THRESHOLD  VALUE  (0-F)?".  The  valid  response  is  a  single 
hexadecimal  digit  0  (zero)  thru  F. 

The  "B"  option  requests  a  matrix  size.  The  valid  response  is  the  odd  number 
3,  5,  or  7. 


j 
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TG  MEM.  TEST: 

This  function  is  invoked  by  pressing  the  "M"  key.  The  timing  generator  memory 
test  is  executed.  If  no  errors  are  detected,  no  response  Is  given  and  the  pro¬ 
gram  returns  to  the  "SELECT  FROM  LIST"  level.  If  an  error  is  found,  the  bad 
memory  address  and  its  contents  are  printed.  Successful  completion  of  the 
test  requires  about  7  seconds  execution  time. 

10  MEM.  TEST: 

This  function  is  invoked  by  pressing  the  *0"  key.  The  10  buffer  memory  test 
is  executed.  The  results  are  similar  to  the  TG  MEM.  TEST  described  above. 

TG  UTILITIES: 

This  function  is  invoked  by  pressing  the  "U"  key.  The  program  responds  with 
an  asterisk  prompt. 

This  utility  package  allows  the  user  to  examine  and/or  modify  the  contents  of 
the  timing  generator  and  I/O  buffer  memories,  as  well  as  debug  timing  generator 
microcode  using  the  utility's  single  step  capability. 

Load  Timing  Generator  Memory  --  Syntax: 

L  i>  <start  addr>Q  <end  addr>]  <cr> 

Tijis  command  loads  timing  generator  memory  with  microinstructions  from 
SBC-86/12.  <start  a ddr>  and  <end  add r>  specify  the  range  of  addresses 
of  locations  in  timing  generator  memory  whose  contents  are  to  be  loaded 
(if  <end  addr>  is  omitted,  then  <end  addr>  *  <start  addr>).  Both  <start 
addr>  and  <end  addr>  must  be  in  the  range  0  to  7FF.  Oata  for  microword 
N  are  taken  from  bytes  at  addresses  K+6*N  through  K+6*N+5,  where  K  is  a 
constant  whose  value  is  determined  by  the  location  of  the  label 
"START_0F_TG_8UFFER" . 

Examples: 

L  0,1FF  loads  and  verifies  timing  generator  address  0  through  IFF 
from  SBC 

L  200  loads  address  200  from  SBC 
Verify  Timing  Generator  Memory  —  Syntax 

V  p  start  addr  ,  end  addr  cr 

This  command  verifies  the  contents  of  timing  generator  memory  against  SBC- 
86/12.  The  parameters  "<start  addr>"  and  "<end  addr>"  have  the  same  in¬ 
terpretations  as  In  the  load  command.  If  a  discrepancy  Is  found  during 
verification,  a  string  of  6  hexidecimal  bytes  is  printed  out  for  the  ad¬ 
dress  found  to  be  in  error. 


■  i 

I 

i 

i 

j- 
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Example: 


V  5,F  verifies  timing  generator  address  5  through  P  against  SBC 

Save  Timing  Generator  Memory  —  Syntax: 

SB  <start  addr>  [,  <end  addr>]  <cr> 

This  command  stores  the  contents  of  timing  generator  memory  in  SBC  RAM. 

Data  are  taken  from  microwords  at  locations  <start  addr>  through  <end  addr> 
{if  <end  addr>  is  omitted,  then  <end  addr>  *  <start  addr>).  The  6  bytes 
from  microword  N  are  written  to  the  bytes  at  RAM  locations  6*N  through 
6*N+5.  Typically,  timing  generator  microinstructions  are  written  to  disk¬ 
ette  using  the  ICE-86  emulator  or  the  SBC  monitor  upload  feature  once  the 
microinstructions  have  been  transferred  to  SBC  RAM  using  this  command. 

Example: 

S  0,7FF  saves  entire  contents  of  timing  generator  memory 

Display  Microinstructions  —  Syntax: 

DMB  <p*art  addr>[,  <end  addr>]  <cr> 

This  command  displays  timing  generator  microinstructions  as  strings  of 
ASCII  hexadecimal  bytes  on  the  serial  terminal.  The  address  and  6  data 
bytes  corresponding  to  each  microinstruction  between  locations  <start  addr> 
and  <end  addr>  are  printed  out  as  ASCII  hexadecimal  constants  (if  <end 
addr>  is  omitted,  then  < end  addr>  *  <start  addr>). 

Example: 

DM  100  displays  microinstruction  at  location  100 

Change  Microinstruction  Byte  —  Syntax: 

CMB  <micro  addr>,<byte  addr>,<data  byte><cr> 

This  command  changes  a  single  byte  in  a  timing  generator  microinstruction 
and  displays  the  microinstruction  once  the  change  has  been  made.  Byte 
<byte  addr>  of  the  microinstruction  at  location  <micro  addr?  is  set  to 
<data  byte>.  Note  that  microinstruction  bytes  are  numbered  0  through  5, 
with  0  corresponding  to  the  most  significant  byte. 

Example: 

CM  2F,0,AA  sets  byte  0  of  microinstruction  2F  to  AA 
Display  I/O  Buffer  —  Syntax: 

DBB  <start  addr>[,  <end  addr>3  <cr> 
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This  command  displays  the  contents  of  the  PIP  I/O  buffer  between  addresses 
<start  addr>  and  <end  addr>  (if  <end  addr>  Is  omitted,  then  <end  addr> 

»  <start  addr>).  <start  addr>  and  < end  addr>  must  be  legal  I/O  buffer 
addresses  (i.e.,  between  0  and  IFFF). 

Example: 

OB  0,10  displays  I/O  buffer  locations  0  through  10 

Change  I/O  8uffer  —  Syntax: 

CBSJ  <start  addr>  ,  <nbytes>  ,  <data  bytel>  ....  <count>  * 

<data  byte.J>  .....  <data  byte  K>  <cr> 

This  command  sets  one  or  more  bytes  of  the  I/O  buffer  to  specified  values; 
<nbytes>  bytes  starting  with  location  <start  addr>  are  set  to  the  ASCII 
hex  values  specified  in  the  data  list  following  the  <nbytes>  parameter. 

Each  entry  in  this  data  list  consists  of  either  a  single  byte  or  a  repli¬ 
cated  byte  of  the  form  <count>  *  <data  byte>,  the  latter  being  equivalent 
to  <data  byte>  appearing  as  a  single  entry  <count>  times.  The  total  number 
of  bytes  specified  in  the  data  list  must  be<nbytes>. 

Examples: 

CB  0,10,0,B,E*FF  sets  byte  0  to  0,  byte  1  to  B,  and  bytes  2  through 

F  to  FF 

OB  0,100,100*0  clears  bytes  0  through  FF 
Change  Macro  Address  Register  —  Syntax: 

CAfcKmacro  addr><cr> 

This  comnand  sets  the  macro  address  register  on  the  timing  generator  to 
a  specified  address  (the  address  must  be  a  legal  timing  generator  memory 
address). 

Example: 

CA  300  sets  macro  address  register  to  300 

Run  Timing  Generator  —  Syntax: 

R#  <clock  rateXcr> 

This  command  enables  the  timing  generator  sequencer  to  run  continuously 
at  a  specified  clock  rate.  Clock  rate  must  be  in  the  range  0  through 
7  and  selects  the  clock  rate  as  follows: 

0  »  8.33  MHz 

1  *  4.16  MHz 

2  *  2.08  MHz 

3  *  1.04  MHz 
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4  »  500  KHz 

5  «  250  KHz 

6  «  125  KHz 

7  *  63  KHz 


Example: 

R  0  runs  timing  generator  at  8.33  MHz 

Halt  Timing  Generator  —  Syntax: 

H<cr> 

This  command  disables  the  timing  generator  sequencer  from  running  continu¬ 
ously. 

Single  Step  Timing  Generator  —  Syntax: 

G  [[t]0[<start  addr>,]<end  addr>j<cr> 

This  command  single  steps  the  timing  generator  sequencer  either  uncondi¬ 
tionally  or  until  a  specified  microinstruction  address  is  encountered. 

In  its  simplest  form,  the  command  steps  the  sequencer  one  instruction; 
this  is  the  form  "G<cr>."  The  form  "G[T]  0[<start  addr>,]<end  addr><cr>" 
specifies  that  the  sequencer  is  to  be  stepped  until  the  microinstruction 
at  location  <end  addr>  has  executed.  If  <start  addr>  is  present,  then 
execution  is  forced  to  start  at  address  <start  addr>;  otherwise,  execution 
starts  at  whatever  address  the  sequencer  is  currently  generating.  The 
*T"  parameter,  when  present,  specifies  that  a  trace  be  built  of  those  ad¬ 
dresses  generated  by  the  sequencer  until  <end  addr>  is  encountered.  Each 
address  generated  by  the  sequencer  after  each  single  step  is  written  to 
the  trace  buffer  starting  with  buffer  location  0.  The  buffer  holds  up 
to  256  addresses;  if  the  buffer  overflows,  single  step  execution  is  halted 
and  a  message  to  that  effect  is  printed  on  the  terminal.  Once  single  step¬ 
ping  has  halted,  either  normally  or  abnormally,  a  word  called  the  trace 
pointer  gives  th?  index  of  the  last  location  written  in  the  trace  buffer. 
The  contents  of  both  the  trace  pointer  and  the  trace  buffer  can  be  examined 
using  the  commands  in  the  Display  Trace  Pointer/Trace  Buffer  command. 

Examples: 

G  steps  once 

G  200  steps  until  address  200  is  executed 

G  300,200  starts  execution  at  address  300  and  steps  until  address 

200  is  executed 

GT  2FF  enables  trace  and  steps  until  address  2FF  is  executed 
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Display  Trace  Pointer/Trace  Buffer  —  Syntax: 

DTP<er>  display  trace  pointer 

DT  <start  addr>Q<end  addr>]  <cr>  display  trace  buffer 

These  commands  allow  examination  of  the  trace  buffer  as  described  in  the 
Single  Step  Timing  Generator  command. 

The  first  command  listed  above  displays  the  contents  of  the  trace  pointer 
so  that  the  user  can  determine  how  much  of  the  trace  buffer  contains  useful 
Information. 

The  second  command  displays  the  contents  of  the  trace  buffer.  In  its  sim¬ 
plest  form,  the  latter  command  is  "DT<cr>,"  which  displays  the  entire  trace 
buffer.  If  <$tart  addr>  and  <end  addr>  are  present,  the  contents  of  the 
trace  buffer  between  locations  <start  addr>  and  <end  addr>  are  displayed 
(if  <end  addr>  is  omitted,  then  <end  addr>  »  <start  addr>). 

Examples: 

OT  0,1FF  displays  locations  0  through  IFF  of  trace  buffer 

DT  50  displays  location  50  of  trace  buffer 

Initialize  Timing  Generator  —  Syntax: 

I<cr> 

This  command  issues  a  reset  pulse  to  the  timing  generator  and  forces  the 
sequencer  to  execute  the  microinstruction  at  address  0. 

Quit  —  Syntax: 

<cr> 

A  carriage  return  on  a  line  by  itself  returns  the  program  to  the  SELECT 
FROM  LIST  level. 

Certain  control  characters  have  special  functions  when  entered  by  the  user  from 
the  serial  terminal.  Control  S,  when  entered  during  any  data  printout  on  the 
terminal,  causes  the  printout  to  temporarily  halt;  entering  control  Q  causes 
the  printout  to  resume.  Control  X,  when  entered  as  part  of  a  command  line  being 
typed  in  by  the  user,  causes  that  line  to  be  deleted.  Control  C  has  the  effect 
of  returning  control  to  the  SBC  monitor  program. 

INTERPRET  AND  EXE: 

This  function  is  invoked  by  pressing  the  "I"  key.  The  program  responds  with 
“ENTER  PRIMITIVES".  The  correct  response  is  one  of  following  2  character  mne- 
Inonics: 
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SR  Shift  Array  Right 

SL  Shift  Array  Left 

SU  Shift  Array  Up 

SD  Shift  Array  Down 

LI  Load  Store  #1  from  X 

L2  Load  Store  #2  from  X 

L3  Load  Store  #3  from  X 

U1  Unload  Store  #1  to  X 

U2  Unload  Store  #2  to  X 

U3  Unload  Store  #3  to  X 

El  Erase  Store  #1 

E2  Erase  Store  #2 

E3  Erase  Store  #3 

M2  Store  #1  minus  Store  #2 

R1  Regenerate  #1 

R2  Regenerate  #2 

OS  Execute  once 

CO  Execute  continuously 

IA  Input  entire  array 

OA  Output  entire  array 

D3  Store  #1  minus  Store  #3 

D4  Store  #1  minus  X 

RF  Rotate  X,  Store  #2,  Store  #3 

SW  Swap  Store  #2,  #3 

RB  Rotate  X,  Store  #3,  Store  #2 

MV  Move  Store  #1  to  Store  #2 

MX  Move  Y  to  X 

After  each  mnemonic  is  entered,  the  code  is  interpreted  and  a  "JUMP  TO  SUBROUTINE" 
2910  microinstruction  is  placed  in  the  next  available  TG  memory  location.  When 
all  primitives  have  been  entered,  a  carriage  return  is  typed  as  the  first  charac¬ 
ter  of  a  line.  This  causes  the  TG  sequencer  to  execute  the  primitives  and  re¬ 
turns  the  program  to  the  "SELECT  FROM  LIST"  level. 

LIST  FRAME  MEMORY: 

This  option  is  not  complete.  Any  response  from  the  user  will  cause  a  return 
to  the  SELECT  FROM  LIST  level.  This  option  is  intended  to  allow  the  user  to 
list  the  contents  of  frame  memory  to  a  printer. 
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